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During the tenure of this grant, we have been looking into the chemistry of the 
nightside ionosphere of Venus with a view toward elucidating the relative roles of 
electron precipitation and plasma transport as sources of the nightside ionosphere. 
Secondary goals have included determining the densities of minor species on the 
nightside, and verifying the relative normalization of the Pioneer Venus orbiter 
ion mass spectrometer (OIMS) and orbiter neutral mass spectrometer (ONMS) in 
the photochemical equilibrium region. Our studies have involved a combination of 
numerical modeling and analysis of the Pioneer Venus UADS data base, specifically 
data from the OIMS, ONMS and electron temperature probe (OETP). We have set 
up a one-dimensional model of the Venus nightside ionosphere, in which downward 
fluxes of atomic ions are introduced at the upper boundary to simulate transport 
of ions from the dayside. Our model shows that the densities of mass-28 ions (C0 + 
+ N^) resulting from an influx of atomic ions from the dayside are quite small, 
due to the high ionization potentials of CO and N 2 that make chemical production 
difficult. A look at the data reveals that the actual densities of mass-28 ions are 
quite variable, from values near 10 to more than 10 4 cm -3 . We suggest the excess 
mass-28 ions are produced by electron precipitation and that the presence of high 
densities of mass-28 ions is a signature of “auroral” precipitation. Details are given 
in the attached reprint of an article published in GRL in September, 1990 by H. 
A. Taylor and myself entitled “A signature of auroral precipitation in the nightside 
ionosphere of Venus”. 

Also attached are two other articles supported partially under this grant. The 
first, a chapter of the book Dissociative recombination: experiment, theory and 
applications is entitled “Dissociative recombination in aeronoiny”. The chapter 
contains a discussion of the atomic oxygen green line in the nightglow of Venus, 
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which is produced mainly by dissociative recombination of 0^ . In this work, the 
0+ is assumed to be produced by reaction of 0 + transported from the dayside with 
CO 2 . In fact, the presence of electron precipitation will increase the emission in the 
green line in two ways: by increasing the density of 0^ and by direct production 
of 0( 1 JD) from electron impact on atomic oxygen. 

The second is a review article that Steve Bougher and I have written entitled 
“Structure, luminosity and dynamics of the Venus thermosphere”, to be published 
in Space Science Reviews soon. This manuscript, while officially a review article, 
contains original calculations of production rates of excited states for models based 
on Pioneer Venus data, that I needed to do for completeness. Tables 2, 3, 5, 6, 7, 
and 11 contain original, unpublished calculations. 
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A SIGNATURE OF AURORAL PRECIPITATION IN THE NIGHTSIDE IONOSPHERE OF VENUS 

J. L. Fox 

Institute for Terrestrial and Planetary Atmospheres 
State University of New York at Stony Brook 

H. A. Taylor, Jr. 

Taylor Enterprises, Beaver River, Nova Scotia 


Abstract. We show here that the densities of mass- 
28 ions measured by the Pioneer Venus Orbiter ion mass 
spectrometer (OIMS) on the nightside of Venus are highly 
variable and show little correlation with the values of the 
0+ densities. We have determined the total production 
rates of mass-28 ions in the chemical equilibrium region 
and find that this production rate cannot be explained 
by known chemical production reactions. We propose 
that the “excess” production is due to precipitation of 
electrons into the nightside thermosphere. 

Introduction 

The source of the Venus nightside ionosphere has been 
disputed since the first radio occultation measurements of 
the nightside electron density profile by Mariner 5 [Kliore 
et al., 1967]. Since then, the viable mechanisms have been 
reduced to two: transport of atomic ions, mostly O , 
from the dayside and precipitation of energetic electrons 
that have been observed in the Venus umbra. Models ap- 
propriate to high solar activity have shown that plasma 
transport is the dominant maintenance mechanism, since 
it alone can produce the observed densities of , with 
electron precipitation playing a possible secondary role 
in enhancing the lower 0^ peak [Spenner et al., 1981, 
Cravens et al., 1983]. Gringauz et al. [1977] maintain 
that precipitation of electrons alone is sufficient to pro- 
duce the main ion peak, as measured by the Venera 9 and 
10 radio occultation experiments at solar minimum. 

We have been investigating the chemistry of the night- 
side ionosphere, with a view toward elucidating the rel- 
ative roles of electron precipitation and plasma trans- 
port in its maintenance. Model calculations based on 
ion transport as the source of the nightside ionosphere 
showed that the reaction 

0 2 + N — ► NO + + O (1) 

was the major source of NO+in the chemical equilibrium 
region. Analysis of data from eight Pioneer Venus (PV) 
orbits chosen for low periapsis and the availability of data 
showed that the NO + densities could be modeled with re- 
action (1) as the sole source for some orbits, but for others 
it was completely inadequate. A close look at those orbits 
showed that they were all characterized by high densities 
of mass-28 ions (N t + CO+). This deviation is easily ex- 
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plained since Nj reacts rapidly with O to produce NO+ 

Nj+0-NO + + N. (2) 

Because of the high ionization potentials of both N 2 and 
CO, the chemical sources are limited. Yet a survey of 
early Pioneer Venus orbits with periapsis below 165 km 
reveals a large variability in the maximum densities of 
mass-28 ions with values ranging from about 10 cm to 
about 10 4 cm -3 . Taylor et al. [1982] had earlier remarked 
that the occasional large densities of NO + and 0 2 mea- 
sured by the OIMS appeared to require a source in ad- 
dition to ion transport; they suggested that the densities 
were enhanced by particle precipitation, but the mech- 
anism of the source was not resolved. We propose that 
the excess NO+ is produced from Nf created by parti- 
cle precipitation, and that high densities of mass-28 ions 
in the nightside ionosphere are an in situ signature of 
precipitation of electrons into the nightside atmosphere. 

In the discussion that follows, we identify the sources 
and sinks of mass- 28 ions in the chemical equilibrium re- 
gion, that below approximately 165-170 km. There the 
ion production and loss rates are equal, so the magnitude 
of the source can be estimated because the total loss rate 
can be computed fairly accurately. From the loss rate 
and the known chemical sources, the “excess” produc- 
tion rate, that we ascribe here to particle precipitation, 
is determined. 

Sources of Nj and CO + 

Reactions, references and rate coefficients used in this 
work are summarized in Table 1, and the reactions num- 
bers cited below refer to the numbers in that table. Be- 
cause of the high ionization potential of N 2 (15.58 eV), 
only a few charge transfer reactions of atmospheric spe- 
cies are energetically capable of producing N 2 , including 
charge transfer from He+, 0++ and O + ( 2 D) (reactions 3- 
5). Reaction (3) is potentially important, since nighttime 
He + densities are larger than daytime values, peaking in 
the the pre-dawn sector [Taylor et al., 1980]. Densities 
of 0 + + were measured by the PV OIMS, but are often 
above the noise level of the instrument only above the 
chemical equilibrium region. In addition, the products of 
reaction (4) are unknown; the large exothermicity sug- 
gests that the dissociative channel probably dominates. 
Charge transfer from 0+( 2 D) to N 2 , reaction (5), is a ma- 
jor source of Nj in the dayside ionosphere [Fox, 1982], 
and the importance of this reaction depends on whether 
the metastable 0+( 2 D) ion could survive transport from 
the dayside. A detailed analysis would be required to 
model the role of reaction (5) in producing the observed 
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Fox and Taylor: Auroral Signature in Venus Ionosphere 
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Fig. 2. Same as for Fig. 1, but for Orbit 529 (inbound). 
In Fig. 2b, the total production rate and excess produc- 
tion rate are essentially superimposed. 

cm' 3 s _1 . For orbit 505, the maximum auroral produc- 
tion rate is about a factor of 40 less than for orbit 529, 
although the production rates due to chemical reactions 
are about the same for the two orbits. 

Other evidence for the precipitation of electrons into 
the nightside thermosphere exists. Fluxes of suprather- 
mal electrons have been measured above the atmosphere 
in the Venus umbra by the PV retarding potential an- 
alyzer (ORPA) (e.g. Knudsen and Miller, 1985) as well 
as by the plasma analyzers on the Soviet Venera 9 and 
10 spacecraft [e.g. Gringauz et ah, 1977). Continuous 
but highly variable emissions of atomic oxygen at 1304 
and 1356 A have been observed in images of the night- 
side of Venus by the PV orbiter ultraviolet spectrometer 
(OUVS) [Phillips et ah, 1986). Fox and Stewart [1990] 
have proposed that these emissions are caused by the pre- 
cipitation of soft electrons into the nightside atmosphere. 
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titude below 165 km, and data at the mass-28 ion peak 
for the densities of 0 + , O^, He + , CO 2 , O, and N 2 - In 
this calculation, only the reactions of He + (reactions 3 
and 6) were included as chemical sources of mass-28 ions, 
since experience with individual orbits showed the other 
sources to be unimportant. The electron temperature 
was assumed to be about 1000 K for all the data points. 
Errors resulting from this assumption should be small, 
since the electron temperature dependence of the disso- 
ciative recombination coefficient is weak, and dissociative 
recombination is not the major loss mechanism. In order 
to reduce the effects of variability due to changes in the 
neutral atmosphere and in the altitude of periapsis, the 
excess production “frequency” was obtained by dividing 
the rate by the factor 3 [N 2 ) + 0.1[CO2). The first term 
in the factor represents ionization of N 2 and CO, and 
was chosen because the CO density is about twice the N 2 
density on the nightside above the homopause [Hedin et 
al., 1983]. The second term represents the contribution 
of dissociative ionization of CO 2 ; the factor 0.10 is the 
high energy ratio of the electron impact cross sections 
for production of CO + from CO 2 and CO. On the day- 
side, the excess production is due to direct ionization and 
chemical reaction other than reactions (3) and (6). The 
production frequency appears large on the dayside also 
because CO is more abundant relative to N 2 there [Hedin 
et ah, 1983). The variability on the dayside is small, be- 
cause all the data were taken near solar maximum, but 
some scatter remains from the variation of periapsis alti- 
tude. The large variability in the source on the nightside, 
ranging over more than two orders of magnitude, itself 
suggests particle precipitation as the source mechanism. 
Some excess production appears to be present in nearly 
all the orbits. This is consistent with the morphology of 
the 1304 A aurora [Phillips et ah, 1986], which is nearly 
always visible on the nightside, but exhibits large spatial 
and temporal variations. 



LOCAL TIME 

Fig. 4. Ratio of [28 + ][0]/[N 2 ] at the mass-28 ion peak to 
the density of 0 + at its peak, as a function of local time 
for the first 600 orbits. The variability on the nightside is 
substantially greater than the variability on the dayside. 


Fig. 3. Production frequency of mass-28 ions in excess 
of that that can be accounted for by chemical reactions 
at the mass-28 ion peak as a function of local time for 
all orbits of the first 600 for which sufficient data are 
available. 

Figure 3 shows the excess production “frequency” of 
mass-28 ions at the mass- 28 ion peak for all orbits among 
the first 600 that meet the following criteria: periapsis al- 


Since models have shown that the measured 0 + densi- 
ties can only be accounted for by transport from the day- 
side, and if mass-28 ions are mostly produced by auroral 
processes, the relative values of the maximum mass-28 
ion and 0 + densities would be indicative of the relative 
importance of ion transport and particle precipitation in 
the nightside ionosphere. Figure 4 is a plot of the ratio of 
the peak mass-28 ion density (divided by [N 2 ]/[0]) to the 
peak 0 + density as a function of local time for all orbits 
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Abstract. We review here observations and models related to the chemical and thermal structures, airglow 
and auroral emissions and dynamics of the Venus thermosphere, and compare empirical models of the 
neutral densities based in large part on in situ measurements obtained by the Pioneer Venus spacecraft 
Observations of the intensities of emissions are important as a diagnostic tool for understanding the 
chemical and physical processes taking place in the Venus thermosphere. Measurements, ground-based and 
from rockets, satellites, and spacecraft, and model predictions of atomic, molecular and ionic emissions, 
are presented and the most important sources are elucidated. Coronas of hot hydrogen and hot oxygen have 
been observed to surround the terrestrial planets. We discuss the observations of and production 
mechanisms for the extended exospheres and models for the escape of lighter species from the atmosphere. 
Over the last decade and a half, models have attempted to explain the unexpectedly cold temperatures in 
the Venus thermosphere; recently considerable progress has been made, although some controversies 
remain. We review the history of these models and discuss the heating and cooling mechanisms that are 
presently considered to be the most important in determining the thermal structure. Finally, we discuss 
major aspects of the circulation and dynamics of the thermosphere: the sub-solar to anti-solar circulation, 
superrotation, and turbulent processes. 


1. Introduction 

C0 2 was first identified as a constituent of the Venus atmosphere by Adams and 
Dunham (1932), who, in their search for evidence of 0 2 and H 2 0 in ground-based 
infrared spectra observed these C0 2 bands at 7820, 7883, and 8689 A. Over the 
following forty years, spectral resolution in the infrared improved and allowed the 
abundances of HCl, HF, and CO in the atmosphere to be determined (e.g., Connes 
et ai, 1967, 1968, 1969). A review of this early period in Venus spectroscopy has been 
presented by Young (1972) and more recently by von Zahn etai (1983). That the 
atmosphere is composed mostly of C0 2 was suggested by Connes et ai (1967), Belton 
et ai (1968), and by Moroz (1968), but confirmation was not obtained in the first in situ 
measurements were performed in 1967 by instruments on the Soviet Venera 4 entry 
probe (Vinogradov et ai , 1968). The experiments were simple: the pressure of a sample 
of atmospheric gas was determined before and after adsorption by a chemical absorber. 
A description of the instruments on Venera 4, and on two similar spacecrafts that were 
launched two years later, Veneras 5 and 6, have been presented by Vinogradov (1971). 
The major constituents and their mixing ratios (used hereafter to mean fraction of 
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particles by number) in the lower atmosphere were reported by Kuz'min (1970, 1971) 
as 95 ± 2% C0 2 , 2-5% N 2# 0.1-1% H 2 0 and less than 0.4% 0 2 . 

Up to the homopause, the composition of the thermosphere of a planet is nearly the 
same as that of the bulk atmosphere, at least with respect to the major constituents that 
are not affected by photochemistry ; above the homopause photolysis of the molecular 
constituents and diffusive separation in the gravitational field causes lighter constituents 
to become increasingly abundant. Exploration of the thermosphere of Venus began in 
October 1967 when two spacecrafts, the U.S. fluby Mariner 5 and the Soviet Venera 4 
carried ultraviolet photometers to Venus. Both instruments were designed to detect 
resonance emission features of H at 1216 A and of O at 1304 A, from which column 
densities and altitude profiles of the H and O densities could be derived. No signal was 
reported at 1304 A by either teh Mariner 5 photometer (Barth etai , 1967) or the 
Venera 4 instrument, which observed the atmosphere in darkness below about 350 km 
(Kurt etal. t 1968). Both detected strong La emission. The Venera La intensities were 
used to derive a nightside density profile of H atoms from 1 to about 4 Venus radii. From 
the scale height of the H densities and the upper limit on the O densities Kurt et al. ( 1 968) 
deduced, more or less correctly, that ‘the Venus nightside atmosphere is cold . . . and 
undergoes a sharp transition to interplanetary space'. The Mariner 5 La channel 
detected emission that was characterized by an abrupt increase in scale height as the 
minimum distance of the photometer line-of-sight to the center of the planet exceeded 
about 9000 km (Barth et al ., 1967; Barth, 1968). From the scale height, the ratio of the 
temperature to the mass of the species responsible for the of the inner component was 
derived to be about 350 Kamu~ I . Exospheric temperatures as low as 350 K were 
considered unlikely, based partly on model calculations (McElroy, 1969) that were 
reinforced by prejudices arising from the analogy with Earth, where the solar fluxes are 
a factor of two lower and the average exospheric temperature is about 1000 K. There- 
fore, species heavier than H, including D and H 2 were considered also as the source 
of the low altitude component, the latter assumed to produce emission from photodisso- 
ciative excitation rather than resonance scattering (e.g., Barth, 1968; Donahue, 1968; 
Wallace, 1969; McElroy and Hunten, 1969). Later work showed that the actual 
exospheric temperature was even lower than the ‘unlikely' value of 350 K (Anderson, 
1975). The measurement of the outer, larger scale height component was the first 
detection of the hot hydrogen corona, although it was not immediately identified as such. 

That the thermospheric abundance of O was low was suggested by the non-detection 
of 1 304 A emission and by profiles of the electron density in the daytime ionosphere that 
were obtained from the radio occultation experiment on Mariner 5 (Kliore et a/., 1967). 
The electron density profile, shown in Figure 1, exhibited a single maximum that was 
interpreted as an F, peak, a Chapman layer consisting of a molecular in (initially believed 
to be C0 2 + ) produced by photoionization and photoelectron impact ionization and 
destroyed by dissociative recombination. The absence of a distinguishable O * (F 2 ) peak 
above the F, layer was interpreted as implying that the mixing ratio of atomic O in the 
Veznus thermosphere was smallk (Stewart, 1968; McElroy, 1969; Kumar and Hunten, 
1974). 
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Fig. 1 . Electron density profile as measured by the radio occultation experiment on Mariner 5. Taken from 

Kliore etai (1967), 

Mariner 10 flew by Venus in February 1974, carrying an objective grating spectrome- 
ter with channel electron multipliers at the wavelengths of 9 atomic and molecular 
emission features between 200 and 1700 A, two zero-order channels and a background 
channel at 430 A (Broadfoot era/., 1977). The spectral features that the Mariner 10 
instrument was designed to detect are listed in Table I. Strong emissions were detected 
in the positions of the helium, hydrogen, oxygen, CO, and carbon features in the 
spectrum (Broadfoot etai , 1974). 

In October 1975, the Soviet spacecraft Veneras 9 and 10 were inserted into orbit about 
Venus. The spacecraft were equipped with hydrogen and deuterium absorption cells in 
front of an NO counter to study the hydrogen corona, and visible spectrometers to 
record nightglow emission features in the 3000 to 8000 A wavelength range (Keldysh, 
1977; Krasnopol’sky, 1983). Four band systems of 0 2 were detected in the Venus 
nightglow with maximum volume emission rates appearing near 100 km. Upper limits 
were placed on other species and processes from the absence of emission features at 
5577 A (the atomic oxygen green line), 5893 A (the sodium D-line doublet) and 3914 A 
(the (0,0) N 2 * first negative band) (Krasnopol’sky, 1978, 1979). Unfortunately, the 
sensitivity of the detectors decreased by a factor of 5 after the first month of operations, 
so the number of orbits yielding useful data was limited. 

In December 1978 ten separate space probes arrived at Venus, including the Pioneer 
Venus (PV) Orbiter and Multiprobe, which included the Bus and four probes that were 
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TABLE I 

Features detected by the photometers on Manner 10 and Veneras 11 and 12 


Species 

Wavelength (A) 

Mariner 10 

Venera 1 1 

Venera 12 

He * 

304 

no 

yes 

yes 

He 

584 

yes 

yes 

yes 

Ne 

740 

no 

no 

no 

O* 

834 

n/a* 

yes 

n/a 

At 

869 

no 

n/a 

no 

At 

1048 

no 

no 

no 

H 

1216 

yes 

yes 

yes 

O 

1304 

yes 

yes 

yes 

C 

1657 

yes 

yes 

yes 

CO 

1480-1500 

yes 

yes 

yes 


• No detector at that wavelength. 


released prior to entry into the Venus atmosphere, and the Soviet Veneras 11 and 12 
flyby spacecraft and landers. The composition of the atmosphere below the clouds was 
measured by a mass spectrometer (LNMS) (e.g., Hoffman etai , 1980a) and gas 
chromatograph (LGC) (e.g., Oyama et al 1980) on the Pioneer Venus sounder probe, 
neutral mass spectrometers on the Venera 1 1 and 12 landers (e.g., Istomin ex a/., 1979; 
Grechnevef a/., 1980) and a gas chromatograph on the Venera 12 lander (GePraan era/., 
1979). A comparison of the results of these measurements has been presented by 
Hoffman et al. (1980b). C0 2 and N 2 are found to make up 99.9% of the atmosphere. 
The range of values of the N 2 mixing ratio was 2,5 to 4.5 % . Von Zahn et al. (1983) have 
reviewed the available measurements of N 2 in the bulk atmosphere and recommend an 
abundance of 3.5 + 0.8% ; that value has also been adopted by the Venus International 
Reference Atmosphere (von Zahn and Moroz, 1986) (see Section 2.2). 

The Venera 11 and 12 flybys carried UV spectrophotometers with detectors at 9 
wavelengths with bandpasses about 12 A wide and zero-order channel that provided 
an integral of the intensity over the wavelength region from 300 to 1700 A. The wave- 
lengths of and species responsible for the emission features that the instruments were 
designed to detect are also listed in Table I. The Venera instruments were similar to, 
but more sensitive than, the Mariner 10 instrument (Bertaux et a/., 1981 ; KrasnopoPsky, 
1986a). Unlike the Mariner 10 instrument, however, during the flyby the slits were not 
oriented parallel to the limb and consequently the spatial resolution was insufficient to 
derive altitude profiles of the intensities; except for La and the He 584 A emission, only 
surface brightnesses were reported (Kurt etai , 1979, 1983; Bertaux etai , 1981). The 
resonance emissions of Ar and Ne were not detected, but the first observations of the 
O* and He + emission features at 834 and 304 A, respectively, were reported. 

The Pioneer Venus bus and orbiter carried a number of instruments that have 
provided information about the structure and composition of the thermosphere. A 
magnetic sector neutral mass spectrometer (BNMS) aboard the bus measured neutral 
densities from about 650 km to below the homopause, near 128 km (e.g., von Zahn, 
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1977; von Zahn et ai , 1980). The orbiter payload included a quadrupole neutral mass 
spectrometer (ONMS) to measure thermospheric densities of C0 2 , 0, CO, N 2 , N, and 
He (e.g., Niemann et ai y 1979, 1980a, b), and an atmospheric drag experiment (OAD) 
to measure the total mass density profile (e.g., Keating, 1977 ; Keating et ai, 1980). The 
orbiter ultraviolet spectrometer (OUVS) was designed to measure intensities of airglow 
emission features of neutral atmospheric atoms, molecules, and ions. A complete 
description of the OUVS has been presented by Stewart (1980). Briefly, the instrument 
consists of 250-mm fj 5 Cassegrainian telescope, an Ebert-Fasie monochromator of 
125 mm focal length, and two photomultipliers, sensitive in the ranges 1 100 to 1800 A 
and 1600 to 3300 A. The OUVS instrument can be commanded to scan the spectrum, 
or the wavelength may be fixed. Images of the planet at a fixed wavelength may be 
obtained from successive scans of the disk that are made as the spinning spacecraft 
approaches the planet from apoapsis. When the spacecraft is nea^periapsis altitude 
profiles of the emissions may be obtained. Among the features that the OUVS has 
detected are the y and & bands in the NO nightglow, the 1304, 1356, and 2972 A features 
of O, H La, the 1561 and 1657 A features of C, and the CO fourth positive bands in 
the dayglow, and the 1304 and 1356 A auroral emission features. 

We present here a brief history of investigations into and our current understanding 
of the structure, emission features, and dynamics of the Venus thermosphere. In 
Section 2, we present a description of the basic structure of the thermosphere, including 
the exospheric temperature and chemical composition, mostly as determined by Pioneer 
Venus in situ measurements. We compare empirical models of the major thermospheric 
constituents. In aurora, both Earth-bsed and from American and Soviet spacecraft. 
Measurements and models of emission features of the species, O, CO, N 2 , 0 2 , C, NO, 
and He and their ions are described. Discussion of emission features of atomic hydrogen 
is deferred until Section 4, where we describe the hot atom coronas and atmospheric 
escape. Section 5 contains a detailed discussion of the thermal structure and heating and 
cooling mechanisms, mostly as given by models. The last part of the paper, Section 6, 
describes the circulation and dynamics of the thermosphere. 


2. Basic Thermospheric Structure 
2.1. Exospheric temperature 

The prejudice that Venus exospheric temperatures, 7^, should not be too much smaller 
than terrestrial values, persisted until about 15 years ago, in spite of early indications 
that the temperatures were lower. The scale height of the H profile derived from Venera 4 
La intensities close the the planet implied temperatures on the Venus nightside near 
100 K, but a specific value was not presented by Kurt et ai (1968), who merely remarked 
that the nightside atmosphere was cold. Stewart ( 1 968) first suggested the interpretation 
that the two-scale height La profile measured by Mariner 5 was H atoms at two different 
temperatures, which would imply that the temperature of the inner component on the 
dayside was about 350 K. Over the next few years, further measurements and models 


SPAC 191/5 



6 


J. L. FOX AND S. W. BOUGHER 


ruled out H 2 and D at higher temperatures as sources of the emissions (see Section 4 
for further details). The higher temperatures, near 700 K, computed by McElroy (1969) 
could be accommodated only if C0 2 * was considered to be the major ion. Fehsenfeld 
et al. (1970) showed that in the presence of a small admixture of O, C0 2 * is converted 
rapidly to 0 2 * through the reaction 

co; + o^o; + co. d) 

Kumar and Hunten (1974) showed that a model ionosphere with an exospheric tem- 
perature of 350 K and as the major ion fit the Mariner 5 electron density profiles. 
Mariner 10 measurements of the He corona gave a T ^ of about 375 ± 105 K (Kumar 
and Broadfoot, 1975). Anderson (1976) reanalyzed the Mariner 5 La data and reported 
for the inner component a day side temperature of 275 K and a nightside temperature 
of 150 K. These low values were later confirmed by the Pioneer Venus in situ measure- 
ments. A value of T ^ of 275 ± 15 K was derived from the He densities measured by 
the Bus Neutral Mass Spectrometer at 08 : 30 local time, 37.9° S latitude (von Zahn 
etal., 1980). Keating et al. (1980) constructed an empirical model from the OAD data 
obtained in the early PV orbits; day side values of T ^ were about 280-300 K. Airglow 
emission profiles measured by the OUVS were consistent with an exospheric tempera- 
ture of about 275 K (Stewart et al., 1979). Hedin et al. (1983) derived a global average 
temperature of 228 K from ONMS data taken between December 1978 and August 
1980 for a narrow latitude range of a few degrees around 16 c N. Variations with solar 
activity during this limited time period were observed to be only 10% of the change seen 
in the terrestrial thermosphere. The variation in the exospheric temperatures with solar 
flux in the model of Hedin et al. (1983) is 0.14 K per unit /* l07 , the 10.7 cm solar flux 
in units of 10 ~ 22 W m ~ 2 Hz " 1 at 1 AU. 

The temperature at the lower boundary of the thermosphere, near 100 km, takes on 
values near 175 K and, unlike the exospheric temperature, varies little with solar zenith 
angle (Seiff eta/., 1980). The thermospheric temperature structure will be discussed in 
detail in Section 5. 

2.2. Models of the Venus thermosphere 

Information about the composition and structure of the thermosphere of Venus, 
obtained largley from in situ measurements, have been summarized in several empirical 
models. Von Zahn etal. (1980) presented a model of the momingside atmosphere 
(hereafter the BNMS model) that is based mostly on measurements of densities of C0 2 , 
He, N 2 , and CO made by the BNMS. The 10.7 cm solar flux index (F l0 7 ) was about 
189 at the time of the bus entry. The densities of O in the model were derived from the 
ratio 0/C0 2 of 3 reported at 167 km by Niemann et al. (1980b) from ONMS measure- 
ments. The mass densities measured by the BNMS agree fairly well with those measured 
by the OAD and merge smoothly near 130 km with densities measured by the 
accelerometers on the entry probes (Seiff etal., 1980). Altitude profiles of the major 
constituents in the model are shown in Figure 2. The exospheric temperature is 275 K 
and the temperature at 100 km, the lower boundary of the model, is 178 K. 
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Fig. 2. Number densities of He, N 2 , CO, N 2 + CO, and C0 2 measured by the BNMS versus altitude, 
interpolated to integer altitudes. Also shown are mass densities measured by the large probe accelerometer 
from 100 to 126 km, those derived from the BNMS between 130 and 180 km, and those from orbiter drag 
data (OAD) between ISO and 170 km. The lines represent number densities, mass density and temperature 
from the BNMS model. Taken from von Zahn etal. (1980). 


The density profiles of He, CO, and N 2 were used by von Zahn et al. (1980) to derive 
an expression for the eddy diffusion coefficient K (in units of cm J s " 1 ): 


where n is the total number density in cm -23 . The constant A was determined to be 
within a factor of 2 of 1.4 x 10 13 . The molecular diffusion coefficient D has a steeper 
altitude dependence, since it is proportional to n~ l . K and D are equal at the 
homopause, which for N 2 is near 136 km. The molecular diffusion coefficient for He is 
larger than that of N 2 and the He homopause is located lower, near 130 km. Extra- 
polating the mixing ratios to the lower atmosphere, the composition of the lower 
atmosphere is derived to be about 95.5% C0 2 , 4.5% N 2 , and 12 ppm He. The mixing 
ratio of He is very sensitive to the uncertainty in the eddy diffusion coefficient and is 
good to only about a factor of three. The probably uncertainty in the N 2 mixing ratio 
is about 25%. 

Hedin etal. (1983) presented a global empirical model of the Venus thermosphere 
(hereafter referred to as the VTS3 model) based on measurements of C0 2 , O, CO, N 2 , 
He, and N made by the Pioneer Venus ONMS. The VTS3 model incorporates 
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variations in the densities with latitude, local tune and solar activity based on three 
diurnal cycles of observations. The altitude mage of the orbiter measurements is from 
about 142 to 250 km, but data from the entry probe accelerometers (Seiff etai , 1980), 
the BNMS (von Zahn et ai, 1980), and theoretical considerations were used to extend 
the mode! to 100 km. Dissociation of C0 2 is the primary source of O and CO in the 
lower thermosphere; the production is balanced by three-body recombination 

0 + 0 + C0 2 -*0 2 + C0 2 (3) 

below 1 10 km (Massie et ai , 1983). Thus the mixing ratio of CO should be equal, to a 
first approximation to the sum of the mixing ratio of O plus twice that of 0 2 . In the VTS3 
model, the mixing ratios of O and CO are assumed to be equal at 100 km, since 
photochemical considerations indicate that the mixing ratio of 0 2 is small (Yung and 
DeMore, 1982). Data from the PV sounder probes, shown in Figure 3 compared to the 
ONMS model densities, indicate that latitude and local time variations of the mass 
densities near i 00 km are small (Seiff et ai , 1980), so pressure gradients at the lower 
boundary in the model are constrained to be minimal. The mass densities measured by 
the ONMS arc about 60% less than those given by the simpler PV OAD, the 
accelerometers on the probes and the BNMS. Ionospheric models also indicate that the 
ONMS densities are low (e.g., Cravens et ai , 1981). It has been suggested that the gas 
flow into the ONMS ion source at satellite velocities may deviate frm that assumed in 
the data reduction, but efforts to quantify the effect experimentally have been successful 
(Hedin et al ., 1983). The model densities presented by Hedin et al. (1983) have, however, 
been normalized by a factor of 1.63 to bring them into agreement with the dag data. 
Since the nightside temperatures are much lower than the dayside, the midnight densities 
above 100 km are less than the noon densities, and above the homopause, diffusive 
separation is more effective, cuasing the nighttime mixing ratios of lighter constituents 
at a given altitude to be larger. For example, in the model for F l0 7 = 200 at the equator, 
the mixng ratio of O at 135 km is 7.3% at noon and 17% at midnight; atomic oxygen 
becomes the dominant component above 155 km at noon and above 140 km at mid- 
night. The noon and midnight models are shown in Figure 4. The variability of the 
neutral densities is found to be larger on the nightside than on the dayside: the standard 
deviations are 1 5-23 % on the dayside and 50-60% on the nightside (Hedin et al. , 1983 ; 
Keating et ai , 1985). Although the VTS3 model is based on obser/ations near high solar 
activity, a full rnage of solar activity variations has been incorporated into the model 
using the 10.7 cm flux at 1 AU(Fi 0 7 ) as a parameter. Thus, the low solar activity models 
should be considered somewhat speculative. It is interesting to note, however, that Kim 
et ai (1989) found that the solar cycle minimum electron density profiles measured by 
the PVO radio occultation experiment (ORO) are well reproduced by photochemical 
equilibrium models that employ the VTS3 model, but efforts to reproduce the solar cycle 
maximum ORO profiels suggest that the VTS3 densiteis are too low below 135- 140 km. 

The BNMS model can be compared to the model for conditions appropriate to the 
PV Bus entry: 39°, F 10 7 = 189 and 08 : 30 local time. The C0 2 densities at 100 km in 
the BNMS model are smaller by about 15% compared to the VTS3 model, but the 
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ALTITUDE (Km) 


Fig. 3. Total mass densities from the sounder probes as a function of altitude. Corresponding density 
profiles from the VTS3 empirical model are shown as solid lines. Taken from Hedin ttai (1983). 

thermospheric temperature are higher below and lower above 145 km. The C0 2 den- 
sities at about 160 km are 15% larger in the BNMS model. The mixing ratio of O at 
135 km is 7.9% in the VTS3 model and 7.0% in the BNMS model, the near equality 
indicating only that the BNMS values are based on ONMS data. The mixing ratios of 
N 2 at 135 km are similar: 6.3 and 8.3% for the BNMS and VTS3 models, respectively, 
but the CO mixing ratios differ by nearly a factor of three: 14.7% (BNMS) and 5.7% 
(VTS3). The He densities in the two models are equal to within 10%, but the mixing 
ratio is smaller by about 30% in the BNMS model. 
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ALTITUDE (Km! 

Fig. 4. VTS3 model densities as a function of altitude at noon and midnight at 0 e latitude. Taken from 

Hedin et ai (1983). 

Massie et ai (1983) constructed daytime (08 : 00 and 16 : 00 hr local time) and mid- 
night models of the densities of C0 2 , O, CO, N 2 , He, and 0 2 in Venus thermosphere 
over the altitude range 100 to 180 km by solving the one-dimensional continuity equation 
constrained by the PV in in situ measured neutral densities. No morning-evening asym- 
metry was included, although a bulge is observed in the light constituents near the 
morning terminator (see Section 4). The downward fluxes of N and O near midnight 
were inferred from the observed NO and 0 2 nightglow intensities (see Section 3), and 
the nightside density profiles of N and O were constrained to fit the measurements. 
Massie et al. preferred the BNMS measurements of C0 2 and N 2 , which were about 1.8 
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and 2.4 times the non-normalized values ONMS values, and the ONMS measurements 
of O and CO. Alternative values of O and CO that were 1.6 times the ONMS values 
were also given. The densities of 0 2 in their davside model are computed from the 
continuity equation assuming that the major production reaction is three-body recom- 
bination (Equation (3)) and loss is by photodissociation and three-body recombination 
of O and 0 2 : 

O + 0 2 + M - 0 3 + M . (4) 

At the lower boundary, the mixing ratio of 0 2 was set equal to 1 x 10 ’ 3 , that computed 
by Yung and DeMore (1982) from mesospheric chemistry. The dayside and nightside 
models of Massie et al. with non-normalized values for O and CO densities are shown 
in Figures 5(a) and 5(b). 

A comprehensive model of thermospheric composition, part of the Venus Inter- 
national Reference Atmosphere (VIRA), has been constructed by Keating etal (1985) 
based on information available through 1985. The thermosphere is divided into two 
regions, the lower region includes altitudes from 100 to 150 km and the upper region 
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is that above 150 km. The models for noon and midnight are shown in Figures 6(a-d). 
Diffusive equilibrium is assumed for all species in the upper region; the composition is 
based on the measurements of the ONMS, normalized to the OAD data. The normali- 
zation factors used in the VIRA models, 1.83 for C0 2 and 1.58 for O, are slightly 
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Fig. 6. Neutral number densities as a function of altitude (r) from the VIRA model, (a) Dayside 100 to 
150 km. (b) Dayside 150 to 250 km. (c) Nightsidc 100 to 150 km. (d) Nights idc 150 to 250 km. Taken from 

Keating et al (1985). 


different from the factor 1.63 used for all species in the VTS3 model. Normalization 
factor for CO and N 2 could not be derived and were arbitrarily assumed to be the same 
as for O. 
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The VIRA models nominally apply to moderately high solar activity F i01 = 150, but 
Keating el al. (1985) give an empirical expression for the change in the temperature and 
in the log of the number densities of C0 2 , N 2 , CO, 0, He, and N as a function of F l01 
and the 81 day mean value of F l01 . The expressions are based on the VTS3 model and 
on the drag data, but they are simpler than those used in the VTS3 model. The use of 
the 10 j 7 cm solar flux as an index of ultraviolet fluxes and thermospheric densities and 
emissions is widespread, but not firmly grounded. For example, Tobiska et al. (1988) 
show that upper thermospheric heating in the terrestrial atmosphere is dominated by 
absorption of chromospheric emissions, whereas F )0 7 correlates better with transition 
region and cool coronal emissions. 

The VIRA model can be compared to the VTS3 model appropriate to 16° N latitude 
“Hid F |0 7 = 150. The total number densities at noon and 135 km are comparable, but 
the O and CO mixing ratios are about 5% in the VTS3 model and somewhat larger, 
about 6% in the VIRA model. The VTS3 model for high solar activity (F l0 7 = 200), 
however, has larger mixing ratios of about 7% for O and CO at 135 km. Possibly the 
solar activity variations in the O abundance assumed in the VTS3 model are larger than 
those implied by the VIRA model, although there have been no in in situ measurements 
at low solar activity. In the VIRA model, the O density profile exhibits a maximum near 
110 km, as the photochemical model of Massie et al. (1983) shows, whereas the O 
densities continue to increase with decreasing altitude in the VTS3 model. In situ 
measurements of O are not available below the altitude of minimum periapsis of the PV 
Orbiter, near 140 km, and the theoretical profiles of O are preferred below 140 km. At 
midnight, the mixing ratios of O and CO at 135 km are 26 and 8%, respectively, in the 
VIRA model compared to 17 and 5%, respectively, in the VTS3 model. Both models 
give N ; mixing ratios near 100 km of about 4%, but because of different assumptions 
about the temperatures and eddy diffusion coefficients, at 135 km the noon mixing ratio 
of N 2 is 7.3% in the VTS3 model and 5.3% in the VIRA model; at midnight, the N 2 
mixing ratio is 9.8% in the VTS3 model and 8.6% in the VIRA model. 


3. Emission Features 


3.1. Introduction 

Atmospheric emission features are usually categorized as dayglow, nightglow, or aurora. 
Dayglow emissions are produced by the direct interaction of solar radiation or photo- 
electrons with atmospheric gases, and by prompt chemiluminescent reactions. Night- 
glow is luminosity produced by chemiluminescent reactions of species that are produced 
during the day or transported from the dayside. Auroral emissions are usually con- 
sidered to be those that are produced by impact of particles other than photoelectrons. 
Although Venus does not have an intrinsic magnetic field, ultraviolet emission features 
attributed to particle impact have been detected on the nightside of the planet in images 
from the PV orbiter ultraviolet spectrometer (Phillips et al., 1986; see Section 3.2). 
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The sources of dayglow emission features include photoionization and excitation 
X + hv-* X* * + e , ( 5 ) 

electron impact excitation 

X + e -*X* + e, (6) 

and electron impact ionization and excitation 

X + e -*X + * + 2e , ( 7 ) 

In the foregoing expressions, X is either an atomic or molecular species, and the asterisk 
(*) represent an excited state of the species that can decay to a lower state with the 
emission of radiation. In addition, dayglow emission features attributable to fragments 
or fragment ions of a molecular species A B may be produced by photodissociative 
excitation 


AB + hv^A +fl*, 
electron impact dissociative excitation 
AB + e-*A + B * , 

photodissociative ionization and excitation 


( 8 ) 

( 9 ) 


AV + hv-* A + B* * + e 

and electron impact dissociative ionization and excitation 


( 10 ) 


AB + e-*A + B* * + le . 


( 11 ) 


Excited species that may emit to the ground state by a dipole-allowed transition can 
be produced by resonance scattering of solar photons by an atomic species X 

X + h v — * X* — ► X + h v (12) 

and fluorescent scattering by a molecular species AB 

AB(v) + hv-+AB*(v')-»AB(v") + hv' . ( 13 ) 

The probability of an atom resonantly scattering a solar photon per unit time is called 
the g factor, which is related to the oscillator strength / by the expression 


g = (nF) — X 2 /, 
me* 


(14) 


where nF is the photon flux in units of photons cm ' 2 s “ 1 A ~ 1 and A is the wavelength. 
In these units, the constant ne 2 /mc 2 has the value 8.829 x 10“ 13 cm. For fluorescent 
scattering in a molecular band system, it is convenient to define an excitation rate q v 
that is the probability per unit time of producing the excited state in the vibrational level 
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v' by absorption of a photon 

q v = (kF)X 2 — f 0 . 
me 


(15a) 


or 


co' k* 

q v . = (ref) — — A v 0 , 
c o 8ttc 


(15b) 


where a>' and w" are the statistical weights of the upper and lower states and A v 0 is 
the absolute transition probability of the transition from vibrational level v* of the upper 
state to v = 0 of the ground state. The g factor for production of emission in the (u\ v") 
band is then given by 


Sv’v" 



(16) 


A pedagogical discussion of atmospheric emissions and their measurement can be found 
in Barth (1969). 


3.2. Atomic oxygen emissions 

The atomic oxygen emissions at 1304 and 1356 A arise from the transitions 

0( 3 S°) - 0( 3 />) + Av(1302, 1304, 1306 A) (17) 

and 

O( 3 5 0 ) -► 0( 3 7 > ) + A v(1356, 1358 A), (18) 

respectively. The 0( 3 5°) state is connected to the ground state by a dipole-allowed 
transition with an oscillator strength of 0.048 (Doering etal., 1985) so resonance 
scattering of solar radiation is an important source of the emission in the dayglow. Other 
sources include electron impact excitation of O, and electron-impact and photo-disso- 
ciative excitation of CO and C0 2 . The sources of the 1356 A emission are the same, 
except that, since the transition from the ground state of O to the s S° state is spin 
forbidden, resonance scattering is not important. 

Mariner 5 and Venera 4 searched for the 1304 A resonance triplet of O, under the 
assumption that photodissociation of C0 2 and diffusive separation should cause O to 
become the major constituent of the atmosphere at high altitudes. The 1304 A intensity 
is a measure of the integrated O density above the altitude of C0 2 absorption, although 
because the emission is optically thick, a radiative transfer model is necessary to 
interpret the measured brightness. No emission was initially reported for the Mariner 5 
spectrophotometer (Barth et a!., 1967) or the Venera 4 instrument (Kurt et al., 1968). 
Venera 4 observed the atmosphere below 300 km in shadow. Strickland (1973) 
suggested that the thickness of the airglow layer was not within the resolution of the 
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Mariner 5 instrument. According to von Zahn et ai (1983), the three channels of the 
Mariner 5 spectrophotometer were allowed to go off scale as they swept across the 
bright disk of the planet, and the oxygen channels did not recover before the bright limb 
had been crossed. Several years later, a reanalysis of the Mariner 5 airglow data by 
Anderson (1975) revealed the presence of the 1304 A triplet in the region of the 
terminator with an intensity of 670 R. By scaling to the subsolar point by analogy with 
terrestrial data, a nadir intensity of 1 1 kR was inferred and an O mixing ratio of 10% 
O was derived. 

In 1967 Moos etal. (1969) detected emission from Venus near 1300 A using a 
low-resolution rocket-borne spectrometer. Moos and Rottman (1971) later measured 
the FUV spectrum (1200-1900 A) of Venus using a moderate resolution spectrometer 
and reported preliminary values for the intensities of the O emissions. The psectrum is 
shown here in Figure 7. Rottman and Moos (1973) refined the analysis and reported 



Fig. 7. Rocket spectrum of Venus (points) compared with a best-fit synthetic spectrum. The ‘mystery’ 
feature at 1400 A was identified by Durrance et ai (1980) as the (14, 5) CO fourth positive band From 

Rottman and Moos (1973). 

disk intensities of 5.5 + 0.5 and 2.7 + 0.5 kR for the 1304 and 1356 A emissions, 
respectively. They estimated a mixing ratio for O of 10% at a C0 2 column density of 
4 x 10 16 cm' 3 , the C0 2 column density above the altitude of absorption reported by 
Strickland et al. (1972) for Mars. Their production compares well to values found from 
Pioneer Venus in situ measurements. For example, in the VIRA dayside model, the 
0/C0 2 ratio is 13% at a column density of 4 x 10' 6 cm * 2 . The agreement is, however, 
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largely fortuitous. The analysis of Rottraan and Moos was based on the Venus model 
of Strickland (1973), except for the cross sections for electron impact excitation of O 
which were adopted from Zipf and Stone (1971). These cross sections have recently 
been shown to be too large, by a factor of about 7 for the 1304 A emission and a factor 
of2forthe 1356 A emission (Stone and Zipf, 1974; Zipf and Erdman, 1985; Zipf, 1986). 
Rottraan and Moos found electron impact excitation to be the major source of both of 
the emissions. Strickland (1973) computed intensities for the emissions using a model 
atmosphere of Me Elroy (1969) with an exospheric temperature of 700 K. He could not 
reconcile the 1304 A intensity for the ratio of 1304 to 1356 A intensity of about 2 with 
the model even with several tens of percent O, based on excitation parameters (g values 
for resonance scattering and ‘effective’ g values for electron impact excitation) derived 
from analyses of terrestrial airglow and Mariner 6, 7, and 9 airglow from Mars 
(Strickland eta!., 1972). Stricfend’s (1973) assumed impact excitation rate for the 
1304 A emission was a factor of 3 to 4 less than that derived from the cross section 
reported by Zipf and Stone (1971), but only about a factor of 2 larger than the currently 
accepted value. Strickland suggested, among other possibilities, that another source of 
emission might be present in the 1356 A signal. Durrance et al. (1981) later showed that 
the (14, 4) CO fourth positive band contaminates the 1356 A signal. 

The Mariner 10 spectrophotometer also detected strong emission at 1304 A 
(Broadfoot et al., 1974); the brightness was 17 kR, about 3 times the Moos andRottman 
(1971) rocket values and an order of magnitude larger than the values reported from 
Manners 6 and 7 for Mars (Strickland etal . , 1972). Although the O abundance was 
difficult to model accurately, the conclusion that the atmosphere of Venus contained 
more O than that of Mars was inescapable. 

In 1978, the Venera 11 and 12 spectrophotometers measured disk brightness at 
1304 A of 6.4 kR, from which they inferred an atomic oxygen column density of 
1.9x10 cm . They also concluded that electron impact on O must be an important 
source of O( 3 S 0 ), since the intensity was not flat across the disk. The importance of 
electron impact can be evaluated by reference to Table II, where we present computed 
excitation rates of various excited states of atmospheric species produced in electron- 
impact processes for a model appropriate to high and low solar activity based on the 
neutral VTS3 models computed by Fox (unpublished calculations; cf. Fox, 1982). 
Table III shows the production rates due to photoprocesses for the same models The 
integrated electron-impact excitation rate of O( 3 S 0 ) for the high solar activity model is 

only 3.4 x 10® cm 2 s ', but the apparent emission rate will be enhanced by radiative 
transfer. 

Stewart etal. (1979) reported the first measurements of the UV airglow from the 
Pioneer Venus OUVS, including the oxygen emission lines; the noted that the 1356 A 
ermssion was limb brightened and ascribed this effect to the production of energetic 
( ) atoms in dissociative excitation of CO. Cross sections are not available for 

production of 0( 5 ) in electron impact or photodissociative excitation of CO. 
Table III shows that electron impact and photodissociative excitation of CO, are 
negligible sources of the 1304 A emission. In their study of the CO fourth positive band 
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TABLE II 

Integrated column production rates for excited electronic species due to 
electron impact for model Venus thermospheres based on Pioneer Venus data. 
The high solar activity (High SA) model is that of Hedm era/. (1983) for 
F - 200 and the low solar activity (Low SA) model for F w 7 - 74. The solar 
fl^es are the SC No. 21REFW and F79050N spectra of Hinteregger (private 
communication; see also Torr er at., 1979) for the low and high solar acuvtty 
models, respectively. 


Source Column production rate 

(10‘cm' I s"') 


High SA Low SA 


2.1(4) 7.0(3) 

1.1(4) 2.1(3) 

5 . 5 ( 2 ) 1 - 8 ( 2 ) 

4.0(3) 7.5(2) 

3.7(3) 1.2(3) 

2.2(3) 7.2(2) 

1.1(3) 2.0(2) 

3.5(2) 6.5(1) 

1.4(4) 2.1(3) 

6.2(3) 2.1(3) 

1.1(3) 17(2) 

4.0(1) 1-3(1) 

3.0(2) 4-6(1) 

2.0(0) 4.0(- 1) 

4.0(1) 1-3(1) 

8.8(2) 14(2) 

2.6(2) 7.9(1) 

8.8(1) 2.7(1) 


CO(a 3 /T) 

co 2 

CO 

CO(-4'/I) 

co 2 

CO 

co ;(A 2 nj 

co 2 

co ;(B 2 £;) 

co 2 

CO* (A 2 17) 

co 

CO*(B 2 Z*) 

CO 

0(‘D) 

0 

0('S) 

co 2 

0 

0(’S°) 

co 2 

0 

CO 

O CS°) 

co 2 

o 

n i(A 2 n.) 

N 2 

n;<b 2 i;) 

n 2 


system Durrance et ai (1981) found that the limb brightening at 1356 A is the result 
of contamination of the 1356 A emission by the (14, 4) CO fourth positive band at 

1352 A. w , . 

Meier et at. (1983) modeled the atomic oxygen emissions using a Monte-Carlo tecn- 

nique for solution of the radiation transport equations. They found that the incorpo- 
ration of partial frequency redistribution also enhanced the predicted limb brightening. 
Their model calculations reproduced the measurements if the O densities from the 
BNMS model were reduced by about 40%, but their analysis was based on cross 
sections for electron impact production of emission at 1304 and 1356 A reported by 
Stone and Zipf ( 1974) that were later reduced by a factor of 2.8 (Zipf and Erdman, 1985, 
Zipf, 1986). A model of the terrestrial dayglow oxygen emissions (Meier etal . , 1985) 
suggested that the Stone and Zipf cross sections should be scaled down by 40%. A 
reanalysis of the Pioneer Venus data using the renormalized cross sections showed good 
agreement between the O densities necessary to explain the measured emission and the 

BNMS model densities (Paxton and Meier, 1986). 

Remote sensing can be a reliable method of determining thermospheric densities, as 


SPAC 191/19 



20 


J. L FOX AND S. W. BOUGHER 


TABLE III 


Column production rates of some excited electronic states of species in the 
Venus thermosphere due to photodissociative xcitation or photoionization and 
excitation. The models used are as described in Table II. 


Species 

Source 

Column production rate 



( 10 * cm ~ 2 s ~ 

') 



High SA 

Low SA 

CO (a J /7) 

co 2 

18(4) 

7.9(3) 

CO(A'/7) 

co 2 

7.8(2) 

3.1(2) 

co;(A 2 nj 

co 2 

9.0(3) 

3.6(3) 

CO,* (it 2 !.* ) 

co 2 

1.0(4) 

4.1(3) 

CO*(.4 2 i7) 

CO 

3.4(3) 

7.4(2) 

co*(a 2 i*) 

CO 

13(3) 

25(2) 

0('2>) 

co 2 

1.6(6) 

6.5(5) 

O('S) 

co 2 

2.1(5) 

7.6(4) 

0( 3 S«) 

co 2 

1-2(1) 

6.1(0) 

n 

N 2 

2.0(3) 

7.8(2) 

N,* (•»!?) 

N, 

4.1(2) 

15(2) 


Meier and Anderson (1983) have suggested, but in practice, the use of 1304 and 1356 A 
intensities to infer atomic oxygen densities both in the terrestrial atmosphere and in 
planetary atmospheres has been plagued by uncertainties. The major sources of error 
in models are in the electrons impact cross sections and in the solar fluxes, although, 
if limb scans of an optically thick emission are available, the altitude distribution of the 
scatterer can be obtained from the shape of the intensity profile independent of the solar 
flux or the calibration of the instrument (e.g., Anderson et al., 1987). In this case, the 
solar flux can be derived from the intensity of the emission. The accepted values for the 
electron impact cross section have changed by nearly a factor of 3 in recent years, 
although there is now substantial agreement between the excitation cross sections 
measured using electron energy loss methods (Vaughan and Doering, 1986; Doering 
and Vaughan, 1986; Gulcicek and Doering, 1988; Doering and Gulcicek, 1989) and the 
renormalized emission cross sections. There is also a wide range in the values of the 
solar photon fluxes used in or derived form resonance scattering calculations, even 
allowing for solar activity variations. Both the rocket flights of Moos and Rottman 
(1971) and the initial PV OUVS measurements were during period ofhigh solar activity. 
Strickland (1973) required a photon flux integrated over the multiplet of 
5 x 10* cm “ 2 s _ whereas the calculations of Meier et al. (1983) were consistent with 
a total flux of 1.4 x 10 ,o cm- 2 s-». Link et al. (1988) required a 1 304 A irradiance of 
7 x 10 9 cm 2 s 1 to explain the terrestial dayglow experiments at high solar activity. 
Experimental measurements of the total photon flux in the 1300 to 1310 A range, of 
which the 1304 A triplet accounts for about 85-90% (Mount and Rottman, 1981), also 
vary widely. At high solar activity, reported values range from 9.5 x 10 9 cm ‘ 2 s ~ 1 from 
the F79050N spectrum of Hinteregger (private communication; see also Torr et al.. 
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1979) to 1.54 x lO 10 from rocket measurements (Mount and Rottman, 1981); at low 
solar activity from 4.8 x 10 9 from the SC No. 21REFW spectrum of Hinteregger to 
7.1-7. 3 x 10 9 (Mount and Rottman, 1983, 1985) to 1.21 x 10 lo cm’ 2 s" 1 from the 
Solar UV Spectral Irradiance Monitor (SUSIM) data (Van Hoosieref al , 1987). There 
are, however, independent reasons to believe that the Hinteregger fluxes are too low. 
Ogawa and Judge (1986) used a rocket-borne Ne ionization chamber to measure the 
integrated solar flux between 50 and 575 A, with an estimated error of 7.3%. Their 
results indicate that the Hinteregger fluxes in this wavelength region are low by a factor 
of 2 at low solar activity and by 30% at high solar activity. Given the uncertainties in 
the parameters used in modeling the oxygen emissions, densities derived from these 
measurements should probably not be considered reliable to better than a factor of 2. 

Emissions at 1304 and 1356 A have been detected by the Pioneer Venus OUVS in 
images of the nightside of Venus. The emissions appear in bright patches that vary in 
size and intensity; the sheer spatial and temporal variability of the emissions suggests 
that they are produced by particle impact. The morphology of the aurora has been 
described by Phillips et al (1986). Figure 8 hows a series of brightness images recorded 



Fig. 8. Series of brightness images from the PV OUVS for the orbital sequence of June 25-July 1, 1982 
with apo apsis near midnight. Note expansion and intensification of emissions to fill the nightside in orbit 
1301. Missing orbits were not imaged at 1304 A. From Phillips et al. (1986). 


by the OUVS for our nearly contiguous orbits. The intensity of the 1304 A emission is 
usually about 10 R, but values near 100 R have been recorded; the 1304/1356 ratio is 
about 6. Fox and Stewart (1990) have presented evidence that the emissions are 
produced by impact of very soft (~* a few eV) electrons. More energetic electrons would 
penetrate too deeply and produce high intensities of the Cameron bands of CO, which 
are not observed, and ion densities that exceed observed values. Both radiative recombi- 
nation of O * and precipitation of O * , the sources in the terrestrial equatorial nightglow, 
are ruled out for Venus because the 1304 and 1356 A emissions are produced with 
comparable intensities by these mechanisms and because the predicted intensities due 
to radiative recombination are less than 10" 3 R (Julienne et al, 1974; Abreu et al, 
1986). In addition, soft electrons have been observed in the Venus wake by both the PV 
orbiter retarding potential analyzer (ORPA) (e.g., Knudsen and Miller, 1985) and the 
plasma analyzers on board Veneras 9 and 10 (e.g., Gringauz et al, 1979). Fox and 
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Stewart have modeled the intensities that would be produced by precipitation of elec- 
trons with the same energy spectrum and fluxes as the downward traveling portion of 
the ORPA spectrum. The total fluxes required were determined to be between 8 and 
28 % of the ORPA values. 

The energies of the electrons producing the aurora are constrained to be low so that 
they lose their energy in the upper thermosphere where O is the dominant constituent 
of the atmosphere. If the electrons penetrate down to the part of the atmosphere where 
C0 2 is the most important constituent (below about 140 km) unobserved emissions, 
such as the CO Cameron bands, would be observed with large intensities and the O 
emissions would be excited less efficiently, producing a larger ratio of ionization to 
excitation for a given 1304 A intensity. An alternative explanation has been proposed 
by J* Grebowsky (private communication). He has suggested that the same effect would 
be realized if the depth of penetration of the electrons were limited by the presence of 
a horizontal magnetic field in the nightside ionosphere, rather than by their energy. 

The excitation of ground state oxygen atoms by electron impact, photodissociation 
and electron-impact dissociation of C0 2 and CO, and dissociative recombination of 0 2 + 
can produce O atoms in the l S and l D states, which can either be quenched or decay 
by the emission of photons. 0( l D) produces the atomic oxygen Ted line* in the process 

o CD) -> OCP) + Av(6300 A) . ( 19 ) 

Ninety-five percent of 0( 'S) decays on O CD) producing the atomic oxygen ‘green line’ : 

0('5)-0( , Z)) + Av^5577/l). (20) 

Only 5 0 /„ decays to the ground state producing an ultraviolet photon of wavelength 
2972 A 

0(‘S)-* 0( 3 /*) + Av(2972 A) . (21) 

Dayglow intensities of 730 kR, 48 kR, and 2.4 kR were predicted for the 6300, 5577, 
and 2972 A emissions, respectively by Fox (1978) and Fox and Dalgamo (1981) in a 
low solar activity pre-Pioneer Venus model. The major source of both 0( '£>) and 0( 'S ) 
is photodissociative excitation of C0 2 . Since pre-Pioneer Venus models differed from 
current models mostly in the mixing ratios of O and CO, the intensities should not be 
highly model-dependent. 0('£») is strongly quenched by C0 2 , so most of the O('D) 
produced does not radiate. The rate coefficient for quenching of O('S) by C0 2 

0( l S) + C0 2 + 0('Z>) + C0 2 (22) 

is 3.3 x 10 exp(- 1320/T) cm 3 s~ l (Atkinson and Wedge, 1972) and the quenching 
altitude, where the lifetime against collisional deactivation is approximately equal to the 
radiative lifetime, is about 100 km. The altitude profile for production of O('S), shown 
in Figure 9 is characterized by two peaks, the upper one due mainly to absorption of 
solar photons in the 1 100-1140 A range and the lower peak, near 1 15 km is due to 
photodissociation of C0 2 by Lx 

LeCompte et al. (1989) have reported intensities of the 2972 A emission measured by 
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VOLUME EMISSION RATE (ph cm 3 s 1 ) 

Fig. 9. Volume production rate of 2972 A emission from photodissociation of C0 2 by La (dashed line), 
from the rest of the solar spectrum (dotted line) and total (solid line). From LeCocnpte ft erf, (1989). 

the Pioneer Venus OUVS. The zenith intensity reported for orbit 187, for which high 
solar activity was appropriate, is 7.0 kR, compared to 2.4 kR predicted by the low solar 
activity model of Fox and Dalgarno (1981). LeCompte et al analyzed limb profiles from 
two orbits. They showed that photodissociation is the major source of 0( 2 S) below 
175 km. Above that altitude, dissociative recombination of 0 2 + and a contribution from 
the Cl 2967 A ( 5 S -► 2 P) emission dominate the production profile. Fairly good agree- 
ment was obtained between the measurements and the model at high altitudes, although 
the model intensity deviates from the measurements at altitudes below 130 km. They 
suggested that the difference arises from oscillations in the temperature profile that result 
from the passage of gravity waves. Gravity waves have been shown that affect terrestrial 
airglow emissions (e.g., Porter etai , 1974; Taylor etal y 1987). Emission from 0( l S) 
is sensitive to temperature because the rate coefficient for quenching of O('S) by C0 2 
(reaction 22) depends exponentially on temperature. Oscillatory structures appear in the 
temperature profiles inferred from neutral densities measured by several PV instru- 
ments, including the BNMS (von Zahn etal y 1980), the ONMS (Kasprzak etai, 1988), 
and by the accelerometers on the probes (Seiff et a/., 1980). The temperature profile that 
LeCompte et al use for their standard model is the VTS3 model of Hedin et at. (1983), 
which is not based on (ONMS) measurements below 140 km. It has been suggested that 
the temperatures below 140 km in the VTS3 model are too low (Keating etai , 1985; 
see also Section 2). LeCompte et al. tested a constant (positive) offset in the temperature 
profile from the VTS3 model, but did not consider a temperature profile in which the 
deviation increases as altitude decreases. Their assertion that the temperature profile 
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must return to the VTS3 profile at the lower boundary is not entirely convincing. As they 
admit, LeCompte et al. test deviations of the VTS3 temperature profile without changing 
the neutral densities in a consistent way. More important, the model intensity profile 
is normalized to the measurements at low altitudes, where the source of the signal is 
mostly Rayleight scattering. At 05 km the 2972 A volume production rate vanishes, so 
the sensitivity of the model to the quenching coefficient and, hence, to the temperature 
near that altitude is small. 

Visible spectrometers on Venera 9 and 10 searched for emissions on the nightside at 
6300 and 5577 A, but found none. Upper limits have been placed on the intensities of 
the green and red lines of 10 and 20-25 R, respectively (Krasnopol’sky, 1981, 1986, and 
private communication, 1988). In the nightglow, the emissions are produced mostly m 
dissociative recombination of 0 2 . Fox ( 1 990a) has combined a model for the vibrational 
distribution of 0 2 * on the nightside with rate coefficients computed by Guberman ( 1987, 
1988) for production of 0( '£>) and 0( 1 S ) in dissociative recombination of 0 2 (t>) from 
different vibrational levels, v. Altitude profiles of the volume emission rates are shown 
in Figure 10. The integrated overhead intensities are 1-2 R for the green line and about 
46 R for the red line. The intensity of the green line is sensitive to the vibrational 
distribution of 0 2 , but the yield of O('Z)) and, therefore, the emission rate of the red 
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line is less sensitive to the vibrational distribution. The predicted intensity of the red line 
exceeds the Venera upper limit for the particular model ionosphere chosen. The night- 
side ionosphere, at least at times of high solar activity, is produced mainly by transport 
of O* from the dayside (Knudsen etai, 1980; Spenner etal . , 1981; Cravens etal . , 
1983), although it has been suggested that at times of low solar activity, electron impact 
is the most important ion source (Gringauz et ai, 1979; Breus etai , 1985; Knudsen, 
1988). The nightside ionosphere is highly variable, sometiems disappearing nearly 
completely (Cravens etal . , 1982). The model ionosphere of Fox (1990a) corresponds 
to undisturbed, relatively ‘gull-up’ conditions for a moderate downward flux of O + of 
1 x 10 s cm ‘ 2 s~ 1 and may not be representative of any single measurement. 

3.3. O* EMISSIONS 

0*( 2 Z>) and O + ( 2 P) are metastable states of 0 + that lie 3.32 and 5.02, respectively, 
above the ground state of O * . Radiation of O + { 2 D) produces a doublet at 3726 and 
3728 A. 0 + ( 2 P) may radiate to the 0*( 2 Z)) state producing a doublet at 6319 and 
7329 A, or to the ground O + ( 2 S) state producing an emission feature at 2470 A. These 
states are produced in the atmosphere of Venus mainly in photoionization and electron 
impact ionization of O. The yields of 0*( 2 D) and 0*( 2 P) are about 37 and 31%, 
respectively, in photoionization and 42 and 22%, respectively, n electron impact 
ionization of O (see Kirby etal., 1979; Burnett and Rountree, 1979). Emission from 
0 + ( 2 Z>) in dipole forbidden and consequently its radiative lifetime is long, about 
2.4 x 10 4 s (Seaton and Ostcrbrock, 1957). Most of the loss of O *( 2 D) is via chemical 
reactions and the importance of this species is, therefore, in its effect on the ion 
chemistry. For example, charge transfer reactions of 0*( 2 D) have been shown to be 
important for the production of N 2 % CO% and N + (Fox, 1982). 

No measurements of intensities of the emission from O + ( 2 D ) and O * ( 2 P) have been 
reported. Altitude profiles of the volume emission rates computed for a Pioneer Venus 
high solar activity model by Fox (1982) are shown Figure 11. The integrated overhead 
intensities are given here in Table IV. 

The O* emission triplet at 834 A arises from the dipole allowed transition 
2s2p 4 4 P-> ls 2 p y 4 5°, so resonance scattering of solar radiation by O* is a potential 


TABLE IV 


Integrated overhead intensities (volume production rates integrated over the layer) of 
emissions arising from metastable states of O* and N (adapted from Fox, 1982) 


Transition 

Wavelength (A) 

Intensity (R) 

O*(V»)-.O*( J 0) + hv 

7319, 7329 

720 

o*(V»)-o*( 4 S) + hv 

2470 

200 

0’( 2 D)~. 0'( 4 S) + hv 

3726, 3728 

4.3 

N( : D)-» N( 4 S) + hv 

5200 

16 

N( 2 /»)^N( 2 D) + hv 

10,400 

2.5 

N( 2 P)~. N^S) + hv 

3466 

160 
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Fig. H. Volume emissions rates of the transitions 0*( 2 D)-0 *(*S\ 0*( 2 P)^0*( 2 D) and 
CT( 2 />)-Cr( 4 S) at 3726, 7319-7329, and 2470 A. From Fox (1982a). 


source, as are electron impact and photoionizaiion of 0( 3 P). The first measurements 
of the intensity of the O * emission line at 834 A were made by Venera 1 1 in December 
1978. Bertaux etal. (1981) reported a maximum brightness on the disk of 156 R, with 
fluctuations of less than 20% across the disk. They interpreted their measurements as 
showing that resonance scattering of solar radiation is not important, as has also been 
found for the Earth. In the terrestrial ionosphere photionization of O produces a low 
altitude source that is reonandy scattered by the F2 region 0 + ions at higher altitudes 
(e.g., Carlson and Judge, 1973 ; Feldman et al., 1981; Klumar et ai, 1983b). Presumably 
the same is true for Venus also. 

3.4. Atomic C and C + emissions 

From their low-resolution rocket spectrum of Venus, Moos etal (1969) reported the 
existence of a signal longward of 1500 A that they attributed tentatively to braodband 
fluorescence with an intensity of 10-100 kR. Marmo and Engleman (1970) proposed 
that the signal could be explained by the presence of 14 kR of emission in the 1657 A 
resonance line of atomic carbon. Although the suggestion was influenced by the 
erroneous belief that the thermosphere was largely decomposed under the influence of 
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threshold of the mass spectrometers on Pioneer Venus. McElroy and McConnell (1971) 
modeled the atomic carbon in the Venus atmosphere, but were unable to explain the 
large intensities using a model atmosphere with a thermospheric CO mixing ratio equal 
to the measured abundance above the clouds, 4.5 x 10 “ 5 (Connes et al 1968). Pioneer 
Venus in situ measurement showed that thermospheric mixing ratios are much larger 
than those in the mesosphere (see Section 2). The first successful models of the atomic 
carbon densities based on Pioneer Venus neutral models were constructed by 
KrasnopoPsky ( 1 982b) and by Fox ( 1 982b). These models showed that the major source 
of atomic carbon in the Venus thermosphere is photodissociation of CO. The most 
important sink is reaction with 0 2 (McElroy and McConnell, 1971): 

C + 0 2 -C0 + O, (23) 

which proceeds with a rate coefficient of 3.3 x 10" 11 cm 3 s" 1 (Braun et al. y 1969). 
There are, however, no measurements of the densities of 0 2 in the Venus thermosphere. 
For treated the 0 2 mixing ratio, / 0j as a free parameter and computed the atomic carbon 
density profiles shown in Figure 13. The profiles exhibit maximum C densities at 
145-155 km of 5 x 10 6 and 3 x 10 7 for 0 2 mixing ratios of 3 x 10" 3 and l x 10 " 4 , 
respectively. Fox also showed that photoionization of C is the major source of C* and 
used the measured C + densities (Taylor et ai % 1980) to infer that the mixing ratios of 



DENSITY OF C (cm“ 3 ) 


Fig. 13. Computed number densities of neutral atomic carbon as a function of altitude for 0 2 /C0 2 mixing 
ratio was adopted by the Venus International Reference Atmosphere (Keating et a/., 1985). Taken from Fox 

(1982b). 
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0, must be closer to the lower value. Krasnopol'sky ( 1 982b) computed C densities using 
an 0 2 mixing ratio of 2 x 10' \ based on the mesospheric photochemical model of 
Krasnopol’sky and Parshev (1981) and obtained similar C densities to those in the larger 
f Q model of Fox. 

Paxton (1985) modeled Cl limb intensities at 1657 and 1561 A made by tne 
PV OUVS, including the contamination of the signal by several CO fourth positive 
bands. Figure 12 shows a comparison of the measured and computed limb intensities 
of the 1657 A emission due to various sources. He reported C0 2 unattenuated photodis- 
sociative excitation frequencies of 2 x 10" 10 and 4 x 10" 10 s' 1 for the 1561 and 
1657 A lines, respectively; for photodissociation excitation of CO the correspndmg 
unattenuated frequencies were 1.8 x 1(T 9 and 2.7 x 10"* s'*. Resonance scattering 
of solar radiation was found to be the most important source of the emissions. Paxton 
used a radiation transport model to compute the resonance scatering intensities and to 
derive the atomic carbon densities. The best fit to the emission rates was obtained for 
a C density profile similar to the lower density profile of Fox (1982b), for which the 
0 2 /C 2 mixing ratio was 3 x 10” 3 . Thus a discrepancy exists between the C profiles 
derived from the C + densities and those derived from the Cl airglow. Paxton has 
suggested that the problem is with the ion chemistry. It is also possible that the rate 
coefficients on which the atomic carbon models are based are inaccurate. In any case, 
the CO photodissociation rates computed in the models above were based, in part, on 
low-resolution photoabsorption cross sections in the wavelength region longward of the 
ionization threshold at 885 A measured by Cook etai (1965). Recent evidence has 
shown that the absorption in this wavelength regin is due to discrete absorption into 
predissociating states rather than continuum absorption (Yoshino et ai, 1988; Letzelter 
et ai , 1987). Fox and Black (1989) constructed high-resolution photodissociation cross 
sections for rotational temperatures appropriate to the Venus thermosphere from a line 
list compiled by van Dishoeck and Black (1988) from information communicated 
privately to them by Stark, Smith, Yoshino, and Parkinson. They showed that the CO 
photodissociation rate, and the corresponding C production rate, are reduced by a factor 
of two over the values obtained using the low-resolution cross sections. This may have 
important consequences for the chemistry of atomic carbon in the Venus thermosphere 
that remain for future work to determine. 

3.5. CO AND CO - * EMISSIONS 

CO was first identified in the atmosphere of Venus from gorund-based infrared observa- 
tions of the first overtone band at 2.35 pm by Sinton (1963) and by Moroz (1964). 
Connes et ai (1968) observed several rotational lines in the CO fundamental (e = 1 -» 0) 
vibration-rotation band and derived a mixing ratio above the cloud tops of 4.5 x 10“ 5 . 
Young (1972) applied new information about the line widths to revise that value to 
5 5 x 10" 5 In situ measurements were carried out only below the clouds by the PV 
sounder probe gas chromatograph (LGC), which measured CO mixing ratios of 
20-32 ppm from 52 to 22 km (Oayama et ai, 1980) and by the gas chromatograph on 
Venera 12, which measured an abundance of 28 ppm (Gel’man et ai, 1979). 
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The first microwave detection of CO in the Venus atmosphere was that of Kakar et al. 
(1976), who observed the 0.26 cm (J = 0 - 1) pure rotational transition of CO. Micro- 
wave measurements are potentially more sensitive to lower pressures and therefore to 
higher altitudes than infrared measurements. Kakar eta/., as well as Schloerb et al. 
(1980) and Wilson et al. (198 1 ) observed a strong diurnal variation in the intensities, and 
concluded that the CO mixing ratio decreases by a factor of 10 from day to night in the 
range of total pressure 2-100 mb (about 65-85 km). Clancy and Muhleman (1985a) 
found that the CO abundance between 80 and 90 km peaks near 08 : 30 local time, but 
above 95 km, the peak was found to occur near the anti-solar point. The magnitude of 
the diurnal variation is a factor of 2-4 for both altitude ranges. Between 100 and 105 km, 
the derived mixing ratios are in the range 4x 1 O '* to 1 x 10 ~ 3 at night and 2x10 
to 5 x 10" 4 during the daytime. 

In situ measurements of thermospheric CO by the PV ONMS and BNMS show that 
the mixing ratio at 140 km is much larger than the microwave measurements indcate 
at 90 km (Niemann et al., 1980b; von Zahn et al., 1980; see Section 2). The thermosphe- 
ric model of Massie et al. (1983) predicts larger CO densities at the lower boundary than 
the measurements of Wilson (1981) or Clancy and Muhleman (1985a) indicate, but 
given the uncertainties the agreement is adequate. That the model does not show the 
observed diurnal variation in the CO densities may be the results of its one-dimensional 
nature. Clancy and Muhleman (1985b) have suggested that the high altitude diurnal 
variation is a result of the sub-solar to anti-solar circulation, the return flow of which 
is in the upper mesosphere. 

The fourth positive band system of CO is the most intense band system in the 1200 
to 1800 A region of the dayglow spectrum of Venus. It arises from the dipole allowed 
transition 

CO(/t'/7; u')-+ CO(A r| X; t?") + hv . (24) 

Potential sources include photodissociative excitation or electron impact dissociation 
of C0 2 : 

C0 2 + (e,hv;E> 13.5 eV)- CO(/t‘/7; r') + O + (e, ), (25) 

electron impact excitation of CO : 

COCVT; v") + e(E> 8.02 eV)- COM 1 /7; ,v') + e, (26) 

dissociative recombination of C0 2 * 

CO; + e - CO(/( 1 77; v) + O , (27) 

and fluorescent scattering by ground state CO, 

COOVT; o) + hv -* CO(A l II; v')-*CO(X l Z, c") + hv' . (28) 

The branching ratio for production of COM 1 IT) in dissociative recombination of C0 2 + 
(reaction 27) is about 5% (Gutcheck and Zipf, 1973). 

The observed signal longward of 1400 A in the rocket spectrum reported by Rottman 
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and Moos (1973) was interpreted as being due to the atomic carbon features and to the 
CO fourth positive band system. The total intensity of the band system was determined 
by fitting a synthetic spectrum to the measurement (see Figure 7). Except for an uniden- 
tified feature at 1400 A (later identified by Durrance et al. (1980), see below) a good fit 
to the measured spectrum was obtained and a brightness of 25 ± 5 kR was reported. 
The process(es) responsible for the observed emission could not be unequivocably 
determined, but the source due to fluorescent scattering was judged to dominate the 
production, and a CO mixing ratio of 10% at a C0 2 column of 4 x 10 16 cm' 2 was 
inferred. The PV ONMS later confirmed the implication that the atmosphere of Venus 
contains larger amounts of CO than the atmosphere of Mars. The VTS3 model for 
F 10 7 = 200 and 12 : 00 hr local time exhibits a CO mixing ratio of 12% at 137.5 km, the 
altitude at which the C0 2 column density is about 4 x 10 16 cm“ 2 . 

The Mariner 10 EUV spectrophotometer recorded a brightness at 1480 A of 55 kR 
that was attributed to emission in the CO fourth positive bands (Broadfoot et al . , 1974). 
The spectrophotometers on Venera 1 1 and 12 also recorded the intensity at 1500 A (see 
Table I). Significantly smaller disk brightnesses of 2. 1 and 2.7 kR for Venera 1 1 and 12, 
respectively, were reported (Bertaux et al . , 1981). The disagreement with the Mariner 10 
data was attributed by Bertaux etai (1981) to contamination of the airglow signal by 
Rayleigh scattering. 

Durrance et al. (1980) presented observations of the Venus dayglow from 1250 to 
1430 A made by the PV OUVS and Durrance et al. (1981) measured the Venus dayglow 
spectrum from 1280 to 1380 A at high (0.4 A) resolution with the International Ultra- 
violet Explorer (IUE) satellite. The strongest features in this region of the spectrum, 
other than the atomic oxygen emission at 1304 and 1356 A, were shown to be due to 
the (14,i?") progression of the fourth positive band system excited by fluorescent 
scattering of solar La. As ponted out by Kassal (1976), this source is larger than 
fluorescent scattering in all the rest of the fourth positive bands for CO column densities 
larger than 1 x 10 17 cm“ 2 . The overlap of the rotational lines of the (14,0) fourth 
positive band with the solar La line are shown in Figure 14. The unidentified feature at 
1400 A in the rocket spectrum of Rottman and Moos (1973) (see Figure 7) was 
identified by Durrance etai (1980) as the (14, 5) fourth positive band. 

The C0 2 absorption cross section is about 8 x 10“ 20 cm 2 at 1216 A (Nakata et al ., 
1965); consequently solar La penetrates to a column density of about 1.25 x 10 19 , or 
to about 1 10 km. The C0 2 cross sections in the region of the CO emissions, 1250 to 
1400 A, are in the range 2-9 x 10“ 19 cm 2 , so the emission bands will be absorbed if 
they are produced below about 120 km. Because the (14,5) bands at 1400 A is 
unblended and produced almost entirely by fluorescent scattering, measurement of the 
intensity of this feature could provide a valuable remote sensing technique for deter- 
mining CO abundances in the Venus thermosphere (Durrance etai . , 1980). 

Durrance (1981) presented spectra with a resolution of 4 A recorded by a sounding 
rocket telescope and spectrometer and constructed synthetic spectra using high-resolu- 
tion solar fluxes. He showed that the (14, 3) and (14,4) bands contaminate the O 
emission features at 1304 and 1356 A; potentially 67% of the 1356 A feature was 
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WAVELENGTH (A) 

Fig. 14. The rotational lines of the (14, 0) CO* fourth positive band showing the overlap with the solar 
La line. Taken from Durrance et al. (1980), 

estimated as arising from CO, thus bringing the ratio of 1304/1356 A intensities to about 
8, in agreement with the terrestrial value. Fox and Dalgamo (1981) presented altitude 
profiles of the sources of the CO fourth positive band system from a pre-Pioneer Venus 
model shown here as Figure 15. Ninety percent of the observed intensity of the band 
system is due to fluorescent scattering, with the other four sources about equally to the 
observed intensity. The computed integrated overhead volume emission rate is about 
2° kR. 

The Cameron band system of CO appears in the 1800 to 2600 A region of the 
spectrum. It arises from the dipole forbidden transition 

CO(o 3 /7; v')->CO(X'Z; v ") + Av. (29) 

The CO(a 3 /7) state is metastable with a relatively short lifetime of about 8 ms (Johnson, 
1972; Lawrence, 1971). The potential sources are the same as those for the fourth 
positive band system, except for fluorescent scattering. The Cameron band system is 
the brightest feature in the Martian ultraviolet airglow measured by the UV spectrome- 
ters on Mariners 6, 7, and 9(e.g., Stewart et al. , 1972; Barth et ai , 1971). The PV OUVS 
longer wavelength channel measured spectra in the region of the Cameron bands, but 
no observations or analyses have been reported. Table V gives intensities predicted by 
a pre-Pioneer Venus model of Fox and Dalgarno (1981) and unpublished calculations 
for a model based on Pioneer Venus data. The source due to dissociative recombination 
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Fig. 15. Altitude profiles of the sources of the CO fourth positive band system. Fluorescent scattering of 
La in the v' » 14 progression is responsible for 90 % of the intensity. The shoulder above the peak in the 
fluorescent scattering source is due to absorption of shorter wavelength photons. Taken from Fox and 

Dalgarno (1981). 


TABLE V 

Computed Cameron band overhead integrated intensities due to 
different sources for high solar activity (high SA) and low solar 
activity (low SA) models 


Source 

Intensities (kR) 

Low SA 1 

High SA b 

C0 2 + hv 

6.7 

20.8 

e + C0 2 

6.7 

20.9 

e + CO 

1.3 

11.1 

CO + e 

5.1 

4.5 

Total 

20 

57 

• From Fox and Dalgarno (1981) (pre-Pioneer Venus model). 

Pioneer Venus model, based on neutral densities from Hedin 

eiat. (1983) (cf. Fox. 1982, 1985). 



SPAC 191 33 


ORIGINAL PAGE IS 

OF POO* quality 











34 


). L. FOX AND S. W. BOLGHER 


of C0 2 " is overestimated in the pre-Pioneer Venus model because the atomic oxygen 
densities were too low and consequently, the densities of C0 2 + are too large. 

The source due to electron impact excitation of C0 2 is uncertain, because the 
normalization of the cross section is a matter of some controversy. Excitation functions 
have been measured by several workers but reported values for the maximum value of 
the cross section range from 3 x 10” 17 to 3 x 10” 16 cm 2 (Freund, 1971; Ajello, 1971; 
Wells etal . , 1972; Erdman and Zipf, 1983). Conway (1981) constructed a synthetic 
spectrum of the Martian dayglow from 1800 to 2800 A, and by comparing the model 
and measured intensities, found that a cross section with a maximum value of 
7 x 10“ ,7 cm 2 was consistent with the data. Taking into account new information 
about the lifetime of the CO(a 3 J7) state, Erdman and Zipf (1983), however, reported 
that laboratory data indicated a maximum value of the cross section closer to 
2.4 x 10“ 16 cm 2 , more than three times the Conway value. The electron-impact exci- 
tation rates in Tables II aand V are based on the cross sections of Conway. Analysis 
of the Pioneer Venus Cameron band data, if they are obtainable, might contribute to 
a resolution of this issue. 

The CO * comet tail and first negative bands arise from the transitions :A 2 I1-*X 2 Z + 
and B 2 l-> X 2 kl&* > respectively. The thresholds for production of the A and B states 
are 2.53 and 5.65 eV, respectively (Herzberg, 1950). Therefore, the comet tail bands 
appear in the visible spectrum and the first negative bands in the ultraviolet. Potential 
production mechanisms include photoionization and electron impact ionization of CO, 
photodissociative and electron impact dissociative ionization of C0 2 , and fluorescent 
scattering of solar radiation by ground state CO * . No emissions have been observed 
that are attributable to either band system. Paxton (1988) has proposed that observa- 
tions of these band systems by instruments on future spacecraft could help to determine 
the fraction of the mass-28 ion densities (which were measured by the PV OIMS (e.g., 
Taylor et al y 1980)) that is due to CO * , rather than to N 2 * . He has predicted a nadir 
intensity of the (0, 0) first negative band at 2191 A of 615 R; for the (3, 0) comet tail 
band at 401 1 A, the predicted nadir intensity is 45 kR. Tables II and III show that for 
high solar activity the predicted integrated overhead intensities due to photoionization 
and electron impact ionization of CO alone are 4.5 kR and 1.6 kR for the comet tail 
and first negative band systems, respectively. 

3.6. C0 2 * EMISSIONS 

The Fox-Duffendack-Barker {A 2 F7 U -+X 2 n g ) band system consists of several narrow 
bands between 2800 and 5000 A. The ultraviolet doublet at 2883 and 2896 A is produced 
by the (0, 0, 0)-(0, 0, 0) transition between the B 2 !* state and the X 2 IJ g ground state. 
Both band systems are produced by photoionization and electron impact ionization of 
C0 2 and by fluorescent scattering of sunlight. They appear as prominent emissions in 
the Mariner 6, 7, and 9 UV spectra of Mars (Barth et ai y 1971, 1972) and should be 
present in Venus spectra as well. No intensities have, however, been reported from the 
Pioneer Venus OUVS. 

Intensities of both C0 2 * band systems were predicted by Dalgamo and Degges 
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(1971). They computed g values for fluorescent scattering at Venus of 4.9 x 10 " 2 and 
5.2 x 10" 3 for the A and B -*X band systems, respectively. The total sub-solar 
zenith intensities were 107 kR and 31 kR, with 87 and 9 kR from fluorescent scattering 
for the Fox-Duffendack-Barker band system and the ultraviolet doublet, respectively. 
The contribution from fluorescent scattering was overestimated because at that time 
C0 2 * was assumed to be the major ion. 

The low solar activity model of Fox and Dalgamo (1981) predicted a total intensity 
of 15 kR for the Fox-Duffendack-Barker band system and 9.2 kR for the ultraviolet 
doublet. The most important source of both emissions is photoionization of C0 2 . The 
branching ratios in photoionization are uncertain, since photoelectron spectroscopy and 
fluorescence measurements imply different values (Samson and Gardner, 1973; 
Wauchop and Broida, 1972; Lee and Judge, 1972; Leach etai , 1978a). It has been 
suggested that a mixing of the A and B states occurs before emission, so that about 50% 
of the ionization into the B state leads to emission in the A state (Samson and Gardner, 
1973 ; Leach et al . , 1978b). The predicted intensities based on the fluorescence branching 
ratios are 19 kR for the Fox-Duffendack-Barker bands and 5.4 kR for the ultraviolet 
doublet. Revised calculations based on a Pioneer Venus model for high solar activity 
are shown in Table VI. The integrated overhead intensity of the Fox-Duffendack- 

TABLE VI 


Computed integrated volume emission rates of C0 2 features in the Venus dayglow for a 
model (high solar activity) based on Pioneer Venus data 


Source 

Intensity (kR) 

A 2 n„->x 2 n t 

b 2 z; ->x 2 n t 

Phoioionization of C0 2 

9.0 (14.1*) 

10.2 (5.1*) 

Electron impact ionization of C0 2 

3.7 

2.2 

Fluorescent scattering 

4.7 

0.25 

Total 

17.4 (23*) 

12.6 (7.5*) 


* Computed assuming a 50% crossover of B to A before radiating. 


Barker band system is 17.4 kR, with 12.7 kR from electron impact ionization and 
photoionizeation and 4.7 kR from fluorescent scattering. For the ultraviolet doublet the 
predicted intensity from photoionization and electron impact ionization is 12.4 kR and 
that from fluorescent scattering is negligible, only 0.25 kR, for a total integrated over- 
head intensity of 12.6 kR. If a 50% crossover from B to A occurs before radiating, the 
total intensities of the B-+X and A->X band systems would be 7.5 and 23 kR, 
respectively. 

3.7. He and He* emissions 

The He abundance in the lower atmosphere has not been measured directly. The large 
probe mass spectrometer (LNMS) had the capability of measuring He densities, but the 
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mixing ratio indicated 0 h«) by the LNMS data, about 450 ppm, was considered 
unreliable because of contamination by the He-N 2 mixture used to fill the probe to 
detect leaks (Hoffman et al ., 1980a; Donahue and Pollack, 1983). Thus bulk atmosphere 
mixing ratios have been obtained only by extrapolation from upper atmospheric values 
obtained from measurements by the PV BNMS and ONMS instruments and from 
airglow analyses. Because the mass of He is much less than the average molecular 
weight, this procedure is extremely sensitive to assumptions about the profiles of eddy 
diffusion coefficients and temperatures. The eddy diffusion coefficients presented by 
von Zahn et al . (1980) were based on a fit to the He profile measured by the BNMS 
above 130 km (see Section 2). The He homopause was lcoated at 130 km, the lower 
boundary of the measurements. The momingside model gives a value for / He at this 
altitude of about 50 ppm. An extrapolation to 100 km yields a mixing ratio of about 
12 ppm. The VTS3 model is based on largely ONMS measurements, but the He profile 
below 130 km is based on BNMS values. Nonetheless, Hedin et al. (1983) found the 
altitude of the He homopause to be lower, about 125 km, and the lower atmosphere 
mixing ratio to be correspondingly lower; values of 1-2 ppm are obtained at the lower 
boundaries of the noon and midnight models. Massie et al (1983) tested the sensitivity 
of their model to the value of the coefficient A in the numerator of the von Zahn et al 
(1980) expression for the eddy diffusion coefficient (Equation (1)); they found that 
changing the value from 0.8 to 2.0 changed / He for the lower atmosphere from 2 to 
9 ppm. For the valued ■ 1.4 x 10 l 3 preferred by von Zahn etal, they derived a helium 
mixing ratio of 5 ppm. 

Thermospheric He densities above 140 km given by the VTS3, BNMS, and VIRA 
models are, however, in fairly good agreement. He densities at 150 km for solar zenith 
angles less than about 60° are 4-5 x 10 6 cm -3 . ONMS data show a strong diurnal 
variation in the He densities that is out of phase with the heavier constituents, with larger 
densities at night than during the day and a distinct bulge in the He densities in the dawn 
sector (Hedin et al , 1983). Figure 16 shows the densities of species included in the VTS3 
model at 100 and 150 km as a function of hour angle. The VIRA model for 04 : 00 hr 
local time shows a density at 150 km of 1.2 x 10 8 cm' 3 ; closer to midnight the density 
is less, about 7 x 10 7 cm~ 3 . 

He resonance radiation at 584 A was detected by the both the Mariner 10 and 
Venera 11 and 12 spectrophotometers (Broadfoot et al , 1974; Bertaux etal, 1981). A 
maximum limb intensity of 600 R was measured by Mariner 10. The Mariner 10 bright 
limb profile from 200 to 800 km was analyzed by Kumar and Broadfoot (1975). They 
used a radiative transfer model that assumed spperical-symmetry and complete 
frequency redistribution to fit the shape of the limb profile. The exospheric temperature 
and the g factor were taken to be free parameters. The ‘best fit’ exospheric temperature 
was 375 ± 105 K. Takacz et al (1980) claim that a reanalysis of the data suggests that 
a temperature of 270 K is a better fit to the data. A look at Figure 17, taken from Kumar 
and Broadfoot (1975), does seem to suggest, with the full benefit of hindsight, that a 
lower exospheric temperature would be an equally good fit. The line center solar flux 
derived from the best fit g factor was 2 x 10 10 cm " 2 s ' 1 A " 1 at 1 AU. A He density 
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HOUR ANGLE (hrs) 

Fig. 16. Neutral number densities as a function of local time at 100 and 150 km, 0° latitude from the VTS3 
model. Note that the dayiime/nighttime ratio increases with increasing mass of the species. For helium (and 
hydrogen) the densities are larger on the nightside and peak in the predawn sector, forming the ‘helium (or 
hydrogen) bulge*. Taken from Hedin et at. (1983). 

of 2 x 10 6 cm “ 3 was derived at the assumed homopause (145 km), where the density 
of C0 2 was assumed to be 2 x 10 11 cm' 3 . Von Zahn et al. (1983) comment that if the 
lower (and more accurate, according to our present knowledge) temperature is assumed, 
the He abundance is almost exactly that determined by PV. We find that the column 
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Fig. 17. Mariner 10584 A brightness as a function of distance from the limb. The points were derived by 
subtracting the background noise level. The dashed lines represent altitude variations based on a single 
scattering model. The solid lines are the result of radiative transfer calculations including multiple scattering. 
The calculations are normalized to the peak observed brightness. From Kumar and Broadfoot (1975). 


density above 145 km, the presumed C0 2 absorption altitude, deduced from the 
Mariner 10 measurements is about 2 x 10 13 cm " 2 , whereas the coresponding value for 
the BNMS model is about 2.9 x 10 13 cm" 2 . Because the exospheric temperature is 
lower than Kumar and Broadfoot assumed, the C0 2 absorption altitude is, however, 
lower in Pioneer Venus models. The C0 2 absorption cross section at 584 A is about 
2.3 x 10" 17 cm" 2 , so the absorption altitude is about 140km in the V1RA dayside 
model and 136 km in the BNMS model. The He column density above 136 km in the 
BNMS model is about 3.7 x 10' 3 cm" 2 . 

The He 584 A disk brightnesses measured by Venera 1 1 and 12 were presented by 
Bertaux et ai. (1981). The maximum intensity on the disk was 270-280 R. Several 
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analyses of the Venera data have appeared. Kurt et al. (1983) modeled the intensity 
across the disk, using the He optical depth as a free parameter. The best fit He optical 
depth, in the range 2-4, is independent of the solar flux and the instrument calibration. 
The derived He column density above the C0 2 absorption altitude, taken to be 140 km, 
is about 1. 9-2.1 x 10' 3 cm" 3 , almost a factor of two less than the value from the 
Pioneer Venus models. Bertaux et al. (1985) presented a summary of similar results; a 
helium optical depth of 2.5, corresponding to a density of 1.6 x 10 6 cm" 3 at 150 km, 
was required to fit the intensity profile across the disk. The implied column density is 
smaller than the PV value. Chassefiere et al. (1986) revised the required optical depth 
to 3.5 ± 1.5, for a number density at 150 km of 2.2 x 10 6 cm -2 . They suggest that the 
discrepancy with the Mariner 10 data may be due to poor absolute claibration of the 
Mariner instrument. Krasnopol’sky (1983a) has suggested that the Mariner 10 instru- 
ment had a more serious problem with stray light than the Venera instruments had. 
Although the hydrogen densities derived from Mariner 5 and 10 La data are a factor 
of 2 larger than the values derived from PV OUVS data (Paxton et al., 1985, 1988) and 
from the in situ PVOIMS H + densities (Brinton etai, 1980), Paxton etal. (1985) 
showed that the fault lay with the use of an isothermal spherically-symmetric radiative 
transport scheme. They also suggested that fly-by observations do not provide enough 
constraints to determine the distribution of scatterers. This may be the case for the 
He 584 A data also. Indeed, given the large solar zenith angle variation of He densities, 
especially the existence of the pre-dawn bulge, one would expect a spherically-symmetric 
radiative transfer model, as all of the above models are including that of Kumar and 
Broadfoot (1975), to be a poor approximation. Furthermore, the analysis of Kurt et al. 
(1983) is isothermal and assumes a temperature in the scattering region of 300-400 K, 
whereas the dayside temperature varies from about 200 K near 140 km to about 275 
above 1 80 km. 

The solar fluxes at the 584 A line center that are assumed or derived are different even 
for similar solar activities, illustrating the inadequacy of the analyses and/or a calibration 
problem with the instruments. Bertaux etal. (1981) determined a value from the 
Hinteregger 79050N solar fluxes, although there is some evidence that the Hinteregger 
fluxes are too low, by up to a factor of 30% at high solar activity (Ogawa and Judge, 
1986; see Section 3.3). If so, the problem would be exacerbated, because a larger solar 
flux implies smaller He densities. Furthermore, Hinteregger presents the photon fluxes 
in the lines as integrated values and it is necessary to assume a lineshape to derive the 
line center fluxes; Bertaux et al. (1981) state merely that an appropriate line shape was 
assumed, but give no details. Doschek et al. (1974) found the shape of the 584 A line 
to be approximately Gaussian with a slightly flattened core. Ogawa et al. (1984) report 
a value for the full width at half maximum of about 0. 1 A from a rocket measurement, 
in good agreement with previously measured values. For a Gaussian shape with a width 
of 0.1 A, the line center solar flux implied by the Hinteregger integrated flux would be 
about 4.9 x 10 lo cm" 2 s“' A "'.Other analyses fit the shape of the intensity profile and 
derive a value for the solar flux at line center that reflects the uncertain calibration of 
the instrument. Kurt etal. (1983) derive a value of 3.3 x 10 ,o cm" 2 s" 1 A"', which 
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they stated— T . ** 

Ss^uo' 1 report that a column dens.ty abmn a .had oHhe 

Sltagiv^ 

fx”lO'C-H sotl activity vanatioa de.ennh.ed fra. .he (high solar 

acthity) F79050N and .he (low tola, acdvity) SC No. 21REFW tola, sp«ra o 
Hinleregger it a (actor of about 3. although d the Hinteregger tores are low by. actor 
”, a ,Tw tola, acdvity and b, 30% at h.gh tola, activity (Ogawa and Judge, 1986), 

«u« cnlar activity variation is reduced to a factor of 2. , , 

The He* line at 304 A was not detected by Manner 10, although one of the chann 
wa ?de"gned to detect emissions at that wavelength. It was, however, measured b 
Venera lfand 12 due to the greater sensitivity of the spectrophotometer. Bertaux el . 
0981) computed a column dens.ty of He * of2 x 10' 1 cm - 2 , which islarger b V 
of ^Le than the column density of about 2 x 10“ cm - implied by thePV OiMS 

measurements. 

3.8. N 2 and N 2 * emissions 

No emission arising from the band systems of N 2 has been detected in the Venus 
dayglow even thought N 2 comprises 3-4% of the atmosphere. The mtensmesofthe 

“ n *Hnnfield (a'IJ -tX'Z*) band systems were predicted by Fox and 
n7™!Sn These Ldl Ihe up'per slates of which are connected to the ground 

tadon of ground state N, in die dayglow, of the terrestrial planets. The V ^" d _, Kap “■ 
tauon u g. vman-Biree-Hopfield bands emit in the ultraviolet. The (1, 9) 

Vegard-Kaplan band at 3199 A, which falls between the (2, 0) and (3. 0) bands of the 
Fcni-Duffendack-Barker band system, and which may have been obs erred n the 
Ma«Undayglow speclr.m(Fox and Dalganto, 1979), was p r'jf?“^““^ 
ollir. and shouldbe detectable. The (3,0) Lyman-Buge-Hopfield band at 334 A 
, . th r\ 1355 a feature; the predicted intensity at 1354 A 

may contaminate the signal at the O 1335 A mature, ui P ra]culalions based 

rtheOUSbA en^sion^cros's Sections rf “one a^d Zipf (1974) indicated Uiat the 
1354 A intensity would contribute little to the observed 1356 A signal. Revised c cu- 

lations, using the renormalized 1356 A electron 

. c ,^on nQRM dve a lower integrated intensely for the 1355 A emission ui <iuu 

150 fallow solm^acdvhy (see Table H). Hence, the (3, 0) Lyraan-Birge-Hopfield band 
may contribute about 15% of the observed intensity yIrMbandsystem scan 

The N," Meinel (A 2 fl u -*X 2 L * ) and first negative (^ 2 ^» ^J'^XZL 
be produced bv photoionization and electron .mpact mmzat.on of N 2 , by fluorescent 
scattering of sol^ radiation, and by electron induced fluorescence (Degen, 1981). The 
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latter source is negligible in the dayglow. The Meinel band system occurs in the 
infra-red; the most intense band in the first negative band system is the (0, 0) band at 
3914 A. No emission at that wavelength was recorded by the visible spectrometers on 
Venera 9 and 10. An upper limit of 45 R has been placed on the overhead integrated 
intensity at 3914 A (Krasnopol’sky, 1978). 

Fox and Dalgamo (1981) predicted the intensities of the Meinel and first negative 
bands due only to photoionization and electron impact ionization. Paxton (1988) has 
suggested that intensities of the N 2 * first negative bands could be used to determine the 
fraction of the mass-28 ion densities due to N 2 * • He has predicted an overhead intensity 
for the 3914 A band of 17 kR, using the neutral Hedin model appropriate to the 
conditions of Pioneer Venus orbit 185, but this prediction must be in error. Fox and 
Dalgamo (1983, 1985) have computed the vibrational distribution of N/ in the 
terrestrial and Martian atmospheres, where vibrationally excited N 2 + is produced mainly 
in fluorescent scattering of solar radiation in the first negative and Meinel band systems. 
The same model has been adapted for the ionosphere of Venus (cf. Fox, 1982a) and 
the emission rates in the N 2 * band systems computed. The results are presented in 
Table VII. The integrated overhead intensity of the entire first negative band system is 


table vii 

Computed overhead integrated intensities of the Meinel (A — X) and first 
negative (B->X) band systems of N 2 * for a high solar activity (F l07 - 200) 

Hedin model 


Source 

Intensities (kR) 
Meinel 

First negative 

Photoionization 

2.0 

0.41 

Electron impact 

0.26 

0.00088 

Fluorescent scattering 

1.7 

1.4 

Total 

4.0 

1.8 


only 1.8 kR. The integrated overhead intensity of the 3914 A band is 790 R; for the (1, 0) 
Meinel band the predicted intensity is 910 R. Fluorescent scattering is the major source 
of both band systems; only about a quarter of the emission of the 3914 A is from 
photoionization and electron-impact ionization. 

3.9. N EMISSIONS 

The N( 2 D) and N( 2 P) states of atomic nitrogen are metastable, with radiative lifetimes 
of about 10 s and 12 s, respectively. They are produced in photodissociation and electron 
impact dissociation of N 2 and in dissociative recombination of N 2 * . The transition 
N( 2 T>) — ► N( 4 S) + hv produces a photon of wavelength 5200 A; the transition ofN( 2 P) 
to the N( 2 D) state leads to emission at 10400 A; the transition to the ground N( 4 S) 
state leads to emission at 3466 A. Like O * ( 2 D) and O *( 2 P), the importance of these 
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metastable states lies in their effect on thermospheric and ionospheric chemistry, ra er 
than in strong emission rates (Fox, 1982a). No measurements of the emissions have 
been reported, but the intensities have been modeled for a high solar activity Pioneer 
Venus model by Fox (1982a). The computed integrated overhead intensities are given 

in Table IV . 

3.10. NO NIGHTGLOW 

THe ultraviolet night airglow was first mapped by the PV OUVS during the I 
period of the mission (Stewart and Barth, 1979; Stewart etal., 1980; Gerard er a/., 
1981) The airglow was detected and identified as the gamma and delta bands of mine 
oxide (Feldman « a/., 1979; see also Stewart and Barth, 1979) excited by the radiative 
recombination of N and O atoms. Since no emissions suggestive of a nightside source 
of N and O atoms were present, it was proposed that the airglow was the result of these 
atoms being transported from their source on the dayside to the mghtside where 
recombination occurs. Thus the morphology and brightness of the airglow can provide 
a sensitive tracer of the thermospheric circulation. Other PV measurements P^ularly 
the hydrogen and helium nightside bulges (Taylor et al, 1984; Hedm et al., 1983) (s 
Figure 16), substantiate this concept of a strong day-to-night Venus thermospheric 
circulation (Mayr et al., 1980). Modelling efforts have been very successful in validating 
the proposed transport mechanism yielding the NO nightglow distribution and intensity. 
However, chemical questions remain regarding the dayside N-atom production 
available for transport to the nightside (Bougher et al., 1990a). .... 

Planet-wide observations by the OUVS at solar maximum initially showed that the 
NO(0 1 ) 6-band airglow typicaly exhibits a bright spot reaching 2.8 kR near 02 : 00 local 
time just south of the equator (Stewart er ai, 1980). This statistical averaged peak 
brightness is roughly four times the dark disk average vertical emission of 780 R as 
calculated for 25-orbits sampled early in the PV mission (Stewart etal 1980). The 
intensities reported by Stewart et al. (1980) have been revised as described by Bougher 
etal (1990a). Briefly, recent recalibration of the OUVS instrument has yielded an 
improved sensitivity for this early PV period, 52.5 counts s' 1 kR- '. An unproved 
algorithm was also developed that conserved photon counts, revealing a factor of two 
errors in previous NO nightglow data processing. Finally, the calculauon of the Stewart 
er al (1980) dark disk airglow intensity omitted a latitudinal weighting factor for area; 
its incorporation further reduced the observed dark disk btens ' 1 ^ u ^^f_ night : 
glow dark-disk average and peak intensities are now 400-460 ± 120 R and 
1.9 + 0.6 kR, respectively, for the 1978-1980 OUVS observational period. Individual 
orbit nightglow patches were seen to vary significantly in intensity and location on 
time-scales of an Earth day or less (Bougher, 1980; Stewart etal., 1980; Bougher etal 
1990a) The mean altitude distribution of the NO airglow was determined from limb 
profiles obtained near periapsis, with a 5 km vertical resolution. The altitude of the peak 
emission was found to be 115 ± 2 km (Gerard etal., 1981). 

The OUVS sensitivity was shown to degrade substantial after orbit insertion ( e- 
cember 1978). Subsequent monitoring of the nightglow beyond the first three dium 
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periods (1979-1981) incorporated this effect. The PV spacecraft periapsis altitude was 
too high during 1984- 1986 for the OUVS to obtain meaningful limb profiles of the solar 
minimum nightglow. However, several measurements during this period from apoapsis 
showed an average nightside bright patch intensity roughly a factor of three smaller than 
observed previously during 1979-1981. A substantial solar cycle variation in the night- 
glow emission is thus confirmed (Stewart, 1989, private communication). 

The No nightglow is particularly valuable in constraining the Venus thermospheric 
circulation in the absence of measured winds (Bougher, 1980; Stewart etal . , 1980; 
Bougher etal., 1990a). The verage brightness distribution map of Stewart etal. (1980) 
(revised in Bougher etal. (1990a)) can be viewed as an approximate map of the 
downward flux of N-atoms being supplied from the dayside. Correspondingly, the 
intensity of the resulting nightside emission is proportional to the dayside column 
production of N-atoms; i.e,, dayside average net column production must balance the 
nightside average column destruction of ground state N-atoms. Since the NO emission 
occurs well below the nightside homopause, vertical transport other than molecular 
diffusion must be important, either due to eddy mixing and/or large-scale winds. In 
addition, the local time position of the average peak emission (and bright patch intensity) 
provides evidence for asymmetries in the otherwise mean sub-solar-to- anti- solar ther- 
mospheric circulation (Stewart et al . , 1980; Bougher et al 1990a). Superimposed retro- 
grade zonal winds having a period of 4-6 days are also inferred from measured helium 
densities (Mayr etal., 1980; Mengel etal., 1989). 

The basics of the Venus dayside odd-nitrogen chemistry are summarized in Figure 18. 
The major processes are little different from those operating on the Earth. C0 2 , in fact, 
plays a similar role on Venus as 0 2 does in the lower thermosphere of the Earth; i.e., 
both species react rapidly with N( 2 Z)) but not with N( 4 S). The initial steps leading to 
the production of odd-nitrogen through N 2 bond-breaking are accomplished through the 
same processes. The principal sources are (1) the dissociation of N 2 by electron impact, 
(2) the dissociative ionization of N 2 by electrons, and (3) the predissociation of N 2 by 
80-100 nra photons. An effective branching ratio (/ = 0.5) is typically chosen for the 
partitioning of N( 2 Z)) and N( 4 S ) atoms that are produced (Cleary, 1986). The sequential 
reactions N ( 2 D ) + C0 2 producing NO followed by NO + N( 4 5) yielding N 2 are by far 
the most important in reducing the total dayside production of excited and ground state 
N-atoms to the net available for nightside transport. The reaction producing NO* 
followed by its dissociative recombination serve to redistribute the forms of atomic-N. 
A standard value of g = 0.75 for branching to N( 2 £>) is usually chosen (KJey etal., 
1977). In addition, partial recovery of N( 4 5' ) atoms is obtained by quenching of N( 2 Z) ) 
by O and CO above the N( 4 S*) peak. The precise value of the N( 2 D) + O quenching 
coefficient is very important for understanding the terrestrial odd-nitrogen chemistry (cf. 
Fesen etal., 1989). 

A preliminary study of the neutral and ion odd nitrogen chemistry was made by Rusch 
and Cravens (1979) soon after the first Pioneer Venus results were available. Their 
model predicted that ground state N( 4 S) atoms were the major neutral odd nitrogen 
species, with a peak concentration of about 2 x 10 7 cm " 3 near 130 km. The integrated 
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Fig. 18. Simplified Venus dayside odd-nitrogen chemical scheme. Taken from Bougher etal. (1989). 
column production of atomic nitrogen produced in their model (for 70° SZA) was 

4.5 x 10 9 cm' 2 s' '. Fox (1982) predicted a slightly larger maximum N( 4 S) density of 

4 x 10 7 cm - 3 . Later mapping studies and the nightside chemical-diffusive modeling of 
Stewart et al. (1980) showed that the originally reported dark disk nightglow intensity 
(780 R) required a downward flux of N-atoms of 1 x 10'° cm “ 2 s " *. This flux should 
be reduced by a factor of nearly two due to the recent downward revision of the 
intensities reported by Stewart etal. (1980). No correlation was found between the 
nightglow brightness and the 10.7 cm solar fluxes for the observational period. The 
nightside vertical diffusion model, and additional limb profile studies (Gerard et al., 
1981), furthermore, placed the peak volume emission rate altitude at 1 15 km using eddy 
diffusion of the same order as that used by von Zahn et al. (1980) for maintaining 
observed dayside densities. Thus, a determination of the altitude of the emission was 
found to be a sensitive indicator of the required nightside eddy diffusion profile (Stewart 
et al., 1980). Either global dynamics and/or small-scale mixing could be the underlying 
mechanism responsible. An adequate supply of dayside N-atoms appeared to be 
available for sustaining the observed nightglow. 

This conclusion was challenged by Krasnopol’sky (1983), who concluded, based on 
a simplified dayside atomic-N model, that virtually all the dayside produced N-atoms 
would be required on the nightside to account for the NO nightglow. His dayside total 
column production of both ground state (N( 4 S)) and excited (N( 2 Z>)) atoms was 
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1.1 x 10 ,o cm' 2 s“*. Virtually all of the N( 2 D) would have to be quenched to N( 4 S) 
rather than be converted to NO to sustain the required dayside to nightside flux. 
Presumably, this problem has been resolved with the reduction in the observed emission 
intensities. 

Gerard etai. (1988) re-examined this problem using a 1-D photochemical-diffusive 
model that calculated N^S), N( 2 Z)), and NO at selected SZA over the dayside. The 
daytime ionosphere and odd-nitrogen chemistry were coupled to obtain N( 4 S) profiles 
and local time distributions for comparison with PV measurements obtained with the 
neutral mass spectrometer (Kasprzak etaL, 1980; Hedin etai . , 1983). The average 
amount of total N-atoms produced on the dayside was estimated at 
1.3 x 10 lo cm~ 2 s" *. A maximum N( 4 S) concentration of about 8 x 10 7 cm" 2 was 
calculated at 124 km near LT * 17 : 00. Calculated N( 4 S) concentrations near 150 km 
were substantially less than those given by the empirical model of Hedin et al. (1983) 
(see Figure 19). However, a reasonable agreement with the original PV densities 
(Kasprzak et al., 1980) was obtained. Approximately 50% of the dayside source was 
required for nightward transport to maintain the originally reported NO nightglow 



Fig. 19. Vertical distribution of N( 4 S), N( 2 £>), and NO number densities as a function of altitude 
calculated at the equator for a local time of 16:00. Taken from Gerard et al. (1988). 
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intensities; the fraction would be smaller for the revised intensities. This estimate reflects 
a higher branching ratio for the (0, l)-<3 band resulting from NO recombination (0.32, 
after McCoy (1983)), which produces a required N-flux of 7 x 10 9 cm' 2 s~'. The 
model showed little sensitivity to the value of the efficiency of the N( 2 £>) production by 
N 2 dissociation (f = 0.5 to 0.7); also, the calculated N( 4 S) was found to be rather 
insenstive to the N( 2 Z))-0 quenching coefficient for values less than 
2 x 10 12 cm 3 s “ '. This is so because the reaction of N ( 2 D) with C0 2 yielding NO 
dominates all other N-atom destruction at the peak of the N( 4 S) density. 

A comprehensive treatment of this flux problem requires the use of a three-dimen- 
sional model coupling the hydrodynamical flow with photochemical production and loss 
of odd-nitrogen. The results of the Gerard etal. (1988) detailed photochemical model 
were subsequently used in a 3-D Venus thermospheric general circulation model 
(VTGCM) to parameterize the sources and sinks of atomic nitrogen in the Venus 
dayside thermosphere (Bougher et al., 1990a). The VTGCM (see Section 6.2), an adap- 
tation of the terrestrial NCAR TGCM, calculates global distributions of C0 2 , CO, and 
O consistent with the 3-D model by day-night temperature contrasts and corresponding 
large-scale winds. Calculations of the observed nightglow distribution and intensity 
provide a further means to validate the basic model circulation and thermospheric 
structure. The model N-atoms produced are transported to the nightside by the 
VTGCM circulation and destroyed primarily by radiative recombination with O atoms, 
and mutual destruction with NO. NO, N( 4 S), and N ( 2 D) are incorporated as minor 
species having no impact on the VTGCM major species, temperatures, or winds. NO 
and N( 2 Z>) are treated assuming photochemical equilibrium, while N( 4 S) was subject 
to diffusion and large-scale transport, in addition to chemistry. 

Results for solar maximum conditions indicate that the recently revised dark-disk 
average NO(0, 1) 6-band intensity at 198.0 nm (400-460 ± 120 R), based ons statisti- 
cally averaged PV OUVS measurements, can be reproduced with minor modifications 
in chemical rate coefficients and global mean eddy diffusion. The calculated average 
dayside production of total N-atoms is about 9.4 x 10 9 cm- 2 s"'; chemical losses 
resulting from NO production and subsequent N( 4 S) destruction reduce this total 
production by a factor of 2 to 4. The nightward hemispheric flux of N-atoms required 
to yield 340 to 580 R of dark-disk average airglow is ~ 2.5-3 x 10’ cm' 2 s' ', which 
is roughly 30% of the total dayside N-atoms produced. This is realized for VTGCM 
cases using standard rate coefficients and branching ratios (Gerard etal., 1988), or 
slightly reduced NO production through a lower limit N( 2 Z>) + C0 2 rate coefficient 
(3 x 10 l3 cm 3 s ') (Piper etal., 1987). A similar nightglow intensity could be 
simulated for somewhat weaker NO production coupled with enhanced global mean 
eddy diffusion (K,Z2x 10 7 cm 2 s' 1 ). A large N( 2 D) + O quenching coefficient, of the 
magnitude measured by Jusinski etal. (1988), is not required for adequate dayside 
N( 4 S) production yielding observed nightglow intensities. The VTGCM peak intensity 
of 1.2 kR, calculated at LT = 01 : 00-03 : 00 near the equator, is in good agreement with 
PV OUVS observations (see Figure 20(a)). The altitude of the nightside peak emission 
layer (113 km) is also in agreement with PV OUVS data, and occurs where chemical 
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Fig. 20. VTGCM NO(0, I ) &-band vertical intensity distributions: (a) a solar maximum simulation, corre- 
s ponding to OUVS data of 1978-1981. and (b) a solar minimum simulation, appropriate to the 1984-1986 
period. The vertical intensity distribution corresponds to the vertically integrated volume emission rate. 

Taken from Bougher etal. (1989). 


transport (global wind plus eddy diffusion) lifetimes are comparable. Use of enhanced 
global mean eddy diffusion (K ~ 2 x 10 7 cm 2 s “ 1 ) also yields atomic O and N densities 
that are consistent with observed PV values. Peak N densities of 6.5 to 9.2 x 10 7 cm - 
are obtained at LT = 17 : 00 at about 130 km. The simultaneous production of observed 
dayside N-densities, corresponding to Hedin etal. (1983) values, and observed NO 
nightglow intensities requires a slightly slower wind system (terminator horizontal winds 
~ 150 m s - ') than currently simulated by the VTGCM. 

Corresponding solar minimum calculations predict that the NO nightglow intensities 
should be very sensitive to the net dayside production of N-atoms. The factor of two 
decrease in the N 2 bond-breaking rate combined with the reduced VTGCM solar-dnven 
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winds causes a reduction of nearly a factor of 3 in the dark-disk intensities. Figure 20(b) 
illustrates the VTGCM nightglow integrated vertical intensity distribution for solar 
minimum conditions; comparison can be made with Figure 20(a). This solar cycle 

variation is similar to that observed by the PV OLVS. . . .. 

These VTGCM calculations confirm that the transport mechanism origin y 
proposed is correct. Since this airglow is a good tracer of the thermospheric circulation, 
it provides a means to study remotely the variation of the Venus winds. The nightglow 
also provides a sensitive constraint on the magnitude and distribution of the global 
winds in the absence of available measurements, much like helium (Mengel et a/., 1 9 b 
Furthermore, the Venus odd-nitrogen chemistry, since it is similar to that of the Earthy 
is also adversely impacted by the use of the large N ( 2 Z» + O quenching rate measured 
by Jusinski et at. (1988). It cannot be used for Venus and yield good dayside N-profiles 
plus NO nig htglow intensities resembling Pioneer Venus observations. On , 

present-day odd-nitrogen chemical schemes (e.g., Roble el a!., 1987 ; Fesen et a., ) 

and observed N ( 2 D) 520 nm airglow (Frederick and Rusch, 1977) cannot be reconciled 
with this large quenching rate. If N( 2 D) + 0 2 is the only important terrestrial NO source 
in the lower thermosphere, then the N( 2 £>) production rate must be increased to 
compensate for the large quenching rate. Recent measurements of this rate coefficient, 
however, have yielded values near 1 x 10" * 2 cm 3 s' \ a factor of 20 less than the 
Jusinski etal. (1988) rate (Piper, 1989; Miller etal, 1988). Further (Venus) progress 
requires future simultaneous measurements of N ( 2 Z», N ( 4 S), and NO in the Venus 
lower thermosphere (Gerard et a/., 1988). NO densities can also be derived from Venus 
nightside ion chemistry for comparison to VTGCM model values. 

3.2. 0 2 EMISSIONS 

The first measurements of 0 2 densities in the Venus atmosphere are those from the PV 
sounder probe gas chromatograph. Mixing ratios of 44 and 16 ppm were measured at 
52 and 42 km; it appears that the mixing ratio decreases with altitude (Oyama etal, 
1980). The gas chromatography experiments on Venera 13 and 14 measured an 0 2 
mixing ratio of 18 ppm for the altitude range 35 to 58 km (Mukhin et ai y 1983). The 
sounder probe mass spectrometer determined an upper limit to the mixing ratio below 

the cloud tops of about 30 ppm (Hoffman et at., 1980a). 

Above the cloud tops only upper limits have been measured. A ground-based search 
for absorption in the oxygen ‘A band’, the (0, 0) transition of the infrared atmospheric 
band system of0 2 (6‘i; -►A' 3 1" ), allowed Traub and 1 Carleton (1974) to placc an 
upper limit on the column abundance of 0 2 of 1 x 10' 6 ; Trauger and Lumne (1983) 

reduced this number fo 3 x 10 ~ 7 . . f 

Explaining these small 0 2 mixing ratios in the mesosphere has been a major goal of 
mesospheric photochemical modeling. Oxygen atoms produced in the photolysis of C0 2 

C0 2 + /iv-C0 + 0 (30) 

recombine very slowly by the process 
CO + O + M -»C0 2 + M . 
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At 300 K, the rate coefficient for reaction (31) is about 4.5 x 10 36 cm 6 s \ but the 
rate coefficient for the recombination of 0 atoms via 

o + o + m-o 2 + m < 32) 

is much larger, about 1 x 10 “ 32 cm 6 s ' 1 (Bauich et aL 1980). Photolysis of C0 2 would 
cause the entire thermosphere to dissociate in a few weeks without effective recombi- 
nation processes. The near-absence of 0 2 in the mesosphere indicates that some 
probably catalytic process occurring in the Venus atmosphere acts to reform C0 2 from 
the O produced in reaction (30) or from the 0 2 formed in reaction (2). On Mars, odd 
hydrogen compounds are implicated in the recombination of O and CO. A simple 
example of an odd-hydrogen catalytic cycle that effects recombination is 


H + 0 2 + M - H0 2 + M , ( 32a > 

H0 2 + O - OH + 0 2 , ( 32b > 

CO + OH - C0 2 + H , ( 32c ) 

with a net reaction of 

CO + o - co 2 ( 32d ) 


(Hunten and Me Elroy, 1970). On Venus, an extreme dearth of water vapor in the upper 
mesosphere precludes the presence of sufficiently large quantities of odd-hydrogen 
compounds. Compounds of chlorine have, however, been known to be present since 
HQ was identified in ground-based infrared spectra by Connes etal. (1967). The 
photochemical models of Sze and McElroy (1975), Krasnopol’sky and Parshev (1981, 
1983), and Yung and DeMore (1982) rely largely on catalytic cycles in volving Cl 


compounds, such as 

C1CO + 0 2 + M - C1C0 3 + M , 

(33a) 

CICO 3 + O -» Cl + C0 2 + 0 2 , 

(33b) 

Cl + CO + M - C1CO + M , 

(33c) 

with a net reaction of 

CO + O - C0 2 , 

(33d) 


to effectively recombine CO and O. 

The model of Krasnopol’sky and Parshev (1981, 1983), which is presented here in 
Figure 21, shows that 0 2 forms a layer 15 km thick near 90 km. The thermospheric 
density of 0 2 is consistent with a mixing ratio of 2 x 10 " 3 at 135 km. The upper 
boundary of the model of Yung and DeMore (1982) is at 1 10 km, where the mixing ratio 
of 0 2 is about lx 10' 3 . The model of Krasnopol’sky and Parshev has been critized 
by DeMore and Yung (1982) as including as a major step in destruction of 0 2 the 
reaction : 

C1CO + 0 2 -» C0 2 + CIO . < 34 ) 
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Fig. 21 Comput^ altitude profiles of the number densities of some neutral species in the Venus 
atmosphere. The C0 2 profile (dashed line) is a factor of 10« larger than shown. From Krasnopol’sky and 

Tomashova (1980). 


Laboratory studies have shown the reaction above does not occur, the C1CO and 0 2 
combining preferentially in the three-body recombination reaction (33d). Krasnopol’sky 
and Parshev (1983) have criticized the model of Yung and DeMore for using a constant 
H 2 0 mixing ratio of 1 x 10 “ 6 throughout the region above the clouds, when observa- 
tions show that immediately above the clouds the mixing ratio is close to 200 times that 
value. 

The chlorine and sulfur cycles have been found to be interconnected. DeMore et al. 
(1985) investigated the reaction of chlorine with S0 2 both with and without 0 2 and they 
also discussed the implications for the Yung and DeMore (1982) model. A new reservoir 
species for chlorine, S0 2 C1 2 (sulfuryl chloride) was proposed, potentially increasing the 
total abundance of chlorine compounds in the mesosphere. Chlorine was found to 
catalyze the oxidation of S0 2 to H 2 S0 4 , thus providing another catalytic cycle for 
destruction of 0 2 . 

All mesospheric models must explain the dayglow and nightglow intensities of 1.5 and 
1.2 MR observed for the 1.27 pm infrared atmospheric band, from ground-based 
measurements of Connes al. (1979). This band, which is the most intense feature in 
the Venus airglow, arises from the transition 

Oiin'Ag-, v' = 0)-» 0 2 ( 3 X ? ~ ; v" = 0) + A v( 1 7 pm) . (35) 
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The intensity of this feature is presumably not indicative of the ambient density of 
ground state 0 2 , but rather of its formation rate by recombination of oxygen atoms. In 
the dayglow, 0 2 (a x A t ) can also be produced by photolysis of ozone, but the near equality 
of the intensities of the day and night airglows shows that this source is minor for Venus. 
The model of Yung and DeMore (1982) required that the yield of 0 2 (a'd r ) be 30% in 
reaction (2) and 67% in the reactions: 

H0 2 + O — * OH + 0 2 , (36) 

a + Oj - a + o 2 , (37) 

CIO + O -♦ a + 0 2 . ( 38 ) 

Subsequently, the yields of 0 2 (a 1 A g ) and 0 2 (b l Z* ) have been measured for the 
reactions above and several other postulated odd-chlorine reactions and the values have 
been found to be very small, gnerall less than 4% (Choo and Leu, 1985; Leu and Yung, 
1987). Ali et al. (1986) have measured the rate for direct production of the ct x \ state 
in reaction (2) and found that the yield is only about 7%, but that in the presence of 
0 2 , the yield increases, indicating that the state is formed efficiently by energy transfer 
from a higher state. This conclusion is supported by an analysis of the oxygen bands 
in the terrestrial nightglow performed by Bates (1988b). Bates showed that the emission 
efficiency of the band systems originating in the A 3 Z * , A ' and cT~ states is small, 
indicating that the states are quenched efficiently, probably to the a 1 and b'l* 
states, producing emission in the atmospheric and infrared atmospheric bands. In order 
to account for the emission rates at 1.27 pm in the terrestrial nightglow it is necessary 
to assume a large efficiency ( — 0.75) for the sum of direct recombination of O-atoms 
and energy transfer from higher states, and an additional, as yet unidentified source (e.g., 
McDade etal., 1987; Lopez-Moreno et al., 1988; Lopez-Gonzales etal., 1989). 

Leu and Chung (1987) show that the total rate of C0 2 photolysis on Venus is about 
8 x 10' 2 cm ” 2 s ~ ', so the maximum emission rate of 1.27 pm emission would be 4 MR, 
if all the recombination reactions proceeded with a quantum yield of 1.0. The measured 
yields above imply an intensity of 0.05 MR of 1.27 pm emission from the chlorine 
reactions, whereas the observed value is 1.2— 1.5 MR. Leu and Yung argue, therefore, 
that the 1.27 p emission cannot arise from recombination of O atoms formed in 
photolysis of C0 2 , and suggest an alternative source of O-atoms: photolysis of S0 2 
introduced into the upper mesosphere by intermittent injections from the region near 
the cloud tops. The discrepancy for Venus is, however, rather more severe than the 
terrestrial nightglow problem, and seems to require a reconsideration of the basic theory 
of production of the 1.27 pm emission. 

Massie etal. (1983) modeled the thermospheric densities of 0 2 , but at the lower 
boundary of the model, 100 km, the 0 2 mixing ratio was constr ain t to the value 
indicated in the model of Yung and DeMore (1982). In the thermosphere, between 100 
and 110 km photolysis of C0 2 is balanced by reaction (2). The computed mixing ratio 
of 0 2 at 1 35 km is about 1 x 10 “ 3 . 

Thermospheric 0 2 densities theoretically can be determined from the chemistry of 
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atomic C. The major loss process for C in the thermospheres of Venus and Mars is the 
reaction 

C + 0 2 -C0 + 0 (23) 

(McElroy and McConnell, 1971; Krasnopol’sky, 1982; Fox, 1982b). Although no 
measurements of C densities are available, since the predicted densities were just below 
the sensitivity of the PV ONMS, the PV orbiter ion-mass spectrometer measured 
densities of C* (e.g., Taylor et al , 1980). Fox (1982b) used the measured C* densities 
to constrain the 0 2 mixing ratio; values of 10” 4 or less were derived. Paxton (1983) 
used the intensities of atomic carbon emissions measured by the PV OUVS (see 
Section 3.4) to constrain the 0 2 densities, and he derived an 0 2 mixing ratio of 
3 x 10 ” 3 . As pointed out in Section 3.4, all of the models above are based on outdated 
photochemistry, and should be revised. 

An outstanding success of the Soviet Venera 9 and 10 missions was the measurement 
of the visible nightglow of Venus in October 1975. The strongest features in the nightglow 
spectrum, shown in Figure 22, were identified by Lawrence et al. (1977) as the v' = 0 



Fig. 22. Visible night airglow spectrum measured by Venera 9 and 10. The dashed line is the instrument 
sensitivity. Taken from Krasnopofsky and Tomashova (1980). 

progression of the 0 2 (c *1* -► X*Z ~ ) Herzberg II band system; an average intensity of 
2.7 kR was assigned to the band system (Krasnopofsky, 1981). Three more band 
systems of 0 2 were later identified in the Venus spectrum: the Herzberg I 
{A y Z: -Jr 3 I”) bands, the Chamberlain (A ,y A u -^a x A t ) bands and the Slanger 
(c 1 T” -► a l A b ) bands. The measured intensities of the Herzberg I and the Chamberlain 
bands were 140 and 200 R, respectively. All of these band systems are formed in 
recombination of atomic oxygen produced by photolysis of C0 2 on the dayside of the 
planet and transported to the nightside (Lawrence e/a/., 1977; Slanger, 1978). The 
maximum of the airglow layer is located near 100 km, the altitude of the maximum O 
density, and is about 15 km thick. Variations of about 10 km in the altitude of the 
maximum were observed. The emission shows a maximum at about 8-28° S latitude; 
diurnal variations of the nightglow intensities, shown in Figures 23, exhibit a maximum 
about 00 : 30 loca time, consistent with transport of O produced on the dayside by C0 2 
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Fig. 23. Diurnal variation of the visible night airglow as measured by Venera 9 and 10. The dashed line 
is the symmetric continuation of the left part of the curve. From Krasnopol’slcy and Tomashova (1980). 


photolysis across the terminator and subsidence at the anti-solar point. A constraint 
similar to that placed on the Venus solar minimum thermospheric structure and circu- 
lation by the measured NO nightglow emissions is given by using the Herzberg II 
nightglow observations of the Venera 9 and 10 instruments. The VTGCM could be used 
to examine this visible airglow, provided that total temperature is incorporated explicitly 
and that its model lower boundary is extended from 105 to ^90 km (Bouger etai , 
1988b, 1989). Such a study could yield information about the thermospheric circulation 
in the vicinity of 100 km. 

Subsequent to the identification of the Herzberg II bands in the Venus nightglow, the 
band system was identified in the terrestrial nightglow, where the major molecular 
features are the Herzberg I bands (Slanger and Huestis, 1981). The dominance of 
different band systems in the terrestrial and Venus nightglows has led to many 
laboratory and modeling studies in an attempt to explain the differences. Slanger (1978) 
and Kenner etai (1979) showed that the presence of C0 2 enhances emission in the 
Herzberg II bands, possibly due to its higher efficiency as a third body. Krasnopol’sky 
ex al (1976) proposed that the non-appearance of the Herzberg II bands in the Martian 
nightglow, as measured by the Mars 5 space probe, suggests that they are suppressed 
by efficient quenching by 0 2 . Subsequently, the rate coefficient for quenching of the 
0 2 (c , r- ; v = 0) state by 0 2 has been measured by Kenner and Ogryzlo (1983); the 
reported value is 3 x 10" 14 cm 3 s“ \ whereas that for quenching by O is a factor of 
200 larger. From studies of the dependences of the band systems on [O], [0 2 ], and [M], 
Kenner et al. (1979) found that the ratio of the intensity of the Herzberg II bands to that 
of the Herzberg I bands is proportional to [M)/[0]. They suggested that the Herzberg II 
bands dominate on Venus because the maximum O recombination takes place in a 
region of higher total density and lower atomic oxygen density. Although this would 
imply an even greater dominance of the Herzberg II bands on Mars, where the total 
density in the region of maximum O recombination is about 10 times that in the Venus 
atmosphere and 10 3 times that in the terrestrial atmosphere (Krasnoporsky, 1981), it 
is possible that quenching renders the total emission rate small. Indeed, from calcu- 
lations of the transition probabilities of the 0 2 band systems (Bates, 1988a; values 
corrected in Bates, 1989), Bates (1988b) has shown that the efficiency of emission 
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relative to quenching in both the Herzberg I and II band systems in the terrestrial 
nightglow is on the order of 1-4%. 

From an analysis of terrestrial and Venusian airglow profiles, Krasnopoi'sky (1986b) 
proposed that the weakly bound 5 /7 f state, produced efficiently in three-body recombi- 
nation of O atoms, is a common precursor for the five metastable states of 
0 2 (A,A', c, a, b), the lower states being produced by quenching or energy transfer to 
another 0 2 molecule. This suggestion was based on calculations of the efficiency of 
production of various states of 0 2 by Wraight (1982) and Smith (1984) that showed that 
the 5 II g state is produced with about a 66% efficiency. Bates (1988b) proposed that the 
fraction be reduced to 50% to account for the smaller potential well depth of the i Il f 
state, about 0.14 eV as computed by Partridge el al. (1990, private communication to 
Bates), compared to the value 0.23 eV adopted by Wraight from the calculations of 
Saxon and Liu (1977). He proposed the fractions of associations into the seven lowest 
electronic states of 0 2 shown in Table VIII. The efficiency of production of the A state 


TABLE VIII 


Fraction of direct associations (reaction 39) into 
various stales of 0 2 . (From Bates, 1988b.) 


State 

Fraction 

x’z; 

0.12 


0.07 

b'z; 

0.03 

c'z; 

0.04 


0.18 

a 3 zj 

0.06 


0.50 


is only about 6%, 50% larger than that of the c state. Bates showed that at temperatures 
near 200 K redissociation of the i IJ f state is efficient and much larger than the potential 
rate for energy transfer to any of the lower states. Furthermore, because of this redisso- 
ciation, the measured value of the three-body recombination rate coefficient underesti- 
mates the actual recombination rate. Temperatures are slightly lower in the Venus 
atmosphere near 100 km, however, so redissociation may not be as important there, and 
the dissociation energy itself is still somewhat uncertain. 

The excited states of 0 2 are relevant also for production of 0( l S) and subsequent 
radiation in the atomic oxygen ‘green line’ at 5577 A in the terrestrial nightglow. It is 
generally agreed that the production of 0(‘S) proceeds by the two-step mechanism 
proposed by Barth (1964): 


(39) 
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This suggests that the precursor 1 P ~ . t for auenching by O and its 

was ruled out based on the large required rate coeffic y eweUyn % e t ai, 1980). 

altitude prohle (e.g.. Slanger , 19 ™ ‘J ^r, the green line 
Krasnopol’sky (1981) showed dial, . observed value, but that the 

emission would be two orders of ma S a '\ u e g precursor did not contradict any 
identification of the vibrationally excit csa terrestrial nightglow consists of 

available data. Indeed, the Herzberg I n system » the terresmal ^ ^ 

bands originating in vibrational leve s w v ^ ’ hi a ssumpUon, Krasnopol’sky 

originating in v * ^f^l^g^^^io^Sj-’Jtheqiienching of thee slate by O: 
derived a rate coefficient of 2.5 x 1U 

o 2 (c , r u -) + o-»o J + o( , s), 

a ..due of 5.9 x 10-on* »as measured by Kenner and Ogryrlo (1983) for « 

reacuon (41b) 

0 2 (c'l;) + O- products . 

It is probable, however that the r * lc ^* C> “* S “brational distributions in the 
depend upon the vibrahonal ^tain. Moreover, Krasnopofsky's 

atmosphere and sometimes in the labora O'. (1984) and was 
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Huestis, 1981) with computed 
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limits for the transfer rate coefficients for the A state, which is marginally consistent, 
and for the A ’ state, which is inconsistent, with the terrestrial required rate coefficients. 
He does note, however, that the rate coefficients may vary with vibrational quantum 
number, so the constraints ar not rigid. The upper limit on the energy transfer rate 
coefficients for the s II g state from Venus data is inconsistent with terrestrial requirement, 
but there is no positive evidence in its favor. 

3.12. 0 2 + EMISSIONS 

The O; second negative band system arises from the transition 0 2 (/f 2 J7„-*A' 2 /7 r ). 

H*) state is produced in the Venus ionosphere mostly in fluorescent 
scattering of solar radiation by ground state 0 2 *, since the ambient densities of 0 2 in 
the Venus thermosphere are too low for photoionization and electron impart ionization 
to be significant sources. In order to determine the role of fluorescent scattering in 
producing vibrationally excited 0 2 \ Fox (1985) computed the intensities of the bands 
in the second negative and first negative (b 4 Z~ ->a 4 /7J band systems. The second 
negative bands are weak and spread over a large wavelength range. The intensities in 
50 A intervals are shown in Figure 24 for two assumptions about the (unknown) rate 
coefficient for the vibrational quenching of 0 2 + {X 2 T1 S ; r) by O: 

(r) + O -» 0 2 + (t> - 1) + O . 


>s 




Wavelength (A) 


Fig. 24. Computed intensities of the second negative band system, for a high solar activity (HSA) model 
summed in 50 A intervals. The solid line is for k„-lx I0‘ ,0 cm J s* '. The dotted line is for 

*43 - 6 X 10" 10 cm 3 s’*. 
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4. Hot Atom Comas and Escape 
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a typical collision cross section is about 3 x 10' 15 cm' 2 , m pracuce the exobase ^s 
located near the altitude above which the column density is about 3.3 x 
the high solar activity (F 10 , = 200) VTS3 model the exobase alutude is ocas 210^ at 
noon and 154 km at midnight. The VIRA model for 16 N and F 107 = 150 places the 
exobase near 201 km at noon and 161 km at 22 : 00 hr local time. In the BNMS model 
which applies to 08 : 30 hr local lime, the exobase is near 191 km. In the region of 
pre-dawn bulge in the densities of H and He, the column density of H 2 is important in 

determining the altitude of the exobase (Kumar et al., 1983). 

Whether the trajectory of a particle traveling upward from the exobase is ballistic 
(bound) or escaping (free) is determined by the total energy £, which is the sum of its 
kin etic and potential energies: 


£ = - mi; 2 + 


f mGH 

J r 2 


dr, 


(44) 


where m and u are the mass and velocity of the particle, G is the gravitational constant, 
M is the mass of the planet, and r is the distance from the center of the planet, the 
subscript c refers to the critical level or exobase. Expression (44) reduces to 


£ - ^mi> 2 - mg c r e , 


(45) 


where g c is the gravitational acceleration at the exobase. If the total energy is negauve, 
the particle is bound; if the total energy is positive the particle may escape. The escape 
velocity, v w , is then defined by the condition t> e , e = (2gr) ,/2 . For Venus, the escape 

velocity is about 10.2 km s" '. . .. 

In the Jeans escape process, particles with velocities greater than the escape velocity 
in the high-energy tail of the Maxwellian distribution may escape if they are onented in 
the upward hemisphere (above the horizon). The escape flux, (pj is given by 


<P j 


2 Jn 


(1 + A e )exp(-<U, 


(46) 


where 


GMm mgr (47) 

X “ rkT “ kT 

X is the gravitational potential energy in units of kT and u is the modal velocity of a gas 
in thermal equilibrium at temperature T, 


u = (2 kT!m) X12 . 


(48) 


Sometimes a correction factor is applied to the expression for the escape flux to account 
for the suppression of the tail of the distribution due to the escape of the en^getic 
particles (Chamberlain, 1963; see also Hunten, 1982). Application of Equation (46) to 
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the escape of H from Venus, for the VI RA model with an exobase near 200 km, gives 
a value for A of about 1 1 and a thermal escape flux of about 22 H atoms cm “ 2 s " 1 at 
noon; at midnight the escape flux is several orders of magnitude smaller, due to the low 
nightside temperatures and the resulting large value of X. These fluxes are insignificant 
compared to the actual escape fluxes that have been inferred from models of the hot 
H corona. 

Hunten (1973a) was the first to recognize that the escape rate of a light species from 
a planetary atmosphere may be controlled by diffusion of the species to the exobase, 
rather than by the escape process itself. The escape flux of H from Venus is limited by 
the rate of transport of H through the middle atmosphere to the upper atmosphere; the 
limiting upward flux, <p t of a species / with mixing ratio f, can be estimated as 



(49) 


where H a is the average scale height of the atmosphere and b t is the binary collision 
parameter (Hunten, 1973a, b; Hunten and Donahue, 1976). b, is defined by b, =* D,n T , 
where D, is the molecular diffusion coefficient, evaluated as an appropriate average of 
the binary diffusion coefficients D v of the species / diffusing through the major 
atmospheric constituents j. D, is usually computed as 


-L* V n ‘! n r 

D, J+i Dy 


(50) 


The expression (49) above is usually evaluated at the homopause, with the mixing ratio 
taken from a suitable altitude in the middle atmosphere, but above the cold trap. The 
limiting flux obtains if and only if the mixing ratio is constant with altitude. The effect 
of photochemistry can be eliminated if all chemical forms of the species are counted in 
the calculation of /,. H is found in the mesosphere in the form of HC1, H 2 0, and H 2 
an din the thermosphere mostly as H. In order to evaluate the limiting flux of H we must 
first determine the mixing ratios of H-containing compounds in the middle atmosphere 
and at the homopause. The abundances of H, H 2 , and H 2 0 will be discussed below. 


4.2. H AND H 2 DENSITIES 

No in situ measurements of H densities in the thermosphere are available, but H * 
densities were reported from PV ion mass spectrometer data (e.g., Taylor et al„ 1980). 
Brinton etal. (1980) showed that, in the altitude regin where H* is in photochemical 
equilibrium, the major source of H + is the nearly thermoneutral charge transfer of O * 
to H 


0 + + H-H* +0, 
while loss is via the reverse reaction 

H * + O-H + CT 


(51a) 

(51b) 


SPAC 191/59 



60 


J. L. FOX AND S. W. BOUCHER 


and by reaction with C0 2 

H + + C0 2 -+ HCO * + O . (52) 

Thus the steady-state density of H can oe computed from ONMS and OIMS data using 
the expression 


[H)-m 


[O 


U([0] 

nv k ila 


+ ~ [C0 2 

la 


Bnnton et al. found that photochemical equilibrium is a valid approximation below 
200 km on the dayside and below 170 km on the nightside. A large variation with solar 
hourangle was found, with maximum densities occurring in the pre-dawn sector near 
04 : 00 hr local time. Figure 25 shows the H densities derived by Brinton et al. at 165 km 



HOUR ANGLE 


Fig. .5. Diurnal variation of atomic hydrogen concentration in the Venus thermosphere, derived from 
« rm, measurements of ion and neutral composition. Data were obtained near 165 km altitude on 25 orbits 
of The Pioneer Venus spacecraft between December 1978 and July 1979. These values are based on 
non-normaltzed ONMS data and should be increased by a factor of about 1.63 (see text). Taken from 

Brinton et al. (1980). 


as a function of hour angle. Since this study was done, a renormalization of the ONMS 
data by a factor of 1.63 has been suggested (Hedin et al., 1983) (see Section 2), so values 
o [HJ presented by Brinton et at. should be multiplied by this factor also. Taylor et at. 
(1984, 1985), obtained similar results in a later analysis of three years of PV data. The 
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The dayside densities are significantly larger and the nightside densities are smaller than 
the values derived by Brinton et al. (1980) from Pioneer Venus data. The large diurnal 
variation in densities shows that the spherically-symmetric radiative transfer model used 
by Anderson is not valid, especially for the nightside. In addition, Anderson placed the 
critical level at 250 km, significantly higher than Pioneer Venus models indicate, 
although, as we have seen, the H densities do not vary rapidly with altitude. 

The Manner 10 spectrohotometer detected a strong signal at La but the intensity was 
less than at the time of the Mariner 5 encounter (Broadfoot et al., 1974). Takacs et al. 
(1980) used a spherically-symmetric radiative transfer model to analyze the emissions 
on the bright limb. The derived values for the temperatures and densities were similar 
to those of Anderson (1976): 275 ± 50 K and 1.5 x 10 5 cm for the temperature and 
density of the dayside cold component, respectively; and 1250 + 100 K and 500 + 100 
for the temperatures and density of the dayside hot component, respectively. Takacs 
et al. found that a spherically-symmetric model was inadequate to analyze the dark disk 
data. 

In addition to the visible spectrometers, Venera 9 and 10 were equipped with La 
photometers to measure the intensities at 1216 A. Bertaux et al. (1978) analyzed the data 
obtained using a spherically-symmetric, isothermal radiative transfer model. They found 
that the data could be fit by a two-temperature model, but the fit was not better than 
that of a one-temperature model at 500 K with an exobase density of 1.5 x 10 4 cm' 3 , 
values which are bracketed by those of the hot and cold components from the Mariner 5 
and 10 analyses. The spacecraft also carried resonance-absorption cells filled with 
hydrogen and deuterium that were heated by tungsten filaments to produce H and D 
atoms, respectively. The reduction in the intensity of the planetary La line by the 
absorption cells would provide an indication of the line width and therefore of the 
temperature. No absorption was observed for the D 2 cell. The H 2 ceU data showed that 
the temperature increases abruptly above about 3000 km, near the altitude above which 
the Manner analyses showed the hot component to dominate. 

The Venera 1 1 and 12 ten-channel UV photometers were also sensitive at La (Kurt 
et al., 1979). Bertaux et al. (1981) reported sub-solar zenith intensities of 38-42 kR on 
the disk, similar to the values reported by Anderson (1976) from Mariner 5 data. 
Anderson showed, however, that the bright disk intensities are due mostly to Rayleigh 
scattering by C0 2 . Since the Venera photometer aperture was not parallel to the limb, 
the altitude resolution in the limb scan was limited. Analysis showed that the data could 
be fit by a two-temperature model with temperatures of 400 and 700 K. 

Reports of Pioneer Venus OUVS La limb scan observations have been presented by 
Paxton et al. (1985, 1988). Paxton et al. (1988) used a radiative transfer model that 
allowed for variations of temperature and density with altitude, latitude, and local time 
to analyze 20 orbits covering the first three years of operation of the PV orbiter. The 
inferred density of H at the exobase averaged 6 x 10 4 cm' 3 ; this compares favorably 
with the daytime average value of about 6.5 x 10 4 cm " 3 determined from the analysis 
of Brinton et al. (1980). They attributed the factor of 2-3 discrepancy with the Mariner 5 
and 10 analyses to the use of a more sophisticated radiative transfer model. The column 
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density of H above 1 10 km, the altitude below which C0 2 strongly absorbs Lot, averages 
3.6 + 1 x 10' 3 cm" 2 . 

4.3. D/H RATIO 

Interest in determining the D/H ratio in the Venus atmosphere began when the 2-scale 
height distribution of H was first measured by Mariner 5 (Barth et al . , 1967; see 
Section 4.2). The interpretation that the inner component was due to D at the same 
temperature as the outer H component required a D/H ratio of about 10% (McElroy 
and Hunten, 1969; Wallace, 1969). Although the derived ratio, which was enhanced 
over the terrestrial value by a factor of about 10 3 , became meaningless after this 
hypothesis was shown to be invalid, the idea of some enrichment of D due to differential 
escape of H and D remained a possibility. 

The PV OIMS detected an ion with a mass of 2, which could be either H 2 of D * 
(e.g., Taylor etal., 1980). Originally this ion was believed to be H 2 , with the H 2 * 
produced mainly by direct photo- and electron-impact ionization of H 2 . H 2 would be 
destroyed by reactions with O * (reaction 54) and with C0 2 : 

H 2 + C0 2 -* C0 2 H * + H . (55) 

Kumar et al. ( 198 1 ) derived an abundance of 1 0 ppm H 2 at the homopause from a model 
of the chemistry of H 2 + . McElroy etal. (1982), Cravens etal. (1983), and Rodriguez 
etal. (1984) have emphasized that a large mixing ratio of H 2 provides a significant 
source of H through the reactions (54) and (55) above, which cannot be balanced by 
the known loss mechanisms and is difficult to reconcile with the measured abundances 

of H. 

McElroy et al. (1982) suggested that the mass-2 ion measured by the PV OIMS was 
D * rather than H-T , the atmosphere having become enriched in D relative to H by 
preferential escape of H. Donahue et al. (1982) were able to measure the D/H ratio in 
the lower atmosphere with the PV LNMS, by comparing the H 2 0 and HDO peaks 
during the time that the inlet to the mass spectrometer was clogged by sulfuric acid from 
cloud droplets. A value of 1.6 x 10 2 was reported. 

Hartle and Taylor (1983) identified the mass-2 ion as D + by comparing the scale 
height of the mass-2 ion to that of H + in the photochemical equilibrium region and the 
diurnal variation of the mass-2 ion to that of H + . They reasoned that the altitude profile 
of the ratio D*/H + should vary as exp(dz/ff(l)), where H( 1) is the scale height of a 
mass-1 species, but the ratio of H 2 + /H* should vary as [CT ]"' exp(dz///(l)). The 
observed altitude variation of the ratio supported the identification of the species as D * , 
as did the observation that both the mass-2 ion and H * exhibit bulges in the pre-dawn 
sector, with densities an order of magnitude larger near 04 : 00 hr local time than during 
the day. From the measured ratio D*/H"" Hartle and Taylor derived a D/H ratio of 
1.7 x 10" 2 between 155 and 160km, which, when projected down to the presumed 
homopause at 132 km yields a value of 2.2 ± 0.6 x 10" 2 , in fairly good agreement with 

the measured value of Donahue etal . (1982). 

Kumar and Taylor (1985) modeled the pre-dawn bulge ionosphere, where the maxi- 
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mum densities of H * and D * (as well as those of light neutral species) are found. By 
analyzing two orbits of PV data, they obtained additional strong evidence that the 
mass-2 ion was D + . The reaction of 0" with H 2 (reaction (54)) would strongly 
suppress the 0 * densities; a depletion in O + at times of enhanced mass-2 ion was not 
observed. Kumar and Taylor derived D/H ratios of 2.5 and 1.4% for PV orbits 117 and 
120, respectively, and placed an upper limit on the H 2 abundance at the homopause of 
0.1 ppm. 

Using the high-resolution mode if the IUE satellite, Bertaux and Clarke (1989a) 
searched for the deuterium La line, which is located 0.33 A shortward of H La, but 
found no evidence for it. An upper limit was determined for its intensity of 300 R, 
compared to 21 kR measured for the H La line. From the non-detection, they placed 
an upper limit on the D/H ratio at 100 km of 3.6 + 1.5 x 10“ 3 , a factor of more than 
4 lower than the values estimated from the LNMS data by Donahue etai (1982) and 
Hartle and Taylor (1983). Bertaux and Clarke suggested that the LNMS result for HDO 
may indicate a chemical fractionation process, and that the mass-2 ion may indeed be 
H 2 + , thus requiring the presence of 10 ppm H 2 . They did not, however, account for the 
lack of suppression of O + in the bulge region, or for the large source of H atoms implied 
by such a large abundance of H 2 . Donahue (1989) pointed out the inconsistency in their 
inference of 10 ppm H 2 a the homopause and the 0.7 ppm H required by the H La data. 
He also challenged their conclusion that a D/H ratio of 1 .6 x 10 " 2 would lead to 2.5 kR 
of D La, as ignoring the upward flux of D due to the global circulation, which will 
produce a different verticle distribution from diffusive equilibrium. Donahue found that 
a predicted intensity of 1 .0 + 0.4 kR is more appropriate and that the upper limit implied 
by 300 kRofD Lai$about5 x 10“ 3 . Bertaux and Clarke (1989b) suggest that the larger 
D/H ratio measured by Donahue etai. (1982) should produce observable limb 
brightening of the D La emission and that a search of the data collected by the 
PV OUVS should reveal such a limb brightening. 

According to Yung and DeMore (1982), the ratio of H 2 to H at the homopause is 
fixed by the photochemistry of H 2 0 in the mesosphere and should have a value of about 
10%. For the VI RA standard noon model, an extrapolation of the H densities down- 
ward to 130 km, the approximate altitude of the H homopause, yields a mixing ratio of 
about 3 ppm, although this extrapolation, like that for He, is extremely sensitive to the 
assumed value of the eddy fissuion coefficient and should not be considered accurate 
to more than a factor of 2-3. The mixing ratio of H 2 at the horaopause would, therefore, 
be about 0.3 ppm. 

H 2 has been measured below the cloud tops by the PV large probe gas chromatograph 
(LGC). The mixing ratio appears to decrease toward the surface; the reported values 
at 52,42, and 22 km are 200, 70, and 10 ppm. The gas chromatographs on Veneras 13 
and 14 had a threshold sensitivity for detection of H 2 of 0.8 ppm. A value of 25 ppm 
was measured within the clouds for the altitude range 49-58 km (Mukhin etai , 1983). 

The mesospheric photochemical models of Yung and DeMore (1982) and Yatteau 
(1983) required much lower mixing ratios of H 2 above the cloud tops than were 
measured in or below the clouds or derived from the chemistry of H 2 * at the homopause 
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by Kumar et al. (1981). The catalytic processes that recombine CO and O in the lower 
atmosphere of Venus are assumed to involve odd-chlorine compounds, rather than odd 
hydrogen as in the Martian atmosphere, where the abundance of water vapor is larger 
(see Section 3.11). H 2 renders Cl ineffective as a catalyst because the reaction 

H 2 + Cl -► HC1 + H (56) 

produces HCI, which is a relatively unreactive ‘reservoir’ species for odd chlorine. Yung 
and DeMore presented three models that differed mainly in the mixing ratio of H 2 ; 
values of 20 ppm, 0.5 ppm, and 1 x 10” 7 ppm were tested. Only the extremely 
hydrogen-deficient model could account for the observed nighttime depletion of CO 
near the cloud tops reported by Schloerb et al. (1980), Wilson et al. (1981), and Clancy 
etal. (1981) from microwave observations. 

Other problems with these models have, however, become apparent. First of all, as 
Krasnopol sky and Parshev (1983) have noted, a constant mixing ratio of 1 ppm H 2 0 
was used over the entire mesosphere, in contrast to the PV orbiter infrared radiometer 
(OIR) measurements, which show that, near the cloud tops at least, the mesosphere is 
not as dry as assumed. Second, some of the 0 2 production reactions involving odd 
chlorine that were proposed to lead to production of0 2 (a'^), explaining the observed 
strong 1 .27 ft emission have been shown, by Yung and coworkers, to produce 0 2 (a'A g ) 
very inefficiently (see Section 3.1 1). It also appears that the diurnal variation of CO in 
the Venus mesosphere is more complicated than early microwave observations 
indicated. Clancy and Muhleman (1985a) showed that the CO mixing ratio in the 
altitude range 80-90 km is depleted by a factor of 2-4, but at 95 km is enhanced by a 
factor of 2-4 on the nightside relative to the dayside. The maximum in the 80-90 km 
region occurs in the early morning rather than at mid-day. Clancy and Muhleman 
(1985b) have, however, proposed that the midnight bulge above 90 km is a result of the 
sub-solar to anti-solar circulation rather than a photochemical effect. Nonetheless, it 
appears that the photochemical processes that voern the recombination of O and CO 
and the destruction of 0 2 are not yet well understood. 

4.4. H 2 0 densities 

The abundance of water is unknown in the thermosphere of Venus and controversial 
in the lower atmosphere. Below the cloud tops, the Venera 9 and 10 narrow-band 
photometry measurements imply an H 2 0 mixing ratio of 300 ppm, accurate to within 
a factor of 2 in the altitude range 20-40 km (Ustinov and Moroz, 1978). The Venera 1 1 
and 12 descent probes carried scanning infrared spectrophotometers that recorded the 
spectrum of the Venus datytime sky. Synthetic spectra constructed by Moroz et al. 
(1980) indicated that the H 2 0 abundance is about 20 ppm near the surface increasing 
to 200 ppm in the clouds. From measurements made by the net flux radiometers on the 
four PV probes, Revercomb et al. (1985) have reported that the water abundance below 
the clouds varies with latitude, with mixing ratio of 20 to 50 ppm at 60° increasing to 
about 500 ppm at the equator. The gas analyzers aboard the Soviet Vega 1 and 2 
spacecraft indicated that the water vapor is concentrated in a 30 km layer centered at 
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40 km, with a maximum mixing ratio on the order of 1000 ppm at about 50 km, 
decreasing to about 20 ppm near the surface (Surkov et al, 1987). 

For the atmosphere above the cloud tops, ground-based infrared absorption measure- 
ments have indicated an H 2 0 mixing ratio of 1 ppm (Fink et al, 1972; Barker, 1975). 
Subsequently space probes have yielded conflicitng results. In the cloud layer, from 4 
to 58 km, the gas chromatography experiments on Veneras 13 and 14 yielded a value 
of 700 ± 300 ppm (Mukhin et al, 1983). The mixing ratios measured by the Vega gas 
analyzers appear to decrease over the altitude range 60 to 75 km, with values in the range 
1-100 ppm (Surkov et al., 1987). The error bars, however, are more than a factor of 10 
near 75 km. Above about 65 km the PV OIR showed a large diurnal variation in the 
water vapor abundance from about 100 ppm in the early afternoon to less than the 
detection threshold of 6 ppm at night (Schofield et al., 1982). The abundance near the 
cloud tops is probably determined by the saturated vapor pressure of H 2 0 over sulfuric 
acid and should decrease with increasing altitude. Krasnopol’sky (1985) has computed 
a sulfuric acid ‘trap function’, which relates the water vapor and H mixing ratios at 
90 km to that in the lower atmosphere. For an assumed lower atmosphere mixing ratio 
of 200 ppm, the mixng ratio of H 2 0 is about 0.5 ppm at 90 km and that of all 
H -containing compounds is about 3 ppm. 

4.5. Diffusion limit for escape of H 

On Venus, the atmosphere above the cloud tops, the mesosophere, has generally been 
chosen as the region in which to evaluate the ixing ratio of H-containing compounds, 
in order to calculate the diffusion limiting upward flux of H (Hunten and Donahue, 
1976; Walker, 1977). The mixing ratio of HC1 has been determined from ground-based 
infrared absorption measurements to be about 0.6 ppm (Connes etal., 1967). In view 
of the influence of the clouds, care must be exercised in choosing the altitude range for 
determining the abundance of water vapor. Krasnopol’sky (1985) asserts that sulfuric 
acid affects the H 2 0 mixing ratio up to about 75 km, so the abundance at the cloud tops 
is not an appmropriate value to use in computing the diffusion limited flux. The 
abundance of H 2 in the clouds has been reported as 25 ppm by Mukhin etal. (1983) 
from Venera 13 and 14 gas chromatography, but mesospheric models cannot accommo- 
date such large mixing ratios (Yung and DeMore, 1982 ; Yatteau, 1983). If it is assumed 
that the abundance of H 2 above the cloud tops is negligible, and that the abundance 
of H 2 0 above 75 km is 1 ppm, then the total abundance of H is 2.6 ppm. This is good 
agreement with the H-mixing ratio extrapolated to the H homopause at about 130 km 
of about 3 ppm, but the latter value is accurate to not better than a factor of 2. 

In order to compute the flux through the homopause, the binary collision parameter 
{b = Dn r ) must be evaluated. We assume here that the major form of H at the 
homopause is atomic hydrogen. The diffusion coefficient of H through the atmosphere 
can be computed as a suitable average of the binary diffusion coefficients for H through 
each of the major atmospheris constituents (cf. Chamberlain and Hunten, 1987). The 
values of the diffusion coefficients assumed are shown in Table IX. The diffusion 
coefficients of H through N 2 is taken from Banks and Kockarts (1973), and is taken 
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TABLE IX 

Parameters A and s used to compute diffusion coefficients for H through the 
major gases in the Venus thermosphere. The diffusion coefficient in cm 2 s ' 1 
is given by D - AT*jn T% if the temperature T is in K and the total density n T 

is in cm " \ 


System 

A (10‘ 7 cm" 1 $-' K~‘) 

s 

H-O* 

7.25 

0.71 

H-COj b 

3.87 

0.711 

H-Nj c 

4.87 

0.698 

H-CO* 

4.87 

0.698 


• Cooper tt ai . (1984). 

b Obtained by scaling the values for H-0 2 from Banks and Kockarts (1973) 
for the larger mass and radius of C0 2 . 
c Banks and Kockarts (1973). 
d Assumed the same as for H-N 2 . 

to be the same as that for H-CO. The H-C0 2 diffusion coeficient was estimated by 
scaling D for H-0 2 for the larger mass and radius of C0 2 compared to 0 2 (cf. Banks 
and Kockarts, 1973). The H-O diffusion coefficient is taken from the calculation of 
Cooper etal, (1984). The average value for b H at 130 km in the VIRA dayside model 
is 1.8 x 10 19 cm" 1 s~ \ and the diffusion limited flux is 9.0 x 10 7 cm” 2 s ~ If the 
mixing ratio of H 2 is as large as 25 ppm in the mesosphere, and we consider this to be 
unlikely, the diffusion-limited flux could exceed 10 9 cm“ 2 s~ ! . 

4.6. Models of hot and escaping hydrogen 

McElroy and Hunten (1969) showed that the Mariner 5 La data implied a Jeans escape 
flux for H of only about 6 x 10 5 cm -2 s" \ for assumed exospheric temperature of 
about 700 K. The exospheric temperature is now known to be less than 300 K and the 
magnitude of the thermal escape flux has been shown to be negligible (see Section 4.1); 
subsequently the bulk of the escape has been assumed to arise from non-thermal 
mechanisms. Kumar and Hunten (1974) evaluated several sources of hot hydrogen, H*, 
including momentum transfer from the solar wind to the atmospheric gases, resonant 


charge transfer from H + to H 


H* + H-H* + H + , 

(57) 

and the near-resonant process 


+ O^H* + O* . 

(58) 


Since ion temperatures are higher than neutral temperatures at high altitudes, the neutral 
H formed in these reactions will be hotter than thermal. For an H 2 abundance of 
1-2 ppm at the homopause, the major source of energetic hydrogen atoms was 
reaction (54), which is followed by dissociative recombination 

OH * + e-0 + H. (59) 


SPAC 191/67 



68 


J. L. FOX AND S. W. BOUGHER 


Chamberlain (1977) did the first detailed model calculations of the Venus exosphere, 
in which he attempted to reproduce the distribution of hot H inferred from Mariner 5 
Lot measurements, based on charge exchange of H * * with thermal H (reaction (57)) as 
the source of H*. Assuming a spherically-symmetric exosphere, he solved Liouville’s 
equation (the collisionless Boltzmann equation) for the distribution of H and found that 
a higher ionopause than observed was necessary to reproduce the inferred densities. The 
ionopause has, however, been found from Pioneer Venus measurements to be highly 
variable (Taylor et al., 1 979). An escape flux of 2 x 1 0 7 H atoms cm " 2 s “ 1 was derived. 

The Venera 9 and 10 La line shape measurements showed that the temperature 
increase abruptly near an altitude of 3000 km, from which Bertaux et al. (1978) inferred 
that the hot H is not generated near the exobase. They preferred charge exchange from 
solar wind protons as the source of hot H. Kumar et al. (1978), however, argued that 
the experiment could not distinguish an increase in temperature from a net Doppler shift. 

Kumar etal. (1978) evaluated the proposed sources of non-thermal H atoms and 
concluded that reactions (54) and (59) were the most important sources. The production 
of C0 2 H * in reaction (55) followed by dissociative recombination of the ion 

COjH + + e - C0 2 + H (60) 

would also provide a significant, but lower altitude source of H atoms. In order to 
reconcile these large sources with the low thermospheric densities of H inferred from 
Mariner 5 and 10 data. Mayr etal. (1978) proposed that the H atoms created in this 
process are swept to the nightside by thermospheric winds. Kumar et al. found that the 
dayside densities could be reduced by a factor of 1000 by this process. An evaluation 
of the escape rates due to the various processes led to the conclusion that the largest 
source of escaping H atoms was charge exchange with solar wind protons followed by 
reacceleration of the newly formed low-energy protons. An escape flux of 
10 7 H atoms cm " 2 s ” 1 was derived. 

Cravens et al. (1980) presented the first model of the Mariner 5 hot hydrogen distribu- 
tion based on atmospheric and ionospheric structure information from Pioneer Venus. 
Using a two-stream approach to evaluate the fluxes and densities of hot H at the 
exobase, they showed that the unexpectedly hot and dense nightside corona results from 
the high ion temperatures and densities observed by the PV orbiter (e.g., Taylor et al., 
1980; Miller etal., 1980; Knudsen et al., 1979). They found that half the source of hot 
H was due to charge exchange from H * and O + (reactions (57) and (58)), and, for the 
H 2 densities adopted, half was due to the reaction of O + with H 2 (reaction (54)). Only 
preliminary data about the pre-dawn H bulge was available, so the importance of 
reaction (57) was underestimated for the nightside; the H 2 densities in their model were 
taken from Kumar and Hunten (1974), so the importance of reaction (45) was overesti- 
mated. Hodges and Tinsley (1981), however, found that charge exchange alone was 
sufficient to explain the Mariner 5 and 10 La profiles. They constructed a 3-D (non- 
spherically symmetric) Monte-Carlo model of the Venus densities and temperatures, 
and found that the most important source is the charge exchange from O - contributing 
only near the exobase. They criticized the hot H source proposed by Bertaux et al. 
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(1978), charge exchange from solar wind protons, as requiring the deceleration of solar 
wind neutralized protons from 400 km s “ 1 to sub-escape speeds and concluded that the 
magnitude of this source of hot H was negligible. The charge exchange sources ade- 
quately reproduced the day/night ratios of the hot H concentrations and scale heights. 

A factor of 2 uncertainty in the hot hydrogen measurements did not, however, exclude 
the presence of another source, such as reaction (54). The computed global average 
escape rate from charge exchange was about 1.8 x 10 7 cm * 2 s ' '. Kumar et al. (1981) 
predicted an escape rate of 1 x 10 8 cm’ 2 s“ 1 from reaction (54) based on the H 2 
densities they derived from the assumption that the mass-2 ion was H 2 . 

The D/H ratio measured by Donahue et al. (1982), and inferred from the D + /H* 
ratio, an enhancement of a factor of 100 over the terrestrial value, was interpreted as 
implying at least a hundred-fold depletion of the initial H inventory of the planet over 
the age of the solar system. Donahue et al. argued that most non-thermal escape 
mechanisms discriminate against D, but when the abundance of H 2 was enhanced by 
a factor of 100, hydrodynamic escape would take over, so the initial abundance of water 
vapor may be much larger than required by the D/H ratio. Hydrodynamic or transonic 
escape, also called Wow-off, is a thermal loss process in which the atmosphere expands 
rapidly outward. The high frequency of collisions renders the process describable by the 
macroscopic theory of fluid dynamics, rather than by the microscopic kinetic theory. 
The dynamics of atmospheres undergoing hydrodynamic blow-off of hydrogen have 
been discussed by Watson et al. (1981) and the salient features of hydrodynamic escape 
have been reviewed by Hunten (1973b, 1982). Mass fractionation in hydrodynamic 
escape may be large, but the process is not as mass-selective as Jeans escape (Zahnle 
and Kasting, 1986; Hunten et al., 1987). Blow-off becomes an appropriate description 
of the escape process above the altitude where X, the reduced gravitational potential 
energy, is about 2. At this altitude, the rate of expansion of the atmosphere reaches the 
speed of sound ; alternatively, the average thermal velocity of a molecule is approximately 
equal to the escape velocity. Hydrodynamic escape thus becomes the major escape 
mechanism if and when the A = 2 level moves below the exobase (Walker, 1982). 

McElroy etal. (1982) suggested another non-thermal mechanism for the escape of 
hydrogen that discriminates almost completely against deuterium: reaction with hot O 
produced in dissociative recombination of 0 2 at the exobase. This process, which was 
first suggested for He and H escape from Mars by Knudsen (1973), can be represented 
by 

0 2 * + e-*0* + O* < 61 ) 

followed by 

O* + H-»0 + H* . ( 62 > 

Using a Monte-Carlo method and assuming isotropic scattering, McElroy etal. com- 
puted the fraction of collisions (reaction (62)) leading to H atoms with energies greater 
than the escape energy as 15% for an O* with initial velocity of 5.6 km s (2.6 eV) 
colliding with a 300 K H-atom. A planetary average H atom escape rate due to 
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Fig. 26. Hydrogen La spectra of bright disk emission from Venus obtained with the IUE large aperture 
and summing the six best images (solid line). The line intensity corresponds to 21 x 10 3 R', the long-dashed 
line shows the line center. The short-dashed line shows the profile calculated for deuterium Lyman alpha 
emissions of 2.5 x 10 J R (corresponding to a D/H ratio of 1.6 x 10' J ). The numbers refer to spectra 
reference numbers and exposure times. AU denotes 10 5 HJE flux units. Taken from Bertaux and Clarke 

(1989). 


reaction (62) of 8 x 10 s cm" 2 s~ 1 was obtained. Much smaller escape rates due to 
charge transfer from H + (reactions (57) and (58)) were obtained, due partly to the 
assumption of a low (2000 K) ion temperature, compared to PV values, especially on 
the nioghtside, where ion temperatures have been shown to exceed 5000 K (Miller et al., 
1980). 

The escape fraction in reaction (62) computed by McElroy et al. (1982) was shown 
to be overestimated by Cooper et al. (1984), who derived a general formula for the rate 
coefficient for the production of a particle with a specific energy colliding with the atoms 
of a thermal gas. The cross sections for O-H elastic scattering were computed using 
semi-empirical and theoretical interaction potentials, and the effect of the anisotropy of 
the scattering was included. Figure 27 shows their calculated frequency of production 
of a hydrogen atoms with various kinetic energies upon collision with an O-atom 
traveling at 5.6 km s “ 1 . The computed fractions of escaping H-atoms were 5.1, 6.9, and 
8.5% for temperatures of 100, 200, and 300 K, respectively, about half the values 
computed by McElroy et al. (1982). In addition, the O-atom velocity of 5.6 km s “ 1 used 
in the calculations of McElroy et al. was computed assuming that the dissociative 




STRUCTURE, LUMINOSITY, AND DYNAMICS OF THE VENUS THERMOSPHERE 


71 



Fig. 27. The frequency with which a hydrogen atom with energy E is produced in collision with a 2.5 eV 
oxygen atom at temperatures of 100, 200, and 300 K. Taken from Cooper et al, (1984). 


recombination of 0 2 + produces one 0( l D) and one 0( 3 P) atom. As pointed out in 
Section 3.2, it has become clear from calculations and measurements that the yields of 
excited stated in DR of 0 2 (i?) depend on the vibrational level, v (e.g., Guberman, 1987, 
1988; Bates and Zipf, 1980). Further discussion will be deferred to the discussion of 
the hot O corona in Section 4.7. 

Rodriguez et al. (1984) computed the exospheric distribution and escape rates of H 
by an approximate numerical solution to the time-independent Boltzmann equation. 
The model parameters were based on PV data, including the pre-dawn bulge in the H 
densities inferred by Brinton et al. (1980), and the high ion temperatures measured by 
the ORPA (Miller et al , , 1980). Alternative 0 2 * densities were taken from the PV OIMS 
measurements (Taylor et al ., 1980) and the PV orbiter radio occultation (ORO) profiles 
(Kliore et al ., 1979). The H 2 abundance was assumed to be less than 0.5 ppm, in accord 
with the mesospheric photochemical models (Yatteau, 1983; Yung and DeMore, 1982). 
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A more realistic cross section for 0*-H elastic scattering was employed, similar to that 
computed by Cooper et al. (1984). The assumed exobase altitudes, 190 km and 150 km 
on the day* and nightsides, respectively, arc somewhat lower than those used by other 
workers and computed here (see Section 4.1). Rodriguez et al. found that collisions with 
hot O (reaction (62)) and H*-H charge transfer (reaction (57)) were comparable 
sources of hot hydrogen, with the latter contributing more to escape. On the nightside, 
if the OIMS measurements were used for the densities of 0 2 + , reaction (62) was found 
to produce adequate hot hydrogen; if the ORO densities, which are significantly smaller 
than the OIMS values, were used, it was necessary to include reaction (57). An average 
escape flux of about 10 7 cm -2 s~ 1 was computed. 

Hodges and Tinsley (1986) investigated the effect of H-H * charge exchange on the 
velocity distribution and escape of exospheric H in an update of their previous three- 
dimensional Monte-Carlo model of the Venus exosphere (Hodges and Tinsley, 1981). 
Their computed exospheric H-atom densities, presented as a function of solar hour 
angle, show that, due to lateral transport, at increasing planeto-centric distances, the 
pre-dawn bulge becomes less pronounced; outward of 8000 km, the distribution be- 
comes nearly spherically symmetric. Hodges and Tinsley (1986) suggest that a measure- 
ment of the line profiles of La due to resonance scattering of solar photons should carry 
a signature of their source; they have computed the La line profiles for the charge 
exchange source of hot H. The profiles are essentially a superposition of a narrow and 
a broad Gaussian, corresponding to the thermal and hot hydrogen components, with 
asymmetries due to escape and lateral flow. An average planetary escape rate due to 
charge exchange of 2.8 x 10 7 cm “ 2 s " ' was computed. 

Bishop (1989) has discussed modifications to the exospheric distribution and escape 
that result from radiation pressure due to resonance scattering of solar La. The effects 
on the distribution of H are largely confined to distances from the center of Venus greater 
than two Venus radii, where densities are small. A small increase in the hydrogen escape 
flux of 2 x 10 6 cm s ” 1 was computed. 

The most comprehensive global model of non-thermal escape mechanisms and their 
evolution through the history of the planet is that of Kumar et al. (1983). The sources 
of ht H that they considered included charge-exchange of H - with H and impact of 
hot O with H. Their calculations for the production rates due to the reaction of O * with 
H 2 were contingent upon the identification of the mass-2 ion as H 2 * rather than D"\ 
the currently accepted identification. The pre-dawn bulge region of the atmosphere was 
used as an indication of the structure of the thermosphere-ionosphere at earlier times, 
when the hydrogen abundance was larger. At present, the most important source of hot 
H is charge exchange of H*-H (reaction (57)), with a global average escape rate of 
1.2x10 cm 2 s ', and most of the flux occurs over the nightside hemisphere. Impact 
of hot^ O with H is slightly less important; the global average escape rate is 
8 x 10 6 cm 2 s '.At earlier times, the importance of this source declines as the H 
abundances increases and the exobase rises above the 0 2 * peak. Figure 28 shows the 
evolution of the H escape flux due to various mechanisms as a function of the atomic 
hydrogen mixing ratio. The calculations of Kumar et al. show that when the M mixing 
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Fig 28. Hydrogen atom escape flux on Venus as a function of H mixing ratio at the homopause is shown 
for various escape mechanisms. This H mixing ratio is equivalent to twice the HjO vapor mixing ratio at 
the cold trap. The H mixing ratio at the homopause for the present atmosphere is 2.5 x 10 Taken from 

Kumar et at (1983). 


ratio reaches 2 x 10 " 3 , the H density determines the exobase level and the escape rate 
of H levels off at a maximum value of 7.5 x 10 9 cm ' 2 s' When the mixing ratio of 
H exceeds a value near 8 x 10 " 3 , a transition from non-thermal escape to hydro- 
dynamic escape takes place. 

Escape fluxes of H may be derived if the density profile of H below the exobase is 
known because the actual scale height is slightly smaller than the diffusive equilibrium 
value. From La profiles recorded by the PV OUVS, Paxton et al. (1988) deduced a value 
of 7.5 + 1.5 x 10 7 cm " 2 s " 1 for the flux of H through the exobase near the sub-solar 
point. The model calculations discussed above show that the escape flux is a quarter 
to a third of this value. The remainder of this flux could be accounted for by the sub-solar 
to anti-solar circulation, and it would be balanced by a downward flux on the nightside. 

Detailed models for the evolution of the Venus atmosphere have been presented by 
Watson etal. (1981), Kasting and Pollack (1983), and Watson etal. (1984). These 
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models assume that the ultimate source of the escaping hydrogen is water vapor, the 
predicted initial abundance of which is. therefore, much greater than at present. Toward 
earlier times, the larger amount of water vapor in the atmosphere increases the infrared 
opacity and causes the tropospheric temperature to nse. A critical value of the insolation 
exists for which liquid water cannot exist on the surface and the evaporation of the water 
produces a run-away greenhouse (IngersoU, 1969). ITte insolation at Venus is proposed 
to be above the critical value. Kasting er al. (1984) showed that an increase in the water 
vapor mixing ratio also causes a decrease in the tropospheric lapse rate, causing the 
cold-trap to move to higher altitudes and lower pressures. The resulting increase in t e 
mixing ratio of water vapor at high altitudes, the parameter that determines the atomic 
hydrogen escape rate, is dramatic. The predicted initial water abundance vanes greaUyn 
from one model to another, depending on such factors as the mixing ratio of onset of 
hydrodynamic escape, the H/D fractionation factor during such escape, the magnitude 
of the hydrodynamic escape flux (which depends on assumptions made about the 
magnitude of the solar flux in the past), and the outgassing rate of volatiles from the 
interior (e.g., Watson et al., 1981, 1984; Donahue era/., 1982; Kasting and Pollack, 
1983; Kumar era/., 1983). 

An alternative scenario has been proposed by Gnnspoon and Lewis (1988). 
current abundance of water vapor in the atmosphere could be in steady state as a result 
of a balance between loss by escape of H and a cometary source. As in the case of 
continuous slow outgassing, an enhanced D/H ratio could result from this exogenous 
source, but the value depends on the assumed rate of impact of large comets, which is 
somewhat stochastic and difficut to estimate. Nevertheless, a large initial abundance of 
water may not be necessary to explain the enhanced D/H ratio. 

A detailed discussion of atmospheric evolution is beyond the scope of this review. The 
interested reader is referred to a review of the evolution of the atmosphere of Venus 
presented by Donahue and Pollack (1983) and to a book on the origin and evolution 
of planetary atmospheres (Atreya el a/., 1989), especially to the articles by Hunten et al. 
(1989) and Kasting and Toon (1989) contained therein. 


4.7. Hot oxygen corona 

Dissociative recombination of molecular ions produces energetic fragments that are 
important to the thermal structure, coronas and escape from planetary atmospheres. A 
short review of the importance of dissociative recombination in aeronomy has been 
presented by Fox (1989a). Dissociative recombination of 0 2 + can proceed according to 
a number of energetically allowed channels, with exothermicities as: 


0 ; + e - O ( 3 P) + O (V) + 6.98 eV , 
-*0(‘D) + O {'P) + 5.02 eV, 
-O('S) + O CP) + 2.79 eV, 
-.0 ('D) + 0('Z>) + 3.05 eV , 
-.0 ('£>) + O ('S) + 0.83 eV. 


(63a) 

(63b) 

(63c) 

(63d) 

(63e) 
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Since the escape energy per unit mass for Venus is 0.54 eV amu ' *, none of the channels 
produce oxygen atoms with enough energy to escape the gravitational field of the planet 
Oxygen atoms produced above the exobase will, however, travel along ballistic trajec- 
tories in the exosphere, producing a corona of hot oxygen atoms, similar to the hot 
hydrogen corona. A hot oxygen corona has been observed to surround the Earth (Yee 
et al 1980; Yee and Hays, 1980) and the density distribution has been modeled by Yee 
etal. (1980). The existence of hot oxygen coronas around Venus and Mars were 
predicted by Wallis (1978). The corona has a significant effect on the interaction of the 
planet with the solar wind (e.g., Biermann etal., 1967; Wallis, 1972, 1978, 1982) A 

n ooo\ of hydrogen ox ygen coronas of Venus and Mars has been presented by Nagy 
(lyoy). 


Bertaux et al (l 98 1 ) reported the existence of ‘a strange feature’ in the intensity profile 

^ t a \ u thC SpeCtr0ph0t0met er on Venera 1 1. Over a range of about 
beyond the bnght limb, a nearly constant intensity of about 500 R was 
measured Because this feature was not observed by Venera 12, it was tentatively 
interpreted as being due to the presence of ‘a sporadic hot component of exospheric 
oxygen , and an average density of 1 .6 x 1 0 3 cm - 3 was inferred 
Lhnb scans of the PV OUVS also showed the existence of an exospheric signal at 
304 A (Nagyer al., 1981). The oxygen atom density profiles derived from the intensities 
exhibit two scale heights, corresponding to an inner thermal component and an outer 
hot component. The hot component, shown in Figure 29, dominates beyond an altitude 
of about 350 km, where the O density is about 5 x 10 4 cm - 3 . The scale height of the 


Hot Oxygen Densities at Venus 



■O Paxton B R 
R & N. B R 
-o- PVObs. 


Fig. 29. Measured and calculated hot oxygen densities at Venus for different branching ratios in thr 

1 °' (react r 63) - ' Pa * t0n B R ' deno,es branchin * from p«wn (1983) 
RAN B.R. denotes branch.ng ratios from Rohrbaugh and Nisbet (.973), The experimental values are 

enved from PVOUVS 1304 A data. From Nagy and Cravens (1989). 
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hot component is approximately 400 km, indicating a "temperature’ of 6300 K. Nagy 
et al. (1981) constructed the first models of the hot 0 corona on both the dayside and 
the nightside using PV data. In addition to dissociative recombination of 0 2 * , they 
considered the charge exchange processes 


O* + O-O* + O* 

(64) 

cr + H-O* + H + , 

(65) 


as sources for the energetic oxygen, but found these reactions to be relatively unimpor- 
tant. A modified two-stream approach, similar to that used in the hot hydrogen model 
(Cravens el al ., 1980) was used to compute the flux of O at the exobase; the exospheric 
distribution was then determined by solving Liouville’s equation. Nagy et al. obtained 
hot O densities a factor of 4-5 larger than the densities derived from the PV measure- 
ments. The discrepancy was later attributed to the use of preliminary ion densities that 
were too high, and was resolved with the use of appropriate 0 2 + densities (Nagy and 
Cravens, 1988). The corrected calculated hot O density profile is shown also in 
Figure 29. In addition, the earlier paper (Nagy et al ., 1981) assumed very low exobase 
altitudes of 172 and 143 km on the day- and nightsides, respectively, compared to the 
values computed here, which would also tend to exaggerate the magnitude of the O 
source. 

Reasonably successful models of the hot O corona were also presented by McElroy 
et al. (1982b) and by Paxton (1983). McElroy etal. (1982b) presented the first calcu- 
lations of the hot O column density as a function of solar zenith angle from 0 to 90°. 
The results, presented here in Figure 30, show that the column density varies from about 
3 x 10 ,2 cm -2 at the subsolar point to 1 x 10 ,2 cm -2 at the terminator. 

Models of the hot oxygen corona depend on assumptions made about the relative 
importance of the channels (63a-e) in dissociative recombination of 0 2 + . No definitive 
measurements of the yields of 0( 3 R), 0(‘Z)), and O('S) in dissociative recombination 
of the lower vibrational levels of 0 2 + (v) are available. It has become clear that the 
channel by which the dissociative recombination proceeds depends greatly on the 
vibrational state of the ion (Guberman, 1983, 1987, 1988; Bates and Zipf, 1980; 
Queffelec et al., 1989). Zipf ( 1970) determined yields of 0. 1 and 0.9 for 0(‘S) and 0('D), 
respectively, but these values were later withdrawn because the vibrational distribution 
of 0 2 in the experiment was unknown. Guberman (1987, 1988) carried out ab initio 
calculations of the rate coefficients for production of O('S) and O('Z)) in dissociative 
recombination of 0 2 in various vibrational levels. Some of the rate coefficients he 
obtained for a temperature of 300 K are shown in Table X. Specifically, he showed that 
the rate coefficient for production of O('S) increases nearly two orders of magnitude 
from i; = 0 to v = 2! Unfortunately, Guberman did not calculate the rate coefficient for 
the channel that leads to two ground state atoms (reaction (63a)), so the yields of excited 
states cannot be derived from his calculations alone. Recently, Queffelec et al. (1989) 
have studied the yields of ground and excited O atoms for (t) in vibrational levels 
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Fig. 30. Column densities of neutral oxygen above Venus plasmapause as a “f n 
« * in 30 atoms. Taken from McElroy et ai (1982b). 


around r = 9 using a plasma flow tube experiment. The reported yields were ^0.44 TO 
and 0 56 for O('S), 0('D) and 0( 3 />), respectively. An average yield of about 0.1 was 
obtained for the yield of 0(‘S) in dissociative recombination of °^ ^^ nkn °^ 
distribution of vibrational levels 0, 1, and 2, in disagreement with Guberman s calcu 


TABLE X 

Dissociative recombination coefficient (cm J s’ ') « for produc- 
tion of O('D) and 0(‘S) from 0,*(») in different vibrauonal 
levels at 300 K From the ab inilio calculations of Guberman 
(1987, 1988). 


K'S ) 

3.0 x 10- 10 22 x 10 ~’ 

7.3 x 10-’ 18 x 10- 

2.4x10-' 1.2x10- 
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lations. The calculations of Guberman and the experimental study of QuefTelec et al. 
(1989) do show, however, that the vibrational distribution of needs to be known 
in order to compute the production rates of excited states or the energy distribution of 
the O atoms. 

Fox (1985) has modeled the vibrational distribution of 0 2 * in the thermosphere of 
Venus; at the exobase 55% of the 0 2 * ions are predicted to be vibrationally excited. 
Using the rate coefficients of Guberman and measured values for the total dissociative 
recombination rate, Fox (1990b) has computed the yields of the channels above and the 
fractions of O atoms with various energies produced at the exobase on Venus. The major 
uncertainty in this calculation is the yield of channel (63a), which is obtained as a 
difference between the measured total rate and the rates of the channels computed by 
Guberman. The total rate for dissociative recombination almost certainly depends on 
the vibrational level of 0 2 * ; that effect was only accounted for only by the choice of the 
rate coefficient of Mul and McGowan (1978) over that of Alge et al. (1983), since the 
merged beam experiments are expected to contain a larger fraction of vibrationally 
excited ions, and are therefore more appropriate to exobase altitudes, than Bowing 
afterglow values, which pertain to thermalized conditions for temperatures of 200 to 
600 K. The branching ratios for the channels (63a-e) derived from the model of Fox 
(1990b) are shown in Table XI, where they are compared to those assumed in other 


TABLE XI 


Yields of various channels in dissociative recombination of Oj computed or assumed in 
calculations of the hot oxygen coronas 


Channel 

Energy (eV) 

Yield 




a 

b 

c 

d 

i P + i P 

6.98 

0.325 

0.22 

0.0 

0.334 

} P * 'D 

5.01 

0.30 

0.55 

1.0 

0.393 

’P+'S 

2.79 

0.05 

0.0 

0.0 

0.0 

'D + 'D 

3.04 

0.275 

0.13 

0.0 

0.233 

'D + '5 

0.82 

0.05 

0.10 

0.0 

0.038 


(a) Derived by Rohrbaugh and Nisbet (1973) from data of Zipf (1970). Used by Nagy et al. 
(1981) and Nagy and Cravens (1988). 

(b) Determined by Paxton (1983) from a best fit to the PV OUVS hot oxygen data. 

(c) Assumed by McElroy et al. (1982a). 

(d) Derived by Fox (1990b) from calculations of Guberman (1987, 1988) and the measure- 
ments of Mul and McGowan (1979). 


models. The branching ratios agree surprisingly well with those derived by Rohrbaugh 
and Nisbet (1973) from the measurements of Zipf (1970), that were used in the models 
of Nagy etal. (1981) and Nagy and Cravens (1988). The agreement is poorer with the 
values obtained by Paxton (1983), who varied the branching ratios to obtain agreement 
with the PV hot oxygen data. McElroy et al. (1982) assumed that all the dissociative 
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recombinations proceed by channel (63b). While this choice of branc hing ratios dis- 
agrees with those of other workers, the average energy of the 0 atoms is approximately 
correct, so fairly good agreement with experiment was obtained. 

4.8. Escape of O 

The large escape flux of H, the source of which is a putative large initial inventory of 
water or a continuous endogenous or exogenous source (see Section 4.6), raises the 
issue of the fate of the O atoms remaining behind, given the dearth of 0 2 in the Venus 
atmosphere. The O atoms produced in dissociative recombination of Oj’ are not 
energetic enough to escape on Venus, as they do above the Martian exobase. Most 
workers have assumed that the O atoms ‘disappear’ by oxidizing either crustal iron, Fe 
or Fe 2 to Fe 3 + , or atomic carbon or CO rom the planet’s initial endowment of volatiles 
to produce some of the C0 2 that is in the atmosphere today (McElroy etal. , 1982; 
Donahue et al., 1982; Lewis and Prinn, 1984). Prinn (1985) has described some of the 
problems with these scenarios, which include the requirement that a fair fraction ( 1-3 % ) 
of the interior of Venus be exposed to the atmosphere over the first 0.3 billion years or 
so of its existence. Models of the evolution of the atmosphere indicate that the 
weathering process may be accelerated by the presence of liquid water on the surface 
(Kastinger a/., 1984) or by the presence of a molten surface due to a runaway greenhouse 
(Watson etal., 1984). 

McElroy et al. (1982b) and Wallis (1982) proposed that O escapes from the Venus 
atmosphere by photoionization and electron impact ionization of the rmal and hot O 
above the plasmapause, and subsequent pick-up of O + by the solar wind. McElroy et al. 
estimated the total quantity of O atoms as 8 x 10 30 and the loss rate about 
6 x 10 6 cm ~ 2 s ” '. They asserted that this escape flux was about half the escape flux 
of H, and that the escape rates of H and O are regulated by the oxidation state of the 
lower atmosphere. This analysis was, however, based on the H escape flux of 
1.2 x 10 7 cm " 2 s _ 1 computed by McElroy et al. (1982a), which is smaller by a factor 
of two or more than subsequent models indicate (see Section 4.6). 

Recently, Luhmann and Kozyra (1990) have shown that approximately 90% of the 
O + ions created by photoionization above the plasmapause and picked up by the solar 
wind reimpact the planet, rather than being swept away. Significant and perhaps even 
enhanced escape still occurs due to sputtering of O atoms by the precipitating O. They 
computed a global loss rate of 2 x 1 0 25 O atoms s" 1 or an escape flux of 
6 x 10 6 cm “ 2 s “ ', which is comparable to the escape flux computed by McElroy et al 
(1982a). 

4.10. C AND N HOT ATOM CORONAS 

Paxton (1983) modeled the production rate of hot carbon due to knock-on from the hot 
O produced in dissociative recombination of 0 2 + near the exobase (reactions 63(a-e)): 

O* + C -* C* + O . (66) 

A smaller contribution, about 10% of that from reaction (66) was determined to arise 
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from dissociative recombination of CO * : 

CO* +e-C* + OV (67) 

The altitude distribution of C* was derived from PV OUVS limb scans of the high 
altitude optically thin intensity at 1657 A, assumed to arise from resonance scattering 
of solar radiation by atomic carbon. A C/O ratio of about 1% at the exobase was 
determined and the scale height of the altitude profile indicates a temperature that is 
approximately twice thermal. 

Limb profiles of the signal at 1493 A, a feature produced by resonance scattering of 
solar radiation by N( 2 Z)), were modeled by Paxton to determine the density distribution 
of hot N( D) (Keating et al. , 1985). The major source was determined to be dissociative 
recombination of N 2 , with lesser contributions from dissociative recombination of 
NO* and knock-on from hot O. 

4.10. He escape 

Thermal escape of He from the atmosphere of Venus is negligibly small, but escape is 
possible by photoionization and electron impact ionization of He atoms above the 
ionopause, followed by pick-up of the resulting He * ions by the solar wind. Kumar and 
Broadfoot (1975) used Mariner 10 data to estimate an escape rate of 
2 x 10 5 He atoms cm ~ 2 s " 1 due to this source. This estimate was based on limited 
observations of the ionopause altitude, which has been found to be quite variable by PV 
observations (Taylor et al ., 1979). Prather and McElroy (1983) have estimated the total 
escape rate as 1 x 10 6 cm~ 2 s“ *, with an accuracy of about + 30%. It is interesting 
that they find that about 20% of the He above the ionopause is produced by knock-on 
from hot O, the process suggested by Knudsen (1973) for escape of He from Mars. 
Chassefiere et al. (1986) found, from Venera 1 1 and 12 observations of the He emission 
at 584 A, that the abundance of He is a factor of two smaller than the value derived from 
PV data and, therefore, the estimated escape rate, 5 x 10 5 cm “ 2 s - \ is also therefore 
smaller by a factor of 2. 

In the terrestrial atmosphere, escape of 4 He is balanced by production in radioactive 
decay of uranium and thorium in the crust. The terrestrial escape flux is about 
2 x 10 6 cm 2 s 1 (MacDonald, 1963). Prather and McElroy conclude that the com- 
parable magnitude of the sources of 4 He on Earth and Venus suggests comparable 
abundances of uranium and thorium and, along with independent data about other 
volatiles, indicates that Venus and Earth had similar origins. 


5. Thermal Structure 

5.1. Global mean heat balance and thermal structure 

The calculation of a reasonable dayside Venus mesosphere/thermosphere heat budget 
yielding observed temperatures has been an ongoing problem for nearly two decades. 
The difficulty is illustrated in Figure 31 (Keating and Bougher, 1987), which compares 
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Fin 3 1 Temperatures of the neutral atmospheres of Earth, Venus, and Mars. Adapted from Schubert et al. 

(1983). 


Venus, Earth, and Mars dayside thermospheric temperature profiles; Venus, although 
closest to the Sun, has the coldest temperatures above 100 km. This implies that, unlike 
the terrestrial thermosphere, a simple balance of EUV heat and molecular conduction 
may not be sufficient to maintain the observed relatively cold ( <> 300 K) temperatures. 
Several modeling efforts have attempted to determine what additional cooling mecha- 
nisms might be responsible. Uncertainties in the role of eddy conduction, the 0-C0 2 
collisional excitation rate, the EUV heating efficiencies, and the EUV solar fluxes 
themselves make the job very difficult. Recently additional constraints have been 
provided by the Pioneer Venus neutral thermospheric data, new laboratory measure- 
ments, and terrestrial modeling and observational efforts. Model simulations can be 
made in which EUV heating efficiency, C0 2 15-pm cooling, and eddy conduction are 
considered to be free model parameters. However, the range of realistic heating efficien- 
cies reduces the possibilities considerably (see below). 

5.2. EUV HEATING EFFICIENCIES 

Heating due to the absorption of solar ultraviolet radiation occurs by photoionization 
and photodissociation. In photodissociation, the energy in excess of the dissociation 
threshold may be converted to internal or kinetic energy of the neutral fragments. Kinetic 
energy can be converted directly to heat through elastic collisions, or to internal energy 
through inelastic processes, such as collisional excitation of rotational or vibrational 
modes of molecules. Excited vibrational or electronic states may radiate, which results 
in cooling if the radiation escapes from the atmosphere, or they may be quenched, 
releasing their internal energies as heat; rotational energy is usually rapidly thermalized. 
In photoionization, the excess energy is carried away by the ejected electron, which can 
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produce further ionization, dissociation or excitation. Some of the chemical energy of 
the ions is released as heat in exothermic charge transfer reactions and ultimately in 
neutralization through dissociative recombination. 

It is traditional and convenient to define the heating efficiency, e, as the fraction of 
solar energy absorbed at an altitude that is converted locally to heat. Heating efficiencies 
for Venus were first computed by Henry and McElroy (1968). Their computed mean 
heating efficiency, about 60%, was based on the (incorrect) assumption that C0 2 + was 
the major ion in the ionosphere. The predicted global mean exospheric temperature was 
about 700 K (McElroy, 1969). Stewart and Hogan (1969) estimated heating efficiencies 
about half that value from studies of the Martian thermosphere in which the heating 
efficiency was varied to reproduce the available data. They reasoned that the long Venus 
day would allow radiative equilibrium to be established on the dayside, so the computed 
(dayside) exospheric tempedrature was the same as McElroy’s mean value (Hogan and 
Stewart, 1969). This contention was criticized by McElroy (1970), who pointed out that 

the Venus atmosphere could not sustain the large day-night density and temperature 
gradients. 

In the years that followed, the exospheric temperature was determined to be much 
lower (^ 300 K) than these early models indicated. Dickinson (1976) and Dickinson 
and Ridley (1977) computed 1-D and 2-D models of the thermal structure in the Venus 
thermosphere; they found that the low exospheric temperatures required very small 
values of the heating efficiency. They proposed that much of the energy involved in solar 
energy absorption would be taken up as vibrational and rotational excitation of C0 2 
which has a larger number of degrees of freedom than the diatomic molecules that 
dominate the terrestrial atmosphere. Fox and Dalgamo (1981) computed heating effi- 
ciencies in the thermosphere for three models with different assumptions about the 
fraction of the excess energy that is converted to vibrational excitration in photodisso- 
ciation and chemical reactions. In the standard model, for the assumption that 50% of 
the excess energy was converted to vibrational excitation, the heating efficiencies varied 
from about 1 8 to 22% over the altitude range from 115 to 235 km. Studies of the thermal 
structure of the Venus thermosphere in light of the knowledge gained from Pioneer 
Venus indicated that even heating efficiencies of that magnitude could not be accommo- 
dated and suggested that the heating efficiencies were on the order of 10% (Dickinson 
and Bougher, 1986; Hollenbach et ai., 1985). 

Fox (1988) considered the likelihood that 10% heating efficiencies could be justified 
from a molecular viewpoint. Both Hollenbach et al. and Dickinson and Bougher 
partially justified the low heating efficiencies by the assumption that 90% or more of 
the exothermicity of the quenching reaction 

O('Z)) + C0 2 -* 0( 3 P) + C0 2 

is converted to vibrational excitation of C0 2 . Fox adduced evidence that reaction (68) 
proceeds via a collision complex, which should produce statistical energy partitioning 
if the lifetime of the complex is long enough. The vibrational distribution that would be 
produced for statistical energy partitioning was computed by applying information 
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theory (e.g., Levine and Bernstein, 1974; Bernstein and Levine, 1975). It was found that 
only 55% of the excess energy appears as vibrational energy. An upper limit was 
determined from a survey of the literature, which indicated that a value of 75 % would 
be high, but not unreasonable. After a consideration of the probable values and upper 
limits for energy deposited as vibrational excitation in photodissociation (about 25%, 
with an upper limit of about 35%) and chemical reactions (50-60%), and energy 
released in dissociative recombination of 0 2 * , including a model of its vibrational 
distribution, altitude profiles of the heating rates, shown here in Figure 32, were 



Fig. 32. Altitude profiles of heating rates due to the major source of neutral heating m the Venusian 
thermosphere from 1 15 to 200 km from the standard model, in which most probable values for the fraction 
of energy deposited as vibrational excitation in various elementary processes were employed (see text). 
^ Taken from Fox (1988), 


computed. The major source of heating above about 130 km is dissociative recombi- 
nation of O; (reaction (63)). Below that altitude, photodissociation and quenching of 
metastable electronic states are about equally important. Heating efficiencies for stan- 
dard and lower limit models were computed and the results are shown in Figure 33. The 
standard model employs the probable values for the fraction of energy deposited as 
vibrational excitation in the processes discussed above and the lower limit model uses 
the upper limits for fraction of energy deposited as vibration excitation. Below 125 km, 
the heating efficiencies for the lower limit and standard models are about 16 and 22%, 
respectively, and above 130 km, they are 22 and 25%, respectively. The heating efficien- 
cies are slightly higher and the difference between the standard and lower limit models 
is smaller than the calculations of Fox and Dalgamo (1981) indicate. Fox (1988) found 
that heating efficiencies on the order of 10% or less were difficult to justify with 
reasonable assumptions about molecular processes in the Venus thermosphere. 

The FI ion density peak at 140 km (t = 1 for EUV), composed mostly of O; .should 
correspond closely to the maximum in per volume (eV cm “ 3 s 1 ) heating if 0 2 disso- 
ciative recombination and related ion-neutral chemical reactions dominate other heating 
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6 

Fig. 33. Computed heating efficiencies for the altitude range 115-200 km. The curve labeled A is the 
standard model and B is the lower limit model. The curve labeled C is taken from Hollenbach et ai (1985). 

Taken from Fox (1988). 


mechanisms (see Figure 32). But the magnitude of the heating at t = 1 is a function of 
the specific EU V fluxes and the heating efficiency used. The height of this level depends 
on the column of absorber and thus on the neutral atmosphere chosen. Some of the 
discrepancies that arise in the temperature calculations and the heating efficiencies may 
be due to the use of different inputs. Aji effort by the authors is underway to compute 
heating efficiencies and model thermal structure using common background atmosphe- 
ric models and solar fluxes. 
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5.3. The dayside heat budget 

One-dimensional modeling efforts for global mean condiuons data b«k u |toe of 
McElroy (1968, 1969) and Dickinson (1972). for a nearly pnre CO, atmosphere, be 
me advent of Venus general circulation models. In retrospect, n appears that the 
of calculating temperatures based upon a balance of rad«»ve 
si and thermal condition is adequate to address the Venus days.de hea. bdmce 
nroblem Large-scale upwelling winds appear to provide little adiabauc cooling affecting 
die dayside thennospheric heat budge, (see Section 6.3). Therefore. 
are no. required to explain die observed dayside temperatures. As a result, 1-D tnode 
simulations for vanous combinations of input parameters prosnde a u«M andr«lau y 
efficient method for examining heating and coolrng processes on the Venus 

Dickinson (1972) chose to address both the Venus mesosphere and thermosphere (65 
to 200 km), where, as for the Earth, a non-LTE (NLTE) treatment is 
d articular his 1-D radiative transfer calculation utilized a two-level NLTE fonnulation 
for the C0 2 vibrationally excited states. The primary radiative transfer was assumedto 
onW in the 15-pm and 4.3-pm band systems, for which relative populations of the 
v bending mode and v 3 asymmetric stretch vibrational levels were painstakmgly calcu- 
lated Near- 1 R heating was realized through the collisional relaxation of the 15-pm 
fundamental band. Several isotopic and hot bands yielding 15-pm emission were also 
considered Dickinson (1972) calculated a global mean temperature profile assummg an 
rl mdcpcudcn, EUV efficiency of W/. , us suggest* b, curly mgumeuts 
(1972) This large efficiency was determined with the assumption ffiat all the 0( D) 
energy is transformed to heat. Below 1 15 km, 1R processes ^ themselves are <^ance ; 
with heating and cooling rates decreasing from 300 K day at 1 15 km to y 

r65km UV absorbed above 160km (*2500 K day- «) was almost entirely 

offset by downward molecular conduction. This conducted heat was deposited between 
120 and 140 km, and was balanced by 15-pm cooling. Temperatures ranged from 
minimum of ,80 K a. .22 km (the metopautc) to 4,5 K at the top. Dickmson found 
that the mesopause occurs at essentially the lowest pressure at which C0 2 pm 

C °^?rlutTn«rtdTty 1 ffieEUV heating efficiency, as weU as other par^eterv 
rations, lead Dickinson (1976) to conduct several sensitivity studies ^ 

offirienev eddv mixing and C0 2 collisional excitauon by atomic-0 (enhancing 5 pm 
2S?’ it^Twete varied. Studies of Dickinson (1976) and Fox and Ddgamo 
(1981) both indicated that 30% is probably an upper limit fc .eating * > c s ° n 
1976) suggested that the very rapid collisional quenching of 0( D) by GO, s 
with the concomitant excitation of upper C0 2 vibrational-rotational states, voca- 
tionally excited states at the low pressures of the Venus thermosphere lose their energy 
to IR airglow, rather than by being quenched. Assuming ffiat ^cant fraction o 
the kinetic energy of hot atoms, ions, and molecu es is so ans * 3 t jqo - 

increased coohng (or a reduced heating efficiency) is obtamed. Effic.enc.es of S 30 
were examined, yielding global mean exospheric temperatures of 250 to 475 K 
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(compared with later Pioneer Venus values of —215 K). Thermospheric temperatures 
can be further reduced in the presence of O-atoms, for which the coilisional excitation 
of the C0 2 v 2 bending mode: 

C0 2 (r,, v 2 , V}C + O* -► C0 2 (v u v 2 + 1, v 3 ) + O (69) 

is more efficient. Dickinson argued tha tO excites the C0 2 15-pm level possibly ten times 
as rapidly per collision as do self-collisions (corresponding to a de-excitation rate of 
*co 2 -o =* 5 x 10" 14 cm 3 s' ') (see Section 5.4). This assumption resulted in a further 
reduction of model exospheric temperatures by 25-50 K. Dickinson considered the 
introduction of eddy cooling to be the least satisfactory method for modifying 
themrospheric temperatures (see Section 5.5). Calculations showed that eddy cooling 
using eddy diffusion coefficients ^ 10 6 cm 2 s “ 1 does not significantly alter thermosphe- 
ric temperatures. However, values as large as 10 7 cm 2 s' 1 have a substantial impact 
on the dayside heat budget. Temperatures near the mesopause are particularly sensitive 
to the amount of eddy heat conduction prescribed. Dickinson made the common 
assumption that the same eddy coefficient applies to eddy diffusion and eddy con- 
duction ; that is, mass and heat were assumed to be transported with the same efficiency. 

Post-Pioneer Venus heat budget studies were first carried out by Hollenbach et ai 
(1985), who stressed the role of eddy cooling in reproducing the observed thermospheric 
temperature structure using a one-dimensional model. They suggested that the changing 
Venus circulation, as simulated by an effective eddy coefficient, adjusts to maintain a 
nearly constant thermospheric temperature, yielding exospheric values over the solar 
cycle that depart little from 300 K. Heating efficiencies of 12-15% were internally 
calculated and used for all simulations. The large O + C0 2 de-excitation rate (the 
inverse of reaction (69)) suggested by Sharma and Nadile (1980, 98 1) and Gordiets et ai 
(1982) (^ 1 x 10“ 12 cm 3 s' *) was found to be incompatible with their model calcu- 
lations. Later circulation model studies of Bougher et ai (1986) showed that Venus 
dynamics has little impact on the dayside thermospheric heat budget; i.e., the circulation 
is not an effective thermostat regulating dayside temperatures (see Section 6.3). 

Venus heat budget studies were also carried out by Gougher (1985) and Dickinson 
and Bougher (1986) using an updated version of the previous NTLE mesosphere/ther- 
mosphere code (Dickinson, 1972, 1976). They stressed the role played by C0 2 15-pm 
cooling in maintaining observed temperatures. Several modifications were made that 
affected the calculated exospheric temperatures significantly. Solar EUV and FUV 
fluxes were taken from the Torr et ai (1979) and Torr et al. (1980) tabulations for 
day 78348, corresponding to near solar maximum conditions (F l 0 7 = 200) of the early 
in situ Pioneer Venus observations (1978-1981). These changes essentially doubled the 
UV heating rates from those previously calculated, making the problem of obtaining low 
exospheric temperatures that much more difficult. In addition, Pioneer Venus measure- 
ments indicated that dayside atomic-0 densities are about 3 times larger than previously 
thought (von Zahn et n/., 1983). This meant that the 0/C0 2 ratio at the dayside 
ionospheric peak (^ 140 km) was 17-20% rather than 5.6%, as calculated earlier by 
Dickinson and Ridley (1977). This improvement increased the efficiency of excitation 
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of the C0 2 vibrational v 2 - 1 state by O collisions, thereby enhancing 15-|ira emission 
and cooling where NLTE conditions prevail. Finally, the O + C0 2 energy transfer rate 
itself was increased from that previously used, based in part on recent studies of the 
Earth’s lower thermosphere as reviewed by Dickinson (1984), and observational esti- 
mates for this rate coefficient provided by terrestrial rocket measurements (Sharma and 
Nadile, 1981; Stair et al., 1985). A value of 5-8 x 10' 13 cm' 3 $ " 1 at 300 K was 
adopted for this rate (Dickinson and Bougher, 1986; Bougher et ai , 1986), in contrast 
to2x 10' 13 used in Dickinson (1984) for Earth, and about 5 x 10“ 14 used in previous 
Venus studies (e.g., Dickinson, 1976). This amounts to a ^ 150-fold increase over C0 2 
self-collisions. 

It was found that the observed Pioneer Venus global mean temperatures above 
140 km could be simulated using strong 15-pm cooling alone for balancing ~10%- 
efficient EUV heating. The maximum EUV heating rate (~ 1100 K earth day" ’) was 
offset by a combination of molecular conduction and 15-pm cooling. The region above 
160 km, where molecular conduction removes most of the heat, is nearly isothermal. 
Thus, the exospheric temperature of Venus was shown to be rather insensitive to the 
precise details of molecular thermal conduction, if C0 2 cooling is indeed so strong. The 
near IR heating rate was nearly doubled from that of Dickinson (1972), with its peak 
altitude rising from 1 15 to 130 km. The global mean model was also decomposed into 
a dayside and nightside model, with solar IR heating partitioned between the two in 
order to simulate hydrodynamical effects. An improved 15-p.m cooling rate was obtained 
by averaging these day and night profiles, yielding a reduced global mean exospheric 
temperature of 220 K. This is consistent with a dayside mean value of 300 K, and a 
corresponding nightside value of 130 K. Very efficient eddy cooling was shown to be 
required to yield observed temperatures in the face of larger EUV heating efficiencies 
(10-18%). However, such a possibility was generally ruled out since Venus circulation 
model studies (Bougher et al., 1986) require that eddy diffusion must be less than half 
that proposed by previous one-dimensional eddy diffusion models (e.g., von Zahn et al . , 
1980). Hydrodynamic model large-scale winds appear to be very important in 
maintaining the observed dayside density profiles (see Section 6.6). Therefore, the 
requirement that the eddy coefficient for diffusion and conduction be the same 
constrains the cooling to be rather small (K £2 x 10 7 cm 2 s " 1 ). 

The low thermospheric heating efficiencies of 10-15% used by Hollenbach et al. 
(1985) as well as by Dickinson and Bougher (1986) have recently been shown to be 
improbable (Fox, 1988) (see Section 5.2). The narrow range of permissible EUV heating 
efficiencies (22-25% above 135 km) requires a re-examination of th Venus dayside heat 
budget. The NLTE model of Dickinson and Bougher (1986) has recently been updated 
to determine what reasonable combination of cooling parameters and EUV fluxes might 
enable observed temperatures to be simulated using EUV heating efficiencies consistent 
with the calculations of Fox (1988) (cf. Bougher etal., 1988a). The model EUV fluxes 
used have been revised according to Torr and Torr (1985), where solar maximum values 
are reduced as a result of a rocket calibration. The effect on the heating rates permits 
an increased EUV efficiency of 12% to be "lopted. If dayside solar maximum tempera- 
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turcs as large as 310 K are permitted, in accord with the Hedin ex al. (1983) empirical 
model, then an additional 1% increase of the EUV efficiency to 13% is possible. Finally, 
recent terrestrial low thermosphere heat balance studies of Roble etal. (1987, 1989) 
suggest that an O + C0 2 de-excitation rate similar to that initially proposed by Sharma 
and Nadile (1980) ( — 1 x 10" l2 ) is not unreasonable. Its incorporation into the Venus 
NLTE code enhances the C0 2 cooling further, so that an EUV efficiency of 14% can 
be used. Figure 34 illustrates the cooling requirements for observed dayside tempera- 


VENUS COOLING REQUIREMENTS FOR 
OBSERVED DAYSIDE TEMPERATURES 
(SMAX) 



Fig. 34. Venus cooling requirements for observed dayside temperatures. Solar maximum conditions 
(FI 0.7 = 200). Based on the NCAR NLTE energy balance code. 


tures; values along the ordinate correspond to the case of no eddy cooling. A further 
increase of EUV efficiencies seems to require cooling by eddy conduction. This, in fact, 
appears valid for moderate eddy coefficients with K ^ 2 x 10 7 cm 2 s“ \ which is con- 
sistent with an upper limit for eddy diffusion incorporated in Venus odd-nitrogen studies 
(Bougher et a/., 1989) (see Section 3.10). The use of EUV efficiencies larger than 15%, 
as Fox (1988) has suggested they are, seems impossible at the present time, unless the 
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requirement that eddy diffusion and conduction coefficients be equal is relaxed. The low 
efficiencies required may, however, also be a due to the uncertainty of the solar EUV 
fluxes and the 0-C0 2 energy exchange rate being used in model simulations. 

5.4. IR HEATING AND 15-pm COOLING 

Both IR heating and C0 2 15-pm cooling appear to be quite important in the Venus upper 
atmosphere. Dickinson and Bougher (1986) showed that IR heating and cooling offset 
one another below 130 km. Between 130 and 160 km, where the temperature increase 
is greatest, EUV heat is balanced by C0 2 15-pm cooling and molecular thermal con- 
duction. Only above 160 km does molecular conduction balance EUV heating to 
maintain the observed thermospheric temperatures (see Figure 35). 

The 15-jim cooling arises mostly from transfer of atomic-0 kinetic energy to C0 2 
vibrational energy; the magnitude of the cooling rate depends upon the 0-C0 2 energy 
exchange rate and the relative amount of atomic-0 in the Venus day side thermosphere. 


HEATING RATES 



Fig. 35. Global average heating/cooling rates (Kday -1 ) for the Venus thermosphere. EUV heating 
efficiency - 10%. C0 2 -0 de-excitation rate - 8 x 10" 13 cm 3 s“ 1 at 300 K. Taken from Dickinson and 
Bougher (1986). Later NCAR NLTE model calculations (Figure 34) make use of revised EUV fluxes, 
stronger C0 2 15-jim cooling, and larger heating efficiencies (~14%), yielding the same basic thermal 

balances shown here. 
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A pedagogical discussion of translation-vibration energy transfer can be found in 
Yardley (1980). The 0-C0 2 excitation rate has not been measured at low temperatures, 
but may be large. In general the measured quantity that is commonly quoted is the rate 
coefficient for quenching or de-excitation (the reverse of reaction (69)). The rate 
coefficient for collisional excitation, k 69 , is related to that for quenching, k _ 69 , by 
expression 

69 = k. 69 x gexp(- AE/kT) , ( 70 ) 

where AE is the vibrational spacing and g is the degeneracy of the upper state, equal 
to u 2 + 1 for the doubly-degenerate C0 2 bending mode. In the high temperature regime, 
the temperature dependence of the probability, P, of de-excitation in a collision, to which 
the rate coefficient is proportional, is given by the Landau-Teller expression (Landau 
and Teller, 1936) 

P ~ exp( - T ~ ,/3 ) . ( 71 ) 

In the Landau-Teller model and its later refinement and extension to polyatomics, the 
SSH model (Schwarz etal., 1952), energy transfer results from a ‘direct’ interaction 
involving the exponential repulsive potential between the species. This strong tempera- 
ture dependence is characteristic of models involving ‘impulsive’ energy transfer of all 
types (e.g., McClelland etal., 1979). At low temperatures, the energy transfer 
probabilities are found to deviate substantially from the Landau-Teller values. The 
deviation is small for noble gases ; the probabilities are found to decrease with decreasing 
temperature over a large range of temperature, but level off at very low ( 1 50 K) tempera- 
tures (Inoue and Tsuchiya, 1975; Simpson etal., 1977). For some species, usually 
efficient quenchers, the deviation from the Landau-Teller model can be several orders 
of magnitude at low temperatures and, below a certain temperature, the probability of 
quenching can increase with decreasing temperature. For example, the rate coefficient 
for quenching of the v 3 mode of C0 2 by N 2 is large and exhibits a minimum at about 
1500 K (Moore etal., 1967; Rosser etal., 1969). This behavior for a near-resonant 
V —V transfer was successfully explained by Sharma and Brau (1969), who incor- 
porated the effect of long-range attractive multipole forces into the model for vibrational 
relaxation. Large rates and negative temperature dependence were also observed for 
de-activation of C0 2 , H 2 0, and D 2 0 at temperatures below 1000 K (Buchwald and 
Bauer, 1972). Significant departures from the Landau-Teller temperature dependence 
are also oibserved in de-activation of the asymmetric stretch mode of C0 2 by several 
molecular species at temperatures below about 400 K. (Bauer el al., 1987). Quenching 
of the bend-stretch manifold of C0 2 by H 2 is rapid, even though the vibrational modes 
have dissimilar energies. Allen etal. (1980) have reported a rate coefficient of 
5 x 10" 12 cm 3 s" 1 at 295 K; the rate coefficient increases with decreasing energy to 
7.5 x 10~ ,2 cm 3 s' 1 at 170 K. The corresponding quenching probabilities are 
7.4 x 10" 3 and 1.4 x 10" 2 . V - T, R transfer in other systems such as NO-NO has 
been reported to be rapid, with quenching probabilities on the order of 10" 3 and 
increasing with decreasing temperature below 400 K (Yardley, 1980). 
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Quenching of small molecules by atoms can be anomalously fast if the atom has an 
open shell (Chu et al. , 1983). H, Cl, and F quench C0 2 rapidly compared to noble gases, 
such as He, Ne, and Xe (Chu etal., 1983; Flynn and Weston, 1986). The large cross 
sections and the negative temperature dependences at low collision energies for such 
systems have been ascribed to the formation of long-lived collision complexes in the 
presence of strong attractive chemical, hydrogen bonding or Van der Waals forces 
(Ewing, 1978; Lin et al., 1979; Parmenter and Seaver, 1979; Gordon, 1981 ; McClelland 
et al., 1979). The lifetime of a collision complex is longer at low collision energies, and 
the longer-lived the collision complex, the greater the probability that the energy in the 
vibrational mode may flow into the weak bond of the complex and cause its dissociation. 
The V - T,R transfer probability for HF-HF, a system with strong hydrogen bonding, 
is about 1% at room temperature and increases with decreasing temperature (cf. 
Poulson et al., 1978). For C0 2 -H, the large rates even at high energies have been 
ascribed to the presence of a chemically reactive channel, CO + OH -» C0 2 + H, and 
to the formation of a stable intermediate radical, HC0 2 (Flynn and Weston, 1986). 

Atomic oxygen is a very efficient quencher of vibrational excitation, probably due to 
its reactivity, especially where de-activation can take place by atom exchange. The 
de-excitation rate coefficients for CO tr (the double dagger denotes vibrational exci- 
tation) and 0 2tt by O are about 9 x 10" 13 cm 3 s~ '; these rates are about 25 times 
faster than that for quenching of N 2Tt , where atom exchange is not possible (McNeal 
etal., 1974; Fernando and Smith, 1969; Eckstrom, 1973). Isotope exchange studies of 
NO with ,s O showed that the rate coefficient for O atom exchange at 298 K is 
3.7 x 10 " 1 1 cm 3 s ” 1 ; if the complex is assumed to break-up statistically, the implied 
rate of complex formation is 7.4 x 10” 11 cm 3 s~ l (Anderson etal., 1986). The rate 
coefficient for vibrational relaxation of NO by O is of the same order of magnitude, 
6.5 x 10 - " cm 3 s' 1 (Fernando and Smith, 1979). Isotope exchange of ,8 0 with 0 2 
proceeds with a rate coefficient of 2.9 x 10 " 12 cm 3 s* 1 and vibrational relaxation with 
a rate coefficient of 5.4 x 10 _ 12 cm 3 s'* (Anderson etal., 1986; Quack and Troe, 
1977). 

All of the measurements or calculations of the 0-C0 2 quenching rate have been at 
high temperatures where the Landau-Teller model might be expected to be valid. Center 
(1973) found O atoms an order of magnitude more efficient at relaxing the bending mode 
than At atoms; a probability of about 1% was measured at 3000 K. Bass (1974) 
computed the probabilities assuming classical scattering with a hard sphere interaction 
and quasi-diatomic model for C0 2 for energies from 1 to 10 eV and found that when 
the C atom was struck by the impinging O atom, the average fraction of energy 
transferred to the v 2 mode was 0.447. Schatz and Redmon (1981) used quasi-classical 
trajectory theory to compute the cross sections for excitation of the beznding mode and 
found a value of 1 .9 x 10 " 15 cm 2 at 2.2 eV, whereas the excitation of the symmetric and 
asymmetric stretch modes were 2 and 4—5 orders of magnitude smaller, respectively. 
Harvey (1982) obtained cross sections 50% smaller than Schatz and Redmon (1981) 
with a quantum mechanical vibrational close-coupling calculation using the same poten- 
tial surface. Extrapolation of these results to low temperatures using the Landau-Teller 
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expression would give very small quenching probabilities, but is not justified. The 
potential used by Schatz and Redmon (1981) and Harvey (1982) is a sum of pairwise 
exponential repulsion terms; no portion of the potential was attractive. At center-of- 
mass energies larger than 1 eV, the repulsive core of the potential dominates (Bass, 
1974), but at low energies, attractive forces are expected to become important. 

A value for the 0-C0 2tt quenching rate at low energies was derived by Sharraa and 
Nadile (1980) from measurements of the 15 pm radiation in the terrestrial atmosphere. 
They recommended a de-excitation rate coefficient of 6.67 x 10“ 10 ' '*T ° 5 cm 3 s" 1 . 
This rate has been refined by Wintersteiner et al. (1988) to 

1.5 x 10 ’ , 3 T 05 + 2.32 x 10“ 9 exp( - 76.75/T I/3 )cm 3 s“ 1 , (72) 

valid in the approximate temperature range 200-700 K. The first term in expression (72) 
is due to attractive forces, which dominate at low temperatures and the second is the 
Landau-Teller expression that accounts for the short-range repulsive forces. The poten- 
tial error in this expression, about a factor of 2 , are mostly due to uncertainties in the 
atomic O density (Wintersteiner, private communication). At 200 K, the temperature 
which prevails in the lower part of the Venus thermosphere, the rate coefficient would 
be 2 x 10 “ 12 cm 3 s “ \ about 1 % of gas kinetic. The rate adopted by Dickinson (1984) 
is less than 40% of the value in expression (72) and other recent models have assumed 
quenching coefficients less than 1 % of the Wintersteiner et al. value. Rates on the order 
of 1 x 10 - 12 cm 3 s“ 1 appear to be consistent with the terrestrial and martian heat 
budgets (Dickinson et al., 1987; Robl tet al, 1987, 1989; Bougherand Dickinson, 1988; 
Bougher et al ., 1990b). The reluctance to adopt the rate derived from terrestrial 15 pm 
observations is not justified. Certainly, however, a low temperature laboratory measure- 
ment of this rate coefficient is badly needed, not only for models of the Venus 
thermosphere, but for the Earth and Mars as well. 

The strong nonlinear (exponential) temperature dependence of the C0 2 15-pm 
cooling insures that the Venus exospheric temperature does not vary by more than 
-65K over the solar cycle (300-235 K) (Dickinson and Bougher, 1986) (see 
Section 7). In addition, the O densities themselves increase with increasing solar 
activity. Calculations show that the production rate of O above 100 km on Venus 
increases from 2.2 x 10 12 to 3.3 x 10 12 cm “ 2 s “ 1 from low to high solar activity (Fox, 
1989b). The VTS3 model predicts an increase of O from 2.7 to 7.3% at 135 km from 
^10.7 — 74 to F jq 7 = 200. This change in the O abundance will act to dampen the 
response of the thermospheric temperature to changing solar activity, increasing the 
cooling at higher solar fluxes. The combination of the effects of the exponential tempera- 
ture dependence of the 0-C0 2 collisional excitation rte and the change in the atomic 
O mixing ratio with solar activity provide the Venus thermosphere with an ‘atomic 
oxygen thermostat’, regulating C0 2 cooling, that is less effective at Mars where the O 
abundance is lower (Dickinson and Bougher, 1986; Bougher and Dickinson, 1988; Fox, 
1989b). 

The nonlinear temperature dependence of the 1 5 pm cooling also means that it cannot 
be adequately addressed by simple linear perturbation (Newtonian cooling) schemes 
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(Dickinson, 1973). Instead, a nonlinear temperature-dependent cooling parameteri- 
zation was developed by Bougher et al. (1986) for use in multi-dunensional model 
simulations. It first requires an exact 1 -D NLTE calculation of a representative reference 
lS-pm cooling and corresponding temperature profile. For Venus, since the day-to-night 
temperature variation is so large (~200 K), separate dayside and nightside mean 
reference cooling rate profiles were generated. These were later used to drive a nonlinear 
NLTE cool-to-space formulation for deviations from these references. The deviations 
depend on departures of total temperature and atomic-O concentrations from the 1-D 
model values. The separate reference cooling for dayside and nightside mean conditions 
is matched across the terminator by interpolation so that discontinuities are minimized. 
Such an interactive scheme is fast, yet it allows the accurate calculation of total 15-pm 
cooling at other temperatures at any local time over the globe. The expression for the 
cooling rate derived by Bougher (Equation (A.5) of Bougher et al. (1986)) gives approxi- 
mately the correct dependence of the IR cooling on temperature, pressure, and O 
concentration, and reduces to reference model cooling for composition and temperature 
profiles of either the day- or nightside reference models. Recent Mars TGCM model 
studies also confirm the usefulness of this NTLE cool-to-space formulation (Bougher 
etal., 1988c, 1990b). 

The one-dimensional model results of Dickinson and Bougher (1986) showed that 
dayside C0 2 cooling is dominated by the 15-pm fundamental above 130 km. Near-IR 
heating peaks at 130 km, with 2.7-pm bands important above 120 km, and <. 2.0-pm 
bands strongly contributing over 100 to 120 km. A calculation of the dayside v 2 = 1 
populations relative to LTE confirms that complete NLTE conditions hold above about 
140 km, while LTE is a good assumption only below 1 10 km. 

5.5. Molecular and eddy conduction 

Molecular conduction serves to transfer EUV heat deposited in the thermosphere 
downward to a lower level, the mesopause, where it can be radiated to space by IR active 
gases. On Venus, C0 2 15-pm emission dominates all others in providing this cooling 
near the bottom of the thermosphere. Temperatures are also affected by the amount of 
cooling eddy processes are assumed to produce (Hunten, 1974; Izakov, 1978). The 
turbulence responsible can heat the thermosphere through the dissipation of its energy; 
conversely, cooling can take place by eddy heat conduction in a process analogous to 
molecular conduction. The relative importance of turbulent heating or cooling in a given 
planetary thermosphere is still the subject of much debate (e.g., Hunten, 1974). 

Izakov (1978) argued theoretically that conditions exist when cooling of a planetary 
thermosphere by turbulence predominates, and conditions when heating prevails. For 
Venus, eddy mixing is estimated to provide cooling (upper limit) or a compensation of 
heating and cooling (lower limit). Hollenbach et al. (1985) approached the problem from 
a more empirical point of view, whereby the relative roles of heating and cooling were 
varied to match the observed temperatures over the solar cycle. The 1-D model calcu- 
lations of Dickinson ( 1 972, 1 976) and Dickinson and Bougher ( 1 986) have all addressed 
the case where vertical eddy mixing gives rise only to cooling. This establishes an upper 
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limit to the role of eddy mixing in contributing to the overall thermospheric heat budget 
(see Section 5.3). 

Questions still remain as to the complementary roles of eddy cooling and diffusion, 
particularly in the context of Venus circulation models. A few key points can be gleaned 
from recent modeling studies. The same eddy coefficient must be used for eddy diffusion 
and heat conduction in the absence of theoretical arguments to the contrary (Bougher 
etai , 1986, 1988b). However, the coefficient for eddy viscosity and diffusion may be 
different (Prandtl No. 1) (Mengel et al. y 1989). This implies that a self-consistent calcu- 
lation of winds, composition, and temperatures is required to address the impact of eddy 
mixing in the thermosphere. The practice of calculating temperatures using eddy heat 
conduction independent of the corresponding mixing effect on model densities is 
improper. Furthermore, observed densities appear to be influenced by a combination 
of global wind and eddy mixing processes; i.e., smaller eddy coefficients are needed in 
Venus circulation models than formerly implied by 1-D (non-dynamical) models (see 
Section 6.6). Thus, it appears that the introduction of eddy cooling is not generally a 
viable method for modeling thermospheric temperature profiles until all other dynami- 
cal, radiative cooling, and heating processes have been carefully examined. The effects 
of small-scale mixing, although difficult to model, need to be incorporated in a self-con- 
sistent fashion. No Venus model yet devised has this capability. 


6. Thermospheric Circulation and Transport 


6.1. Introduction 

The dynamics of Venus’s atmosphere can be separated into two distinct regimes in 
which very different flow patterns dominate: (1) a retrograde, zonal, superrotating flow 
from the surface up to ^ 70 km, and (2) a mean subsolar-to-antisolar (SS-AS) solar- 
driven flow above —95-100 km. In between (70-95 km), zonal winds are generally 
decreasing with altitude while thermospheric-type SS-AS winds are increasing. This is 
a complicated transition region in which changing zonal and meridional flow patterns 
may be superimposed on a possible lower return branch of the SS-AS circulation 
(Clancy and Muhleman, 1985b; Goldstein, 1989). A summary of these basic dynamical 
regions and the relevant supporting data is given in Figure 36 (Goldstein, 1989). 

Solar UV-IR heating clearly drives the SS-AS circulation. An eddy viscosity 
appears to be required to reduce pressure-driven winds toward acceptable values (see 
Section 6.4). A redistribution of angular momentum by a strong cloud-top Hadley 
circulation maintains the retrograde zonal jets observed at mid-latitudes (Hou, 1984; 
Hou and Goody, 1985; Walterscheid etal ., 1985). Recent studies support the idea, 
originally proposed by Fels and Lindzen (1974), that the superrotation at cloud tops is 
maintained by the pumping action of thermally driven tides (Pechmann and Ingersoll, 
1984; Hou, 1984; Hou and Goody, 1985). Baker and Leovy (1987) and Leovy (1987) 
suggest that the equatorial wind speed near Venus’s cloud top is maintained by a balance 
between the pumping effect of the semidiurnal tide and vertical advection by the Hadley 
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Fig. 36. Dynamical regions of the Venus mesosphere-thermosphere. Taken from Goldstein (1989). 


circulation in the thermal driving region. This rapid superrotation (at all latitudes) has 
not yet been adequately explained. 

Zonal circulation seems to exist throughout the atmosphere, although it does not 
dominate SS-AS flow above 100 km (Goldstein, 1989). Venus cyclostrophic balance 
models predict maximum zonal winds of 1 30 m s " 1 near 70 km. The observed poleward 
warming above this level implies decreasing zonal speeds with altitude (70-90 km), 
becoming negligible by 90-100 km. A concomitant increase in the mean meridional 
component with possible 120 ms* 1 speeds near 90-100 km has also been predicted 
(Taylor etal., 1980). Infrared heterodyne observations detect a weak zonal retrograde 
component (20-30 ms -1 ) (equatorial) at 110+ 10 km, superimposed on strong 
100-140 ms -1 (SS-AS) cross-terminator winds (Goldstein, 1989). Finally, a slightly 
larger component of retrograde zonal winds ( ~ 50-60 m s * ') is implied above 1 50 km 
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from helium and night airglow observations (Mavr et al., 1980; Stewart et al., 1980). See 
Section 6.5 for details of superrotation mechanisms and sensitivity studies. 

There are, however, no direct or indirect measurements of wind velocity in Venus’s 
upper atmosphere above 110 km. The circulation of most of the thermosphere must, 
therefore, be inferred from observations of its temperature and density structure (Mayr 
era/., 1980; Schubert et al., 1983; von Zahn era/., 1983) and subsequent modeling 
efforts (Bougher, 1985; Mayr era/., 1985; Bougher era/., 1986, 1988b, 1989; Mengel 
era/., 1989). 

6.2. Thermospheric circulation: symmetris vs zonal 

The most prominant feature of the thermosphere which suggests a strong SS-AS wind 
system is the large contrast in temperatures and densities between day and night. 
Temperatures on the dayside of Venus increase from ~ 170-180 K at the ‘mesopause 
(100km) level (Taylor era/., 1979, 1980; Seiff era/., 1980) to near 300-310 K above 
150km for solar maximum conditions (von Zahn era/., 1979, 1980; Niemann era/., 
1979, 1980b; Hedin era/., 1983; Keating era/., 1980, 1985). The rise of temperatures 
above a 100 km mesopause is similar to that of the Earth’s thermosphere, although the 
temperatures are much colder (see Figure 31). Nightside temperatures are shown to 
decrease from 170 K at the mesopause to values as low as 100-130 K above 150 km 
(Keating et al., 1979, 1980; Niemann et al., 1979, 1980b; Seiff etal., 1980; Seiff, 1982). 
This vertical structure is remarkably unlike the terrestrial nightside thermosphere, and 
has been termed a ‘cryosphere’ (Schubert et al., 1980). In addition, changes in tempera- 
ture across the terminator are very abrupt, with the minimum value observed just beyond 
mi dnigh t at 02:00 LT (Mayr etal., 1980; Keating etal., 1980). 

These thermospheric (above 150 km) kinetic temperatures are inferred primarily from 
the observed scale heights of the measured neutral densities. The diurnal variation of 
the heaviest species (C0 2 , CO, N 2 , and O) show a near symmetry about the subsolar-to- 
anti-solar axis (see Figure 37), with a density maximum at a fixed altitude at local noon 
and minimum at midnight (von Zahn et al., 1983). The slow planetary rotation coupled 
with dayside solar heating results in this expansion of the dayside thermosphere; i.e., 
the height of a given pressure surface rises (falls) with increasing (decreasing) tempera- 
tures. Nightside densities are strongly depleted with respect to dayside values, suggesting 
a significant contraction of the nightside thermosphere with constant pressure surfaces 
decreasing abruptly in altitude across the terminators. Helium and hydrogen, the two 
lightest species, show a diurnal density reversal, with densities higher at night than 
during the day (Niemann et al., 1980b; Brinton et al., 1980). This signature is consistent 
with their small atomic weights and correspondingly large-scale heights. These 
properties increase the transport efficiencies by the large scale winds for these species, 
just as in the atmosphere of the Earth (Mayr et al., 1978). The diurnal variation ofhelium 
at 165 km reveals a distinct (lasting) nightside bulge (see Figures 38(a) and 38(b)), with 
a maximum at 03 : 00 LT roughly 30-40 times that on the dayside (Taylor et al., 1984). 
Results obtained from Brinton etal. (1980) suggest a similar night-day bulge of hydro- 
gen peaking at 05 : 00 LT, with a diurnal ratio of nearly 400 : 1 at 165 km (sec Figure 25). 
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LOCAL SOLAR TIME (hr) 

Fig. 37. Measurements of C0 2 and O densities at 170 km altitude taken by the Pioneer Venus orbiter 
neutral mass spectrometer (ONMS) over nearly three diurnal cycles. Taken from vonZahn ttal 

(1983). 

These bulges imply a ‘wind-induced diffusion’ in which light constituents are preferen- 
tially transported from day-to-night, with little returning to the dayside along the return 
branch of the circulation (Mayr et al ., 1978, 1980). 

This temperature difference between day and night gives rise to horizontal pressure 
gradients, which to first order should drive a strong subsolar-to-anti-solar (SS-AS) 
circulation in Venus’s upper atmosphere (Dickinson and Ridley, 1977; Mayr etal . , 
1980; SeifT, 1982). A simplified flow pattern would be axisymmetric about the 
Sun- Venus line, with dayside upwelling centered on the subsolar point, strong cross-ter- 
minator flow, and subsidence centered on the antisolar point (Schubert era/., 1983). 


SPAC 191 97 



98 


J. L. FOX AND S. W. BOUGHER 



LOCAL SOLAR TIME (hr) 

Fig. 38. Measurements of He densities at 170 km altitude taken by the Pioneer Venus orbiter neutral mass 
spectrometer (ONMS) over nearly three diurnal cycles. Helium has a dawn bulge analogous to that tn 
hydrogen, with some differences in shape. Taken from von Zahn el at. (1983). 

Return flow, possibly between 70-90 km (Goldstein, 1989), would complete the circuit. 
Streamlines for transported species approximately follow constant pressure surfaces, 
which descend in altitude from day to night (Seiff, 1982). Supersonic velocities are 
estimated across the terminators for laminar flow conditions ; turbulent viscosity is likely 
important in modifying the pressure-driven wind speeds (Seiff, 1982; Bougher el at., 
1986). This basic day-night flow pattern was hypothesized and studied in a series of 
numerical models by Dickinson and Ridley (1972, 1975, 1977) prior to Pioneer Venus 
measurements, and by Bougher et at. (1986, 1988b, 1990a), Mayrer al. (1980, 1985), and 
Mengel et al. (1989) afterwards. 
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This single symmetric cell* driven by solar UV-IR heat* must occur. However, the 
actual motions in Venus’s thermosphere are likely a superposition of this SS-AS flow 
and a retrograde zonal wind component. Post-midnight maxima in the helium and 
hydrogen densities (Niemann etal . , 1979; Brinton etal .* 1980; Taylor etal ., 1984)* 
displacement of the minimum diurnal temperature to 02:00 LT (Mayr etal., 1980; 
Keating etal. , 1980), and the NO night airglow maxima at 02:00LT (Stewart etal., 
1980; Bougher etal., 1990a) all suggest a westward superrotation of the upper 
atmosphere (£ 150 km) with zonal winds speeds of 50-100 ms" 1 . This corresponds 
to a rotational period of 4-8 Earth days. This modified flow pattern now includes 
upwelling (divergence) near the subsolar point, and subsidence (convergence) offset 
from midnight, yielding a nightside flow which stagnates near 02:00-03:00 LT (see 
Figure 39). The zonal winds add to the SS-AS flow across the evening terminator* 
giving stronger winds there than at dawn. 


ROTATION 





Fig. 39. A schematic flow pattern for the nightside thermosphere of Venus. Taken from Niemann et at. 

(1980). 

Models of the Venus thermospheric circulation have been under development for 
nearly 18 years. Experience gained from Earth thermosphere studies initially guided the 
model processes and parameterizations considered. Prior to Pioneer Venus, only the 
SS-AS circulation pattern was generally thought to be operating in the Venus 
thermosphere. The Dickinson and Ridley (1972, 1975, 1977) two-dimensional model 
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(DRM) of the temperature, winds, and composition of the Venus thermosphere, 
correctly predicted the gross characteristics of this SS-AS mean circulation. Axial 
symmetry was assumed in choosing altitude and solar zenith angle (SZA) as the two 
model coordiantes, in accord with the dominant SS-AS flow. Later versions of this 
model were fully self-consistent and used linear perturbation theory to simultaneously 
calculate thermal structure, winds, and CO,, CO, and O densities above 100 km. 
Detailed results of the model, however, differed considerably from the Pioneer Venus 
observations (Schubert etal., 1983; Dickinson and Bougher, 1986). Figure 40 shows 


COMPARISON OF DRM (77) PREDICTIONS AND P.V. OBSERVATIONS 
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Fig. 40. Comparisons of Dickinson and Ridley (1977)(DRM) model predictions and Pioneer Venus (PVO) 
observations. Solid lines are for the 2-D model simulations; dashed lines are from PVO observations as given 
by the empirical model of Hedin etal. (1983). Taken from Dickinson and Bougher (1986). 


that the Dickinson-Ridley temperatures were somewhat too warm on the dayside and 
much too warm, by as much as 100 K, on the nightside. Their calculations failed to 
predict the nightside cryosphere entirely. Model atomic -O and CO bulges on the 
nightside were also predicted by the DRM, that were not observed. Correspondingly, 
the dayside O-mixing ratios at the Fl-peak (140 km) were calculated to be 2.5-5.6/ 0 , 
much smaller than actually observed (17-20%) (von Zahn et ai, 1983). No shift in the 
nighttime minimum temperatures from midnight is possible using the axially-symmetric 
framework assumed in the DRM. Effects of zonal flow require a three-dimensional 
model simulation (Mayr etal., 1980; Bougher, 1985; Bougher etal., 1988b). 


SPAC 191/100 


ORIGINAL PAGE IS 
OF POOR QUALITY 







STRUCTURE, LUMINOSITY, AND DYNAMICS OF THE VENUS THERMOSPHERE 


101 


Several explanations were initially forwarded to account for these discrepancies 
between DRM and the data. Two groups, in particular, focused their attention on the 
Pioneer Venus neutral data and its global modeling. The Goddard group (cf. Mayr et al., 
1980, 1985) used a 3-D spectral model to simulate the observed day— night distribution 
of composition and temperatures. Dynamical properties of the Venus theermosphere 
were initially inferred using a simplified linear spectral model and Pioneer Venus ONMS 
measurements (Niemann etal., 1979, 1980; Mayr et ai, 1980). A combination of 
parametrically prescribed strong vertical eddy diffusion (K £ 3 x 10 7 cm 2 s _1 ) and 
retrograde zonal winds (5—10 day period) seemed to modify the fundamental diurnal 
tide sufficiently to yield a reasonable match between observed and modeled phases and 
magnitudes of global temperatures and densities. Strong vertical eddy diffusion was 
required to counteract the mass-transport effects of the large-scale winds. Non-linear 
processes (i.e., wind-induced diffusion), important for helium and hydrogen, were 
estimated to steepen the nighttime density maxima and shift them significantly toward 
dawn (0200—0400 LT). Maximum horizontal winds of — 200 ms -1 were inferred from 
the day-night density (temperature) variations observed by the ONMS. Strong vertical 
eddy diffusion, retrograde zonal winds, and reasonable day— night winds were all sug- 
gested to improve the DRM simulations. 

Further work was described by Mayr etal. (1985), whose spectral model is more 
realistic in that it accounts for nonlinear processes, includes higher order tidal com- 
ponents, and describes the major gases in self-consistent form. Eddy diffusion and 
conduction were modified to vary inversely as the square root of the number density 
(after von Zahn et al., 1980), with a maximum eddy coefficient of 6 x 10 7 cm 2 s' '. The 
EUV heating efficiency was taken to be 30%. The basic diurnal contrasts in heavy (C0 2 , 
CO, and O) and light (He) species and temperatures observed by PV were reproduced. 
Also, superrotating zonal winds (6-day period) were found to have little impact on the 
diurnal distributions of the heavy species and temperatures. However, a significant 
dayside depletion and nightside enhancement in O and CO densities (with respect to 
Hedin etal. (1983) values) was obtained. This occurs because the calculated solar- 
driven thermospheric wind system, with maximum horizontal wind speed of about 
300 m s' 1 , is too strong. The Mayr etal. (1985) lower thermosphere wind speeds 
( ~ 1 10 km s ” ') are not consistent with ground-based measurements (Goldstein, 1989). 

The group at NCAR updated and improved the DRM code, and subsequently 
adopted the finite-difference approach of the terrestrial Thermospheric General Circu- 
lation Model (TGCM) (cf. Dickinson et al., 1984) modified for Venus. Bougher (1985) 
and Bougher et al. (1986) proposed that a global circulation system weaker than the 
DRM originally predicted would enable increased dayside O and CO densities to be 
maintained, while simultaneously reducing the buildup on the nightside. Furthermore, 
since nightside heating is maintained primarily from the circulation (i.e., adiabatic 
compressional heating), slower winds would also result in cooler nightside temperatures 
more indicative of a cryosphere. Suggestions along this line were also made by von Zahn 
et al. (1983) who noted that dayside 0/C0 2 ratios increased for the series of Dickinson 
and Ridley (1972, 1975, 1977) models, due to progressively weaker day-to-night winds. 
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Schubert et ai (1983) also noted that dynamical processes may be required to explain 
the very low nightside thermospheric temperatures. Specifically, low temperatures must 
be maintained in the face of (1) advection and convection of dayside heat across the 
terminator, (2) compressional heating due to nightside subsidence, and (3) upward 
molecular heat conduction from the relatively warm 100 km level. 

This weakened circulation hypothesis seems to give the most self-consistent 
explanation for all the available Pioneer Venus composition, temperature, and airglow 
data. An updated version of the DRM model was constructed by Bougher et ai (1986) 
to examine how SS-AS winds, eddy diffusion/conduction, and strong C0 2 15-jim 
cooling affect the day-night contrasts in Venus densities and temperatures. It was 
learned that symmetrized empirical model fields of Hedin etal. (1983) could be 
simulated largely by adding a wave-drag mechanism, resulting from turbulent dissi- 
pation effects on the horizontal winds, which weakened the day-to-night flow. Eddy 
viscosity and Rayleigh friction were parameterized within the horizontal momentum 
equation to mimic the wave-mean-flow interaction responsible (see Section 6.4). Maxi- 
mum terminator winds of ^ 230 ms -1 were found to be consistent with observed global 
composition and temperature fields. Such a reduction of winds by nearly a factor of two 
(from those of DRM) permitted less O and CO to be carried to the nightside. The 
dayside FI -peak O-mixing ratios were calculated to be ~ 15%. The increased isolation 
of the day and nightside thermospheres also serves to enhance the diurnal temperature 
contrast, giving midnight exospheric temperatures of about 115 IC, characteristic of the 
cold temperatures inferred from density data. This occurs because slower winds support 
less adiabatic compressional heating, the dominant source of nightside heating 
(Dickinson and Ridley, 1977). It was also shown that drag could not be used simul- 
taneously with strong eddy conduction to derive the cold nightside temperatures. Slower 
winds are sufficient to maintain cold nightside temperatures, and are also consistent 
with the observed diurnal density contrasts. Winds calculated near 110 km match 
Earth-based measured values (100-140 ms" 1 ) quite well (Goldstein, 1989). Figures 
41(a-d) illustrate contour plots from the standard NCAR 2-D model simulamtion 
(Bougher et ai , 1988) for solar maximum conditions (F 10 7 = 200). Notice that constant 
pressure surfaces decrease in altitude from the day to the nightside, in agreement with 
the ‘collapsed* nightside cryosphere observed. 

Superrotation effects on thermospheric densities and temperatures can only be 
examined in a three-dimensional coupled chemical-dynamical model. A fully self-con- 
sistent code coupling large-scale and sub-grid scale effects, including wave-drag and 
superrotation, is very difficult to formulate. A prelude to such a code was developed by 
Bougher et ai (1988b), based on the NCAR terrestrial TGCM. The general framework 
of this Earth code (Dickinson et ai , 1984) was used in concert with new inputs and 
physical parameterizations to examine Venus thermospheric processes. The Venus 
TGCM (VTGCM) is internally self-consistent in that it calculates global distributions 
of C0 2 , CO, and O that are consistent with the three-dimensional day-night tempera- 
ture contrasts and corresponding large-scale winds. The model covers a 5° by 5° 
latitude-longitude grid, with 24 constant log-pressure levels in the vertical, extending 
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SOLAR ZENITH ANGLE SOLAR ZENITH ANGLE 

Fig. 41. NCAR 2-D model contour plots. AU frames are presented on an altitude scale for convenience. 
The sharp upper level cutoff (r - 7) corresponds to the top boundary of the model. Solar maximum 
simulation: (a) z - log r (p 0 /p), (b) total temperature (K), (c) horizontal velocity in ms"\ (d) vertical 
velocity in cm s' . Solar minimum simulation: (e)z - log ^PolPh (0 total temperature (K). Taken from 

Bougher el ai (1986). 
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from approximately 105 to 200 km at local noon. Dayside O and CO sources arise from 
C0 2 net dissociation; sinks are provided primarily by advection and convection of the 
large-scale winds. Total temperature is derived as the sum of an adopted global mean 
reference temperature profile (Dickinson and Bougher, 1986) and calculated pertur- 
bation temperatures. Sub-grid scale processes (i.e., eddy diffusion, viscosity, con- 
duction, and wave-drag) are not self-consistently incorporated, but rather parameterized 
using standard aeronomical formulations (cf. Hunten, 1974; vonZahn etal., 1980). 
Rayleigh friction is used to mimic wave-drag effects on the mean flow (see Section 6.4); 
its incorporation is analogous to early terrestrial mesopause studies which sought to 
define an approximate time-scale for the required wave-breaking process (Schoeberl and 
Strobe!, 1978). 

A symmetric form of the VTGCM was first used to examine the SS-AS circulation 
on Venus for solar maximum conditions (Bougher et al., 1988b). Inputs and parameters 
comparable to the NCAR 2-D model of Bougher etal. (1986) were incorporated, 
including a low 10% EUV heating efficiency, an 0-C0 2 de-excitation rate coefficient 
of 8 x 10 13 era 3 s 1 at 300 K, a maximum nightside eddy coefficient of 
1 x 10 7 cm 2 s~ ', and minimal dayside eddy diffusion and cooling 

(K£5 x 10 6 cm 2 s~ '). Prescribed wave-drag (Raleigh friction coefficient <> 10 s~ ') 
was primarily responsible for the relatively weak winds calculated (^ 230 m s " '), which 
is consistent with the observed day— night contrast in calculated densities and tempera- 
tures. Model exospheric temperatures ranged from 309 K (day) to 136 K (night), simu- 
lating the nightside cryosphere seen by Pioneer Venus instruments. 

Once the VTGCM was validated (in its symmetric form) by comparison with the 
NCAR 2-D code, further tests were conducted to examine the asymmetric character of 
the circulation by incorporating prescribed superrotating retrograde zonal winds. The 
adopted zonal wind profile was specified to increase from zero to 60 m s - 1 over 105 
to 140 km. Such a profile is consistent with recent suggestions that zonal flow near 
100-110 km is quite weak (Clancy and Muhlemann, 1985; Goldstein, 1989), while zonal 
flow above 140 km is significant (Mayr et al., 1980, 1985). Results show that the major 
(heavy) species and temperatures are not greatly affected by this superrotation. How- 
ever, the exospheric temperature minimum shifts to 02 : 00 LT, which is consistent with 
data and with the convergence of the horizontal winds after midnight. Figures 42(a) and 
42(b) show pressure levels slices (at an average of 162 km) of VTGCM total tempera- 
tures, horizontal winds, and vertical winds. Also, the VTGCM nightside temperatures 
appear to be larger in the exosphere than in the lower thermosphere (see Figure 43); the 
NCAR 2-D model, formulated explicity for total temperatures, does not show this 
behavior. This discrepancy may be due to the difficulty of generating a representative 
global mean temperature profile for use with VTGCM perturbation values to give total 
temperatures (Bougher et al., 1988b). 

A more sensitive tracer of the thermospheric winds is required to further constrain 
the circulation. The incorporation of minor (light) species, such as helium, can provide 
such a tracer of the SS-AS winds and the superimposed zonal wind component 
(Bougher et al., 1988b; MengeUr al., 1989). The magnitude of the nightside helium bulge 


SPAC 191/104 



STRUCTURE, LUMINOSITY, AND DYNAMICS OF THE VENUS THERMOSPHERE 


105 


(a) UTIME= 0: 0 T-(U.V) 



(b) UTIME-0: 0 FIELD-W 


Z- 5.0 



LONGITUDE 


Fig. 42. VTGCM superrouting model pressure level slices (local time vs latitude) at z ■ 5 (or roughly 
162 km average altitude): (a) T + ( u , v )- superimposed total temperature and horizontal wind vectors, (b) tv 
- vertical velocity in cm s " The maximum horizontal wind vector in (a) corresponds to 230 ms' 1 . Taken 

from Bougher etai (1988b). 


should be very sensitive to the overall strength of the circulation and the amount of 
global eddy diffusion. Likewise, the local-time position of the bulge is a measure of the 
degree of asymmetry of the otherwise symmetric SS-AS flow. Mengel et al. (1989) 
presented calculations of helium that were used to fine tune the Mayr et al. (1985) 3-D 
spectral model circulation and eddy diffusion formulation. This improved model 


SPAC 191/105 




106 


J. L. FOX AND S. W. BOLGHER 



VIRA 

VTS3 

VTGCM — • — 

BOUGHER et al., 1986 • — 

Fig. 43. VTGCM symmetric model (a) noon and (b) midnight temperature profiles above 100 km. 
Calculated profiles are compared to the corresponding empirical model profiles of Hedin et al. (1983) (VTS3) 
and Keating et al. (1985) (VIRA). The NCAR 2-D model profiles of Bougher et al. (1986) are also shown. 

Taken from Bougher et al. (1988b), 


accounted for nonlinear processes and superrotation through a solid-shell approxima- 
tion. They assumed a 6-day superrotation period and used the EUV heating rates of Fox 
(1988) with an efficiency close to 20%. An altitude dependent Rayleigh friction scheme 
(see Section 6.4) was also adopted to weaken the pressure driven winds. 

Standard spectral-model profiles for noon, midnight, and global mean temperatures 
are illustrated in Figure 44. The global mean exospheric temperature of 220 K is consis- 
tent with data, and lies between 325 K (noon) and 155 K (midnight) model values for 
extreme solar maximum conditions. Superrotation was found to have relatively little 
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TEMPERATURE PROFILE 



solar majumum conditions (FI 0.7 * 238), Taken from Mengel ft ai (1989). 


effect on the temperature contrast between day and night (see Figure 45(a)); however, 
it did produce a shift in the temperature minimum to 02 : 00 LT, in accord with observa- 
tions. This dawn-dusk asymmetry in the flow also enabled the nightside helium bulge 
to peak at 05:00 LT (see Figure 45(b)). Its night-day ratio at 170 km was calculated 
to be 45 : 1, in good agreement with the 30 : 1 value seen at 165 km over the three diurnal 
periods analyzed by Taylor el al. (1984). Results show that the longer the superrotation 
period, the larger is the day-to-night buildup in the helium density, and the shorter is 
the time delay in the peak density after midnight. Figure 46 schematically shows the 
combined effects of the SS-AS flow plus the superrotating winds on the diurnal He 
density distribution. Both the horizontal and the vertical He distributions are also 
strongly dependent on the eddy diffusion coefficient. The data are best fitted by a 
SS-AS circulation of ~ 100 m s~ \ a moderate height-dependent global mean eddy 
diffusion (K <> 3 x 10 7 cm 2 s " '), and a superrotation period of 6 Earth days. In this 
model, the Prandtl number, which is defined as the ratio of the kinematic viscosity to 
the thermal diffusivity, and indicates the relative efficiency of mass and heat transport 
by turbulence, was assumed to be unity. Strong superrotating zonal winds appear to 
reduce (smear out) the magnitude of the nightside helium bulge from what SS-AS 
wind-induced diffusion would otherwise produce. 
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Fig. 46. Schematic illustrating the effect of atmospheric superrotation on the magnitude and phase of the 
diurnal variations of helium in the Venus thermosphere. Taken from Mengel etai (1989). 

An additional tracer of the thermospheric circulation is given by the observed NO 
nightglow, which is shown to be produced by radiative recombination of N and O atoms 
transported from their source on the day side to the nightside (Bougher etai , 1990a). 
This work, described in Section 3.10, shows how the maintenance of this nightglow 
constrains the magnitude of the SS-AS winds 150 m s' *) and global mean eddy 
diffusion (K <> 2 x 10 7 cm 2 s” l ) in the face of moderate zonal winds 55 m s " ') 
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having a period of 8 days. This nightglow provides a means to fine tune the VTGCM 
circulation presented by Bougher etai (1988b), 

It is apparent that differences in the mean state of the solar maximum Venus 
thermospheric circulation (and small-scale processes) still exist in the most recent model 
simulations (Bougher etai y 1990a; Mengel etai y 1989). However, the major tunable 
parameters have been identified: ( 1 ) zonal superrotation period, (2) the magnitude of the 
maximum global mean eddy diffusion coefficient, and (3) the magnitude of the peak 
Rayleigh friction coefficient (see Section 6.4) which is used to modify SS-AS wind 
speeds. A reasonable range of peak SS-AS wind speeds is - 100-200 ms" 1 , yielding 
a close match of model fields to available temperature, density, and airglow data. 

6.3. Departures from radiative equilibrium: adiabatic effects 

Radiative heating and cooling processes can be augmented by dynamical effects in 
planetary atmospheres, thereby maintaining temperatures that depart from radiative 
equilibrium values. This is particularly the case for Venus’s nightside, in which a given 
planet-fixed meridian experiences darkness for as much as 1 17 Earth days. The night- 
side temperatures would be expected to be extremely cold without some sort of day-to- 
night heat redistribution. The NCAR dynamical models of Dickinson and Ridley (1977) 
and Bougher et ai (1986, 1988b) have all presented heat budgets thought to maintain 
day and night temperatures. Little difference in mean dayside heat balances is estimated 
between 1-D and 3-D models; i.e., adiabatic cooling (due to upwelling/divergent flow) 
is negligible (Bougher et ai y 1986). 

Thermospheric winds, however, have a controlling effect on the maintenance of 
nightside temperatures. Weakened SS-AS winds have produced a significant change 
in the nightside heat budget that has enabled nightside cryospheric temperatures to be 
explained (Bougher etai y 1986, 1988b). Figure 47 shows heating and cooling rates at 
120 SZA for both the NCAR 2-D and VTGCM model calculations. Notice that 
adiabatic heating (from subsiding/convergent winds) largely balances C0 2 15-pm 
cooling below 140 km. Above, molecular conduction is offset by this same com- 
pressional heating. The magnitude of this heating is reduced by about a factor of three 
from the Dickinson and Ridley (1977) model, permitting the observed cold nightside 
temperatures (Niemann et ai y 1979, 1980b; Keating et ai y 1980) to be reproduced. Early 
ideas that invoked nightside eddy heat conduction to give cold nightside temperatures 
(cf. Niemann etai y 1980b; Schubert etai y 1983; Gordiets and Kulikov, 1985) appear 
to be inappropriate. Observed nightside temperatures (and the diurnal density con- 
trasts) can be maintained by weakened thermospheric winds alone (Bougher et ai y 
1986). 

6.4. Drag processes 

It is evident by the observed day-night density and temperature contrasts that some type 
of deceleration mechanism is necessary to slow the otherwise pressure-driven neutral 
winds (Seiff, 1982; Bougher etai , 1986). Early in the Pioneer Venus mission, thought 
was given to ion-drag and its possible influence on neutral winds (Mayr etai , 1980). 
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COMPARISON OF MODEL THERMAL BALANCES 



Fig. 47. A comparison of nighuide mean heat balances for the NCAR VTGCM and 2-D models at 
120 SZA. Units in K Earth day ~ ! . (a) VTGCM symmetric model heat balances at the equator, and (b) 2-D 
model heat balances. Taken from Bougher et al. (1988b). 


Since Venus does not have a significant intrinsic magnetic field (Russell etal ., 1979), 
ions are not generally confined to field lines, thereby eliminating a ‘grid* for neutrals to 
encounter. Rather, large horizontal ion-drift velocities, measured by the PV ion mass 
spectrometer (Taylor et al ., 1980) and the retarding potential analyzer (Knudsen et al 
1980) (on the order of ^ 1 km s" 1 below 200 km), may carry along neutrals to the 
nightside. This effect was estimated to be negligible since ion-neutral collision fre- 
quencies below 200 km are small (Mayr etal , 1980). The situation may be different 
during times of high solar wind dynamic pressure, when the region below the minimum 
ionopause ( ~ 220 km) may be penetrated by interplanetary magnetic fields (Luhmann 
et al., 1987). Under these conditions, temporary fields may generate a possible ion-neu- 
tral drag for ion velocities in excess of neutral velocities. A definitive analysis awaits a 
3-D coupled ionosphere/thermosphere calculation of global winds. 

Alternatively, there is ample evidence for turbulence or wave perturbations in the 
upper atmosphere of Venus (Seiff et a/., 1980; Seiff, 1982; von Zahn etal , , 1983; 
Kasprzak etal,, 1988; Mayr etal., 1988). All available observations, however, are 
snapshots of wave structure; horizontal (vertical) wavelengths in the range of 
100-600 km (5-15 km) have been observed in the major neutral densities. There seems 
to be more wave activity during nighttime than during the day, with considerable activity 
occurring in the vicinity of the predawn and postdusk terminators. Comparison with the 
model of Mayr etal. (1988) shows that the observed perturbations can be interpreted 
as being due to gravity waves propagating upward from source regions in the lower 
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thermosphere (130 km), or below 80 km. Most important, the data imply that significant 
wave energy arrives in the thermosphere from lower altitudes on a continuing basis. The 
effect of such waves on the thermospheric densities, temperatures, and winds is signifi- 
cant. Specifically, nightside neutral densities vary by as much as + 50% from day-to- 
day, likely due to upward propagating waves (Hedin et al. % 1983 ; Kasprzak et at . , 1988). 

Eddy conduction, eddy or wave drag, and eddy diffusion are typically parameterized 
in planetary models to mimic turbulent effects on temperatures, densities, and winds. 
For the Earth, it has been shown that the heat and momentum budgets near the 
mesopause cannot both be simultaneously satisfied without a strong zonal drag force 
(Lindzen, 1981; Holton, 1982, 1983). Such a frictional force is required to simulate the 
reversal of the latitudinal temperature gradient observed at the mesopause, with the 
lowest temperatures at the summer pole (due to adiabatic cooling from upwelling) and 
the wannest at the winter pole (due to adiabatic heating from subsidence). A modified 
meridional circulation is responsible. The underlying mechanism has long been thought 
to result from the unstable breakdown of tides or gravity waves near the terrestrial 
mesopause. Waves generated in the troposphere are assumed to attain sufficient ampli- 
tudes in the mesosphere to produce superadiabatic lapse rates, giving rise to convective 
overturning and turbulence. This causes attenuation of further growth and a conver- 
gence of the vertical momentum flux, resulting in a net zonal drag force (Lindzen, 1981; 
Holton, 1982, 1983). 

Rayleigh friction is an ad hoc scheme commonly used to mimi c this ‘wave-drag’ force 
in-lieu of self-consistently calculating momentum stress terms based on wave properties 
and the large-scale flow. Specifically, Rayleigh friction is formulated as a damping term 
added to the horizontal equation of motion which is proportional to the horizontal 
velocity. The constant of proportionality is assigned a fixed peak value (s' ') at some 
reference level, and prescribed to vary as a function of altitude (or pressure) below. The 
shape of the profile is usually chosen to be that of the corresponding eddy coefficient 
used in the model study. 

Models of this terrestrial wave drag progressed from crude Rayleigh friction schemes 
(Schoeberi and Strobel, 1978), to current wave drag models patterned after Lindzen 
(1981) and Holton (1982), which describe self-consistent methods for parameterizing 
the stress and diffusion due to gravity wave breaking in terms of the mean zonal wind, 
static stability, and specific wave parameters. A height-dependent Rayleigh friction is 
a first choice when little is known about wave parameters and their source. The initial 
approach for Venus has been guided by these terrestrial mesopause wave-drag schemes 
and the underlying wave-breaking mechanism as well. Bougher etai (1986) assumed 
that small-scale gravity waves arc responsible for the observed wavelike features in the 
altitude profiles of temperature and density (Seiff etai , 1980; Schubert, 1983). A 
deceleration of the SS-AS winds is proposed to become important at a level where 
wave-breaking maximizes; i.e., also where eddy diffusion peaks on the dayside 
(135-140 km). A profile for drag is prescribed that increases with height up to this level, 
consistent with a vertical propagating spectrum of waves each of which breaks at a 
different level. A Rayleigh friction time-scale of ^1 x 10 ~ 4 s “ 1 is specified. This leaves 
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SS-AS winds below 1 10 km largely unaffected, while those near 150 km are reduced 
from Dickinson and Ridley (1977) values by a factor to two. 

The day-night contrasts in observed Venus thermospheric temperatures and den- 
sities are closely matched by Bougher et ai (1986) using this wave-drag mechanism (see 
Section 6.2). It is, however, somewhat unsatisfying that an ad hoc formulation, rather 
than a self-consistent scheme, for wave-induced drag was used. Nevertheless, Rayleigh 
friction is the best that can be used until further wave parameters (and measured 
thermospheric winds) are available. The Bougher et ai (1986, 1988b) and Mengel et ai 
(1989) benchmark models, all of which rely on Rayleigh friction, serve as starting points 
for later studies. 

6.5. Thermospheric superrotation 

The observed asymmetries in the helium (and hydrogen) densities, as well as the 
observed NO nightglow intensity (Niemann et ai , 1980b; Brinton etai , 1980; Stewart 
et ai, 1980), suggest that the Venus thermosphere superrotates (above 150 km) with a 
period of 5 to 10 days (Mayr et ai, 1980, 1985). The origin of this superrotation is not 
known. Thermospheric superrotation could be a remnant of the cloud-top 4-day zonal 
winds observed in the lower atmosphere (cf. Schubert and Walterscheid, 1984). A 
decrease in the superrotation rate between the cloud tops (4-day period) and the 
thermosphere is related to the temperature reversal observed above 70 km (where it 
increases toward the poles). In principle, thermospheric superrotation could also be 
generated in situ , presumably through zonally symmetric pressure gradients producing 
a geostrophic or cyclostrophic balance or by nonlinear interactions involving the diurnal 
tides (Mengel etai , 1989). 

All circulation model studies to date have chosen to circumvent this question of the 
origin of this zonal flow, in order to examine its effects on the thermospheric structure 
(Mayr etai , 1985; Bougher etai , 1988b; Mengel etai., 1989). Mayr etai (1985) and 
Mengel etai (1989) both specify a solid-body rotation in which SS-AS winds are 
calculated with respect to a rotating reference frame. This has the effect of adding zonal 
wind speeds of ~ 75 m s “ 1 at the equator. A considerable shift in the nightside helium 
bulge (to 05:00 LT) is produced. However, any altitude dependence of a zonal wind 
profile cannot be given by assuming solid-body rotation. Bougher etai (1988b) pre- 
scribed an altitude dependent zonal wind profile increasing from zero to 60 m $ “ 1 over 
105 to 140 km. This vertical variation is consistent with weak zonal winds at 
100— 1 10 km (Clancy and Muhlemann, 1985b; Goldstein, 1989) and strong zonal winds 
above 150 km (Mayr et ai , 1980). Since tracer species are thought to generally follow 
constant pressure surfaces (streamlines) in their path from the day to the nightside, zonal 
wind shear could be important in modifying the nightside distribution of helium densities 
and NO nightglow intensities. A self-consistent model formulation (and mechanism) 
needs to be developed to explain Venus thermospheric superrotation. These initial 
sensitivity studies serve as useful starting points for future work. 
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6.6. Concept of homopause revisited 

One-dimensional coupled continuity-diffusion models are a reasonable first step in 
modeling the vertical density profiles observed in the Venus dayside thermosphere. 
Von Zahn et al. (1980) developed a morningside model at 60 SZA for examination of 
BNMS densities above 130 km. The success of the model fit to observations is in large 
part due to the chosen estimate of the altitude dependance of the eddy coefficient profile 
used for eddy diffusion. The analytical representation of the eddy coefficient profile was 
taken from Lindzen (1971) as K = AjJH, where N is total number density, and A is 
a tuning factor chosen to be 1.4 x 10 13 cm 0 5 s - The helium data taken by the BNMS 
were used to constrain the eddy profile and provide this specific factor. The upper level 
cutoff was set at K - 5 x 10* cm 2 s" \ despite the fact that unique values above the 
homopause (where eddy and molecular diffusion are equal) are not easily determined 
or meaningful. Homopause altitudes were established at 136, 130, and 134 km for N 2 , 
He, and mean composition, respectively. A similar approach was taken by Massie et al. 
(1983), who incorporated detailed chemistry and von Zahn et al. (1980) eddy diffusion 
coefficients (K£ 5 x 10 7 cm 2 s" ') in an improved dayside model extending down to 
100 km. A nightside mean model was also constructed assuming dayside fluxes of O- 
and N-atoms as upper boundary conditions for the nightside equations. One 
dimensional models of this type are subject to several limitations, the principal one of 
which concerns the neglected or improperly parameterized effects of large-scale winds 
on the calculated densities (Bougher, 1985). 

By contrast, Venus circulation model studies of Bougher et al. (1990a) and Mengel 
et al. (1989) suggest that vertical small-scale eddy diffusion is only a partial contributor 
to the maintenance of observed day and nightside thermospheric densities. Eddy 
coefficients required are about 2-3 times smaller than values used by previous one- 
dimensional composition models (von Zahn et al., 1980; Stewart et al., 1980; Massie 
et al., 1983). This implies that the vertical eddy coefficient dependence derived by 
von Zahn etal. (1980) within a one-dimensional model is not solely a signature of 
small-scale vertical mixing. Rather, it provides a reasonable description of compo- 
sitional variations brought about by a combination of large-scale horizontal and vertical 
winds and small-scale mixing. The conventional homopause, a well-defined level where 
molecular diffusion and eddy mixing processes balance, might better be described as a 
statistically average level where gravitational separation is offset by global wind plus 
small-scale mixing effects. The changing Venus thermospheric circulation has a modest 
influence on this level at any time. In general, one-dimensional eddy mixing models are 
valuable since they quantify the magnitude of the total dynamical effect expected; 
however, no information on the relative contribution of winds or turbulence can be 
obtained (von Zahn etal., 1980; Bougher et al., 1986, 1990a). 

Current Venus modeling efforts (Bougher et al., 1990a; Mengel et al., 1989) disagree 
as to the magnitude of this residual eddy diffusion required to supplement global 
circulation in maintaining observed vertical and diurnal density distributions at solar 
maximum. The bougher etal. (1989) VTGCM cannot tolerate global eddy diffusion 
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exceeding 2 x 10 7 cm 2 s~‘ and still yield atomic-N and O densities that reasonably 
match Hedin et ai (1983) empirical model values. Furthermore, the altitude of the NO 
nightglow peak volume emission rate (115 + 2 km ) is a very sensitive indicator of the 
nightside eddy diffusion possible within this global circulation model. Global winds for 
this code reach a maximum of — 100-200 ms -1 across the terminators for a super- 
position of SS-AS and prescribed zonal winds having a period of 8 days. Conversely, 
the Mengel et ai (1989) spectral model uses eddy diffusion reaching 3 x 10 7 cm 2 s “ \ 
in concert with maximum horizontal winds of — 100 m s ~ *, to obtain observed diurnal 
helium distributions. Superrotating winds are included with a 6-day period. 

The discrepancies in eddy, wind (wave-drag), and superrotation parameters in these 
two models reflect the difficulty in obtaining a unique solution without further observa- 
tions, namely measurements of thermospheric winds. The effectiveness of the circu- 
lation in building up the nightside densities of atomic-0 (and lighter species) decreases 
as the eddy coefficient increases. Likewise, the magnitude of the circulation regulates 
the amount that can be transported from day to night. Furthermore, superrotating winds 
not only shift the local time position of the nightside helium, hydrogen and NO airglow 
bulges, but they also serve to ‘smear out’ the magnitudes of these bulges. Therefore, 
future progress in Venus thermospheric modeling will likely require either remote or 
in situ wind measurements for convergence on a unique combination of wind, eddy, and 
superrotation parameters. 


7. Solar Cycle Variation 

Very little thermospheric data exists for Venus at solar minimum. A range of exospheric 
temperatures (275 + 50 K) was obtained from Mariner 5 and 10 La data (Anderson, 
1976; Takacs et ai , 1980) taken in October 1967 and February 1974, respectively. The 
latter corresponds to a near minimum in the solar cycle (F 10 7 = 68). The Venera 9 and 
10 orbiters also visited the planet in October 1975 through February 1976. Remote La 
(Bertaux et ai % 1978) and visible airglow (Krasnopolsky et ai y 1976) observations were 
made (F l07 = 75). No definitive estimate of neutral temperatures could, however, be 
derived. Lastly, the Pioneer Venus orbiter took in situ neutral data for nearly three 
diurnal periods (1978-1980). Periapsis altitude was actively controlled to remain within 
an altitude range of about 145-160 km until June, 1980. Afterward, this altitude was 
allowed to rise to a maximum mean value of about 2200 km in 1986. Presently, periapsis 
is moving downward toward an entry and bum-up in 1992. The only in situ neutrai/ion 
observations available are for near solar maximum conditions (F 107 = 150-240). 

Nevertheless, remote information regarding the ionosphere/thermosphere at solar 
minimum was obtained from various radio occultation measurements (Kliore and 
Mullen, 1988). Most of this data in fact comes from the Pioneer Venus Orbiter, although 
a few profiles come from Venera 9/10 and Mariner 10 (Kliore and Mullen, 1987; Kim 
etai , 1989). Photochemical equilibrium conditions hold for the major dayside ion 
densities below — 1 80 km (Cravens et ai , 1981 ; Kim et ai , 1989). Therefore, the most 
useful electron density measurements for inferring neutral temperatures are just above 
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the ion peak (£140 km). A simple model (Kliore and Mullen, 1987) suggests that 
neutral temperatures during 1984-1986 at 150 km (55-75 SZA) are about 215 K, nearly 
60 K lower than values at solar maximum (275 K). Mars derived neutral temperatures 
near its dayside ionospheric peak, also based on radio occultation electron density data, 
similarly show a rather weak dependence on F ]0 7 (Bauer and Hantsch, 1989). Solar 
cycle exospheric temperatures, however, vary by a larger amount than ionospheric peak 
values, especially for Mars (Bougher and Dickinson, 1988; Bougher eta!., 1988a). 

Bauer and Taylor (1981) sought to explain the weak solar cycle response of the Venus 
thermosphere compared to that of the Earth as arising from the stronger temperature 
dependence of the conductivity of C0 2 (K ~ T' 2i ) compared to that ofO (K ~ T° 7I ). 
Dickinson and Bougher (1986) showed, however, that conduction is the dominant 
cooling mechanism only above 160 km on Venus, where O is the dominant constituent. 

The Venus electron/ion density calculations of Kim et al. (1989) indicate that the 
Hedin et al. (1983) model predicts the neutral densities reasonably well for dayside solar 
minimum conditions; i.e., the resulting ion peak is given at the observed altitude 
(139-140 km). The empirical models of Hedin et al. (1983) and Keating etal. (1985), 
based primarily on in situ solar maximum data, predict a solar cycle exospheric tempera- 
ture variation (F l0 7 = 200 to 70) of about 55 to 70 K (310 to 255 K and 310 to 240 K, 
respectively). The well monitored 27-day exospheric temperature variation of ± 1 1 K 
(cf. Keating and Bougher, 1987) also suggests that a larger, non-negligible, solar cycle 
variation exists. However, without additional Venus in situ density (or temperature) 
measurements, a firm confirmation of this rather weak thermospheric response to the 
solar cycle is still lacking! Figure 48 illustrates a collection of possible dayside Venus 
solar minimum neutral temperature profiles. The only data on the plot comes from the 
PV-ORO radio occultation measurements of Kliore and Mullen (1987). 

The true test of a useful predictive model is its ability to reliably simulate the observed 
fields under a wide range of conditions. The NCAR 1-D (Dickinson and Bougher, 

1986) , 2-D (Bougher et al., 1986) and VTGCM (Bougher et al., 1990a) models have been 
modified to predict the Venus solar minimum thermospheric circulation (dynamical 
models) and structure. The primary solar cycle change occurs in the EUV fluxes 
shortward of 105 nm; i.e., a reduction of peak EUV heating and net C0 2 dissociation 
is given of approximately a factor of 3. This is consistent with a 21% smaller peak ion 
density (55-75 SZA) for solar minimum than solar maximum (Kliore and Mullen, 

1987) . Indirectly, through cooler temperatures, C0 2 1 5-pm cooling is also reduced, since 
it is strongly nonlinearly temperature dependent. The EUV fluxes of Heroux and 
Hinteregger (1978) were chosen, corresponding to anf l01 = 74. 

Dickinson and Bougher (1986) confirmed that the Venus exospheric temperatures 
should be relatively insensitive to changes in the solar flux; dayside values were shown 
to change by less than ~ 65 K over the solar cycle (300 to 235 K). The strong nonlinear 
temperature dependence of the C0 2 cooling and the variation of atomic oxygen densities 
with solar activity (‘the atomic O thermostat’) effectively buffer against solar pertur- 
bations, thereby reducing the temperature variation from that which would occur if a 
balance were present between peak EUV heating and molecular conduction. Corre- 
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VENUS SMIN NEUTRAL TEMPERATURES (x~60°SZA) 



Fig. 48. Venus dayside solar minimum neutral temperature profiles (60 SZA). VTS3 and VIRA profiles 
are only reasonable extrapolations from near solar maximum data. The one true data point is inferred from 
ORO measured electron density profiles just above the ion peak (150 km). The NCAR 1-D (V1D) and 
VTGCM model predictions are also shown for comparison. 


sponding nightside mean calculations suggested that little change in temperatures with 
solar activity should be expected, so that global mean values drop by only 30 K from 
solar maximum. 

Hydrodynamic effects were carefully examined in subsequent Venus circulation 
models for an improved description of solar cycle response. Figures 4 l(e-f) illustrate 
NCAR 2-D solar minimum contour plots, for comparison to corresponding solar maxi- 
mum fields of Figures 41(a) and 4 1(b). The dayside exospheric temperatures are reduced 
from 300 K to — 230 K, giving a dayside solar cycle temperature variation of at most 
70 K. These newly calculated dayside solar minimum temperatures are somewhat low, 
yet the values fall within the range of exospheric temperatures estimated from Mariner 5 
and 10 Lx data (Anderson, 1976; Takacs etai , 1980). Also, values at 150 km (over 
55-75 SZA) are close to 215 K, quite similar to those inferred from 1984-1986 radio 
occultation data (Kliore and Mullen, 1987). Adiabatic cooling serves to reduce dayside 
temperatures slightly below previous radiative equilibrium model values (Dickinson and 
Bougher, 1986). The solar driven winds are also 20-30 m s “ 1 weaker; this slight change 
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causes the nightside thermospheric major densities and temperatures to remain largely 
the same. Both the NCAR 2-D and VTGCM solar minimum models suggest that the 
Venus thermospheric major densities and temperatures are relatively insensitive to the 
factor of 3 reduction in EUV inputs over the solar cycle (Bougher et ai , 1986, 1990a). 

This does not appear to be the case for tracer species and the UV nightglow. 
Odd-nitrogen studies using the VTGCM (see Section 3.10) predict a variation of nearly 
a factor of 3 in the NO (0, 1) 6-band intensities over the solar cycle (see Figures 20(a) 
and 20(b)). This sensitivity is due to the net change in dayside production of N-atoms, 
which varies strongly with solar cycle. A detailed analysis of recently obtained 
PV-OUVS data taken from 1984-1986 will provide an opportunity to validate this 
VTGCM solar cycle nightglow variation. A reasonable match of the model and OUVS 
imaged dark-disk intensity variation over the solar cycle would be a substantial 
validation of the odd-nitrogen chemistry adopted and the solar cycle response predicted 
by the VTGCM thermospheric winds. Also, since the NO nightglow is a good tracer 
of the Venus thermospheric circulation, its monitoring over the solar cycle should 
provide a remote means to study the variation of the winds once the dayside N-produc- 
tion is calculated. 
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